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Preface 



The area: Chaos in Partial Differential Equations, is at its fast developing 
stage. Notable results have been obtained in recent years. The present book aims 
at an overall survey on the existing results. On the other hand, we shall try to 
make the presentations introductory, so that beginners can benefit more from the 
book. 

It is well-known that the theory of chaos in finite-dimensional dynamical sys- 
tems has been well-developed. That includes both discrete maps and systems of 
ordinary differential equations. Such a theory has produced important mathe- 
matical theorems and led to important applications in physics, chemistry, biology, 
and engineering etc.. For a long period of time, there was no theory on chaos in 
partial differential equations. On the other hand, the demand for such a theory 
is much stronger than for finite-dimensional systems. Mathematically, studies on 
infinite-dimensional systems pose much more challenging problems. For example, as 
phase spaces, Banach spaces possess much more structures than Euclidean spaces. 
In terms of applications, most of important natural phenomena are described by 
partial differential equations - nonlinear wave equations. Maxwell equations, Yang- 
Mills equations, and Navier-Stokes equations, to name a few. Recently, the author 
and collaborators have established a systematic theory on chaos in nonlinear wave 
equations. 

Nonlinear wave equations are the most important class of equations in natu- 
ral sciences. They describe a wide spectrum of phenomena - motion of plasma, 
nonlinear optics (laser), water waves, vortex motion, to name a few. Among 
these nonlinear wave equations, there is a class of equations called soliton equa- 
tions. This class of equations describes a variety of phenomena. In particular, 
the same soliton equation describes several different phenomena. Mathematical 
theories on soliton equations have been well developed. Their Cauchy problems 
are completely solved through inverse scattering transforms. Soliton equations are 
integrable Hamiltonian partial differential equations which are the natural coun- 
terparts of finite-dimensional integrable Hamiltonian systems. We have established 
a standard program for proving the existence of chaos in perturbed soliton equa- 
tions, with the machineries: 1. Darboux transformations for soliton equations, 2. 
isospectral theory for soliton equations under periodic boundary condition, 3. per- 
sistence of invariant manifolds and Fenichel fibers, 4. Melnikov analysis, 5. Smale 
horseshoes and symbolic dynamics, 6. shadowing lemma and symbolic dynamics. 

The most important implication of the theory on chaos in partial differential 
equations in theoretical physics will be on the study of turbulence. For that goal, 
we chose the 2D Navier-Stokes equations under periodic boundary conditions to 
begin a dynamical system study on 2D turbulence. Since they possess Lax pair 
and Darboux transformation, the 2D Euler equations are the starting point for an 
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analytical study. The high Reynolds number 2D Navier-Stokes equations are viewed 
as a singular perturbation of the 2D Euler equations through the perturbation 
parameter e = 1/Re which is the inverse of the Reynolds number. 

Our focus will be on nonlinear wave equations. New results on shadowing 
lemma and novel results related to Euler equations of inviscid fluids will also be 
presented. The chapters on figure-eight structures and Melnikov vectors are written 
in great details. The readers can learn these machineries without resorting to other 
references. In other chapters, details of proofs are often omitted. Chapters 3 
to 7 illustrate how to prove the existence of chaos in perturbed soliton equations. 
Chapter 9 contains the most recent results on Lax pair structures of Euler equations 
of inviscid fluids. In chapter 12, we give brief comments on other related topics. 

The monograph will be of interest to researchers in mathematics, physics, engi- 
neering, chemistry, biology, and science in general. Researchers who are interested 
in chaos in high dimensions, will find the book of particularly valuable. The book is 
also accessible to graduate students, and can be taken as a textbook for advanced 
graduate courses. 

I started writing this book in 1997 when I was at MIT. This project continued 
at Institute for Advanced Study during the year 1998-1999, and at University of 
Missouri - Columbia since 1999. In the Fall of 2001, I started to rewrite from the 
old manuscript. Most of the work was done in the summer of 2002. The work was 
partially supported by an AMS centennial fellowship in 1998, and a Guggenheim 
fellowship in 1999. 

Finally, I would like to thank my wife Sherry and my son Brandon for their 
strong support and appreciation. 



CHAPTER 1 



General Setup and Concepts 



We are mainly concerned with the Cauchy problems of partial differential equa- 
tions, and view them as defining flows in certain Banach spaces. Unlike the Eu- 
cli(k;an space M", such Banach spaces admit a variety of norms which make the 
structures in infinite dimensional dynamical systems more abundant. The main 
difficulty in studying infinite dimensional dynamical systems often comes from the 
fact that the evolution operators for the partial differential equations are usually 
at best C° in time, in contrast to finite dimensional dynamical systems where the 
evolution operators are smooth in time. The well-known concepts for finite di- 
mensional dynamical systems can be generalized to infinite dimensional dynamical 
systems, and this is the main task of this chapter. 



The types of evolution equations studied in this book can be casted into the 
general form. 



where t G (time), x — (xi, . . . , a;„) € R", Q = {Qi, ■ ■ ■ ,Qm,) and G = 
{Gi, . . . ,Gm) are either real or complex valued functions, and £, m and n are 
integers. The equation (1.1) is studied under certain boundary conditions, for ex- 
ample, 

• periodic boundary conditions, e.g. Q is periodic in each component of x 

with period 27r, 

• decay boimdary conditions, e.g. Q ^ as x ^ oo. 

Thus we have Cauchy problems for the equation (1.1), and we would like to pose 
the Cauchy problems in some Banach spaces H, for example, 

• H can be a Sobolev space H'', 

• H can be a Solobev space Hj^p of even periodic functions. 
We require that the problem is well-posed in 7i, for example, 

• for any Qq G H, there exists a unique solution Q = Q{t, Qo) G C°[(— oo, oo); H] 
or C°[[0;oo),W] to the equation (1.1) such that Q{0,Qo) = Qo, 

• for any fixed to G (—00,00) or [0, 00), Q(toj<3o) is a C function of Qo, 
for Qo gH and some integer r > 0. 

ExEimple: Consider the integrable cubic nonlinear Schrodinger (NLS) equa- 
tion. 



1.1. Cauchy Problems of Partial Differential Equations 



(1.1) 



dtQ = G{Q,d,Q,...,dlQ) 



(1.2) 



iqt = + 2 [I g 1^ -w^] q, 
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1. GENERAL SETUP AND CONCEPTS 



where i = t G M^, a; G K^, g is a complex- valued function of {t, x), and w is a 

real constant. We pose the periodic boundary condition, 

q{t,x + l) = q{t,x) . 

The Cauchy problem for equation (1.2) is posed in the Sobolev space of periodic 
functions, 

H = |g = {q, q) q{x + 1) = q{x), q G HI^^^^ : the 

Sobolev space over the period interval [0, 1] | , 

and is well-posed [38] [32] [33]. 

Fact 1: For any Qq g Ti, there exists a unique solution Q = Q{t,Qo) G 

C'^[(— oo, oo), 7i] to the equation (1.2) such that Q(0, Qo) = Qo- 

Fact 2: For any fixed to G (—00,00), Q{to,Qo) is a function of Qo, for 



1.2. Phase Spaces and Flows 

For finite dimensional dynamical systems, the phase spaces are often M" or C". 
For infinite dimensional dynamical systems, we take the Banach space H discussed 
in the previous section as the counterpart. 

Definition 1.1. We call the Banach space H in which the Cauchy problem 
for (1.1) is well-posed, a phase space. Define an operator labeled by t as 

Q{t,Qo) = F\Qo); 

then F* : H ^ H is called the evolution operator (or flow) for the system (1.1). 

A point p G H is called a fixed point if F* (p) = p for any t. Notice that here 
the fixed point p is in fact a function of x, which is the so-called stationary solution 
of (1.1). Let q £ H he a, point; then £g = {F*(g),for all t} is called the orbit with 
initial point q. An orbit £q is called a periodic orbit if there exists a T G (—00,00) 
such that F^ (q) = q. An orbit Ig is called a homoclinic orbit if there exists a point 
G H such that F*{q) — » as | t |— » 00, and g« is called the asymptotic point 
of the homoclinic orbit. An orbit £q is called a heteroclinic orbit if there exist two 
different points q± G H such that F*{q) q±, as t — > ±00, and q± are called the 
asymptotic points of the heteroclinic orbit. An orbit £q is said to be homoclinic to 
a submanifold W oiH if infg^w \\ F*{q) — Q \\^ 0, as j t j— > 00. 

Example 1: Consider the same Cauchy problem for the system (1.2). The 
fixed points of (1.2) satisfy the second order ordinary differential equation 

(1.3) fc + 2[|g|2-a;2]g = 0. 

In particular, there exists a circle of fixed points q = we*''', where 7 G [0, 2tt]. For 
simple periodic solutions, we have 



(1.4) 



q = ae'^W^ e{t) = - [2{a^ - uj'^)t - 7] ; 
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where a > 0, and 7 G [0, 27r]. For orbits homoclinic to the circles (1.4), we liave 
cos 2p — sinp sech r cos 2'kx — i sin 2p tanh r 



(1-5) Q = T 



A = 1 + sinp sech r cos 27rx , 

, p G (—00, 00) is the Backlund 



where r = 47r\/ — tt^ i + p, p = arctan 



parameter. Setting a = w in (1.5), we have heteroclinic orbits asymptotic to points 
on the circle of fixed points. The expression (1.5) is generated from (1.4) through 
a Backlund-Darboux transformation [137]. 

Example 2: Consider the sine-Gordon equation, 

utt - + sin u = , 

under the decay boundary condition that u belongs to the Schwartz class in x. The 
well-known "breather" solution, 

, . tani/cosf(cosj/)tl 

(1.6) w(t,.T) = 4 arctan 



cosh[(sin v)x\ 

where is a parameter, is a periodic orbit. The expression (1.6) is generated from 
trivial solutions through a Backlund-Darboux transformation [59] . 

1.3. Invariant Submanifolds 

Invariant submanifolds are the main objects in studying phase spaces. In phase 
spaces for partial differential equations, invariant submanifolds are often subman- 
ifolds with boundaries. Therefore, the following concepts on invariance are impor- 
tant. 

Definition 1.2 (Overflowing and Inflowing Invariance). A submanifold W 
with boundary dW is 

• overflowing invariant if for any t > 0, W G o W, where = U dW, 

• inflowing invariant if any t > 0, F*^ o W C W, 

• invariant if for any t > 0, F* oW = W. 

Definition 1.3 (Local Invariance). A submanifold W with boundary dW is 

locally invariant if for any point q G W, if (J F*{q) W, then there exists T G 

te[o,(x)) 

(0, 00) such that (J F*{q) C W, and F'^{q) G dW; and if [j F*{q)g^W, 

te[0,T) t6(-oo,0] 

then there exists T e (-00, 0) such that (J F\q) c W, and F'^{q) e dW. 

*G(T,0] 

Intuitively speaking, a submanifold with boundary is locally invariant if any 
orbit starting from a point inside the submanifold can only leave the submanifold 
through its boundary in both forward and backward time. 

ExEimple: Consider the linear equation, 

(1-7) iqt = (1 + i)qxx + iq , 

where i = a/^, f G M^, a; G R^, and g is a complex- valued function of (t, x), under 
periodic boundary condition, 

q{x + l)= q{x) . 
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1. GENERAL SETUP AND CONCEPTS 



Let q — e*'^j*+*'^3^; then 

where ^ 2j7r, (j S Z). f^o = 1, and when | j |> 0, Re{^j\ < 0. Wc take the 
space of periodic functions of period 1 to be the phase space. Then the submanifold 

Wo = {q & q = cq, cq is complex and || g ||< 1 

is an outflowing invariant submanifold, the submanifold 

Wi 

is an inflowing invariant submanifold, and the submanifold 

q = Cq + Cie*'^^^, Co andci are complex, and || q ||< 1 
is a locally invariant submanifold. The unstable subspace is given by 

q = Co, Co is complex 



riant sub 

|g e -ff^ q = cie*'"'^, ci is complex, and || g ||< l| 



W=!^qGH^ 



and the stable subspace is given by 



q= CjC^'^^^, Cj 's are complex 

i62/{0} 

Actually, a good way to view the partial differential equation (1.7) as defining an 
infinite dimensional dynamical system is through Fourier transform, let 

q{t,x) = Y,c,{t)e"'^-; 
jez 

then Cj{t) satisfy 

cj=[{l-k^)+ik]]cj, j€Z; 
which is a system of infinitely many ordinary differential equations. 

1.4. Poincare Sections and Poincare Maps 

In the infinite dimensional phase space 7Y, Poincare sections can be defined in 
a similar fashion as in a finite dimensional phase space. Let Ig be a periodic or 
homoclinic orbit in H under a flow F*, and g* be a point on Ig, then the Poincare 
section S can be defined to be any codimension 1 subspace which has a transversal 
intersection with Ig at g*. Then the flow will induce a Poincare map P in the 
neighborhood of g* in Sq. Phase blocks, e.g. Smale horseshoes, can be defined 
using the norm. 



CHAPTER 2 



Soliton Equations as Integrable Hamiltonian PDEs 



2.1. A Brief Summary 

Soliton equations are integrable Hamiltonian partial differential equations. For 
example, the Korteweg-de Vries (KdV) equation 



where u is a real-valued function of two variables t and x, can be rewritten in the 
Hamiltonian form 

ou 

where 

'1 



H = 



2^1 - 



dx 



under either periodic or decay boundary conditions. It is integrable in the classi- 
cal Liouville sense, i.e., there exist enough functionally independent constants of 
motion. These constants of motion can be generated through isospcctral theory 
or Backlund transformations [8]. The level sets of these constants of motion are 
elliptic tori [178] [154] [153] [68]. 

There exist soliton equations which possess level sets which are normally hy- 
perbolic, for example, the focusing cubic nonlinear Schrodinger equation [137], 



iqt = Qxx + 2|g 



2, 



where i = v— T and g is a complex- valued function of two variables t and x; the 
sine-Gordon equation [157], 

Utt = Uxx -I- sin u , 

where m is a real-valued function of two variables t and x, etc. 

Hyperbolic foliations are very important since they are the sources of chaos 
when the integrable systems are under perturbations. We will investigate the hy- 
perbolic foliations of three typical types of soliton equations: (i). (l+l)-dimensional 
soliton equations represented by the focusing cubic nonlinear Schrodinger equation, 

(ii) . soliton lattices represented by the focusing cubic nonlinear Schrodinger lattice, 

(iii) . (l-}-2)-dimensional soliton equations represented by the Davey-Stewartson 11 
equation. 

Remark 2.1. For those soliton equations which have only elliptic level sets, 

the corresponding representatives can be chosen to be the KdV equation for (1+1)- 
dimensional soliton equations, the Toda lattice for soliton lattices, and the KP 
equation for (l-|-2)-dimensional soliton equations. 
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2. SOLITON EQUATIONS AS INTEGRABLE HAMILTONIAN PDES 



Soliton equations are canonical equations which model a variety of physical 
phenomena, for example, nonlinear wave motions, nonlinear optics, plasmas, vortex 
dynamics, etc. [5] [1]. Other typical examples of such integrable Hamiltonian 
partial differential equations are, e.g., the defocusing cubic nonlinear Schrodinger 
equation, 

iqt = ixx - 2|g|^g , 

where i = ^/~^ and g is a complex- valued function of two variables t and x; the 
modified KdV equation, 

Ut = ±6u'^Ux - Uxxx , 

where m is a real- valued function of two variables t and x; the sinh-Gordon equation, 

Utt = Uxx + sinh u , 

where u is a real- valued function of two variables t and x; the three- wave interaction 
equations, 

duj _^ ^ duj 
dt * dx 

where fc = 1,2,3 are cyclically permuted, aj and 6j are real constants, Ui are 
complex-valued functions of t and x; the Boussinesq equation, 

where u is a real-valued function of two variables t and x; the Toda lattice, 

d'^Un/dt^ = exp{-(u„ - Un-i)} - exp{-(u„+i - u„)} , 
where u„'s are real variables; the focusing cubic nonlinear Schrodinger lattice, 

i-^ = {qn+i - 2qn + Qn-i) + |g„P(g„+i + , 
where q^s are complex variables; the Kadomtsev-Petviashvili (KP) equation, 

{Ut + 6uUx + Uxxx)x = ±3Uj,j/ , 

where u is a real-valued function of three variables t, x and y; the Davey-Stewartson 
II equation, 

id,q=[dl-dl]q+[2\q\-'+Uy]q, 

[dl + dl]u = -Ady\q\' , 

where i = \/— T, g' is a complex-valued function of three variables t, x and y; and 
M is a real- valued function of three variables t, x and y. For more complete list of 
soliton equations, see e.g. [5] [1]. 

The cubic nonlinear Schrodinger equation is one of our main focuses in this 
book, which can be written in the Hamiltonian form, 

5H 
Sq 

where 

H= [[-\qx\^±\q\yx , 



2.2. A PHYSICAL APPLICATION OF THE NONLINEAR SCHRODINGER EQUATION 7 



under periodic boundary conditions. Its phase space is defined as 



q = 



r= -q, q{x + l) = q{x), 

q e ^^[0,1] • Sobolev space H'' over the period interval [0, 1]| . 

Remark 2.2. It is interesting to notice that the cubic nonlinear Schrodinger 
equation can also be written in Hamiltonian form in spatial variable, i.e., 

Qxx = iqi ± 2\qfq , 
can be written in Hamiltonian form. Let p = qx', then 



d_ 

dx 



( q\ 

q 
P 

\P J 



/ ML \ 

I Sq \ 



5H 
Sq 



SH 

5p 



\ Sp / 



where 



J = 



/ 








1 





\ 













1 






-1 













V 





-1 








/ 



H = j[\p?T\q\^-\{qtq-qtq)]dt , 



under decay or periodic boundary conditions. We do not know whether or not other 
soliton equations have this property. 



2.2. A Physical Application of the Nonlinear Schrodinger Equation 

The cubic nonlinear Schrodinger (NLS) equation has many different applica- 
tions, i.e. it describes many different physical phenomena, and that is why it is 
called a canonical equation. Here, as an example, we show how the NLS equation 
describes the motion of a vortex filament - the beautiful Hasimoto derivation [82] . 
Vortex filaments in an inviscid fluid are known to preserve their identities. The 
motion of a very thin isolated vortex filament X = X{s, t) of radius e in an incom- 
pressible inviscid unbounded fluid by its own induction is described asymptotically 

by 

(2.1) dX/dt^GKb, 

where s is the length measured along the filament, t is the time, k is the curvature, 
b is the unit vector in the direction of the binormal and G is the coefficient of local 
induction, 

G=^[ln(l/e) + 0(l)] , 
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2. SOLITON EQUATIONS AS INTEGRABLE HAMILTONIAN PDES 



which is proportional to the circulation F of the filament and may be regarded as 
a constant if we neglect the second order term. Then a suitable choice of the units 
of time and length reduces (2.1) to the nondimensional form, 

(2.2) dX/dt = Kb. 

Equation (2.2) should be supplemented by the equations of differential geometry 

(the Frcnet-Scrct formulae) 

(2.3) dX/ds ~ t , dt/ds — ku , dn/ds = rb — nt , db/ds = —rn , 

where r is the torsion and t, n and b are the tangent, the principal normal and the 
binormal unit vectors. The last two equations imply that 

(2.4) d{n + ib)/ds = -iT{n + ib) - Kt , 
which suggests the introduction of new variables 

(2.5) N = {n + ib) exp ^^i J rds 
and 

(2.6) g = Kexp'|i J rds 
Then from (2.3) and (2.4), we have 

(2.7) dN/ds = -qt, dt/ds = Re{qN} = ^{qN + qN) . 

We are going to use the relation = to derive an equation for q. For this 

we need to know dt/dt and dN/dt besides equations (2.7). From (2.2) and (2.3), 
we have 

d^X 

dt/dt = ~ d{nb)/ds = {dn/ds)b— ktH 

= K Re{{— dn/ds + iT){b + in)} , 

K 

i.e. 

(2.8) dt/dt = Re{i{dq/ds)N} = ^i[{dq/ds)N - {dq/ds)' N] . 
We can write the equation for dN/dt in the following form: 

(2.9) dN/dt = aN + + -ft, 
where a, (3 and 7 are complex coefficients to be determined. 

a + a = ^[dN/dt-N + dN/dt-N] 

= ^d{N-^)/dt = , 
i.e. a = iR where R is an unknown real function. 

13 = ^dN/dt ■ N = ^d{N ■ N)/dt = , 
7 = -iV • dt/dt = -idq/ds . 
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Thus 

(2.10) dN/dt = i[RN - {dq/ds)t] . 

From (2.7), (2.10) and (2.8), we have 



dsdt 



-{dq/dt)t- qdt/dt 
-{dq/dt)t- ^iq[{dq/ds)N - [dq/ds)- N] , 



52 A/- 
dtds 



= i[{dR/ds)N - Rqt- {d'^q/ds^)t 
-^{dq/ds){qN + q'^)] . 

Thus, we have 

(2.11) dq/dt = i[d'^q/ds^ + Rq] , 
and 

(2.12) ^qdq/ds = dR/ds - ^{dq/ds)q . 

The comparison of expressions for from (2.7) and (2.8) leads only to (2.11). 
Solving (2.12), we have 

(2.13) R=^Qq\^+A), 

where A is a real-valued function of t only. Thus we have the cubic nonlinear 
Schrodinger equation for q: 

-idq/dt = d'^q/ds^ + i(|gp + A)q . 
The term Aq can be transformed away by defining the new variable 



1 /■* 

g = gexp[--i A{t)dt\ . 



CHAPTER 3 



Figure-Eight Structures 



For finite-dimensional Hamiltonian systems, figure-eight structures are often 
given by singular level sets. These singular level sets are also called separatrices. 
Expressions for such figuro-cight structures can be obtained by setting the Hamil- 
tonian and/or other constants of motion to special values. For partial differential 
equations, such an approach is not feasible. For soliton equations, expressions 
for figure-eight structures can be obtained via Backlund-Darboux transformations 
[137] [118] [121]. 

3.1. ID Cubic Nonlinear Schrodinger (NLS) Equation 

We take the focusing nonlinear Schrodinger equation (NLS) as our first example 
to show how to construct figure-eight structures. If one starts from the conserva- 
tional laws of the NLS, it turns out that it is very elusive to get the separatrices. On 
the contrary, starting from the Backlund-Darboux transformation to be presented, 
one can find the separatrices rather easily. We consider the NLS 

(3.1) iqt = Qxx + 2|(?p(? , 

under periodic boundary condition q{x + 2tt) = q{x). The NLS is an integrable 
system by virtue of the Lax pair [212], 

(3.2) = Uip , 

(3.3) = V<fi, 



where 



U = iA(T3 -|- i 



V = 2i A^CTs -|- iqras + 



q 
-r 



2i\q 



-2iXr + 

where denotes the third Pauli matrix = diag(l,— 1), r = —q, and A is the 
spectral parameter. If q satisfies the NLS, then the compatibility of the over deter- 
mined system (3.2, 3.3) is guaranteed. Let M = M(x) be the fundamental matrix 
solution to the ODE (3.2), M(0) is the 2x2 identity matrix. We introduce the 
so-called transfer matrix T = T(A, q) where q = {q, —q), T = M(27r). 

Lemma 3.1. Let Y{x) be any solution to the ODE (3.2), then 

Y{2n-K) = T" y(0) . 
Proof: Since M{x) is the fundamental matrix, 

Y{x) = M{x) Y{0) . 
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Thus, 

y(27r) = T Y{0) . 

Assume that 

Y{2lTr) = T' Y{0) . 
Notice that Y{x + 2ln) also solves the ODE (3.2); then 

Y{x + 21tx) = M{x) Y(21tt) ; 

thus, 

Y{2{1 + l)7r) = T Y{21it) = T^+^ Y{Q) . 
The lemma is proved. Q.E.D. 

Definition 3.2. We define the Floquet discriminant A as, 
A{X,q) = trace {T{X,q)} . 
We define the periodic and anti-periodic points A^^") by the condition 

|A(A(P),g)|=2. 

We define the critical points A*^^) by the condition 

dA{X,q) 



= 

A=A(=) 



dx 

A multiple point, denoted A^'"^ is a critical point for which 

|A(A('"),g)|=2. 

The algebraic multiplicity of A^"*' is defined as the order of the zero of A(A) ± 2. 
Ususally it is 2, but it can exceed 2; when it docs equal 2, wc call the multiple point 
a double point, and denote it by A^''). The geometric multiplicity of A^™) is defined 
as the maximum number of linearly independent solutions to the ODE (3.2), and 
is either 1 or 2. 

Let q{x.t) bo a solution to the NLS (3.1) for which the linear system (3.2) 
has a complex double point u of geometric multiplicity 2. We denote two linearly 
independent solutions of the Lax pair (3.2,3.3) at A = z/ by (0+, </>"). Thus, a 
general solution of the linear systems at {q, u) is given by 

(3.4) (j){x,t) = c+(l)+ + c-(l)~ . 
We use (j) to define a Gauge matrix [180] G by 

(3.5) G = G{X;v;cI^)=n(^^-'' ^l_^N-^, 
where 

(3.6) iV=^ 
Then we define Q and 4" by 



"2 91 



(3.7) Q{x, t) = q{x, t) + 2{v - v) /^^ 
and 

(3.8) ^{x,t\X) = G{X;u;(l)) tp{x,t;X) 
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where ip solves the Lax pair (3.2,3.3) at {q,v). Formulas (3.7) and (3.8) are the 
Backlund-DarboTix transformations for the potential and eigenfunctions, respec- 
tively. We have the following [180] [137], 

Theorem 3.3. Let q{x,t) be a solution to the NLS equation (3.1), for which 
the linear system (3.2) has a complex double point v of geometric multiplicity 2, 
with eigenbasis for the Lax pair (3.2,3.3), and define Q{x,t) and ^{x,t; A) 

by (3.7) and (3.8). Then 

(1) Q{x,t) is an solution of NLS, with spatial period 2Tr, 

(2) Q and q have the same Floquet spectrum, 

(3) Q{x,t) is homoclinic to q{x,t) in the sense that Q{x,t) — > qg^{x,t), 
expontentially as eyip{—a^\t\), as t — > ±oo, where qe^ is a "torus trans- 
late" of q,(Jv is the nonvanishing growth rate associated to the complex 
double point u, and explicit formulas exist for this growth rate and for the 
translation parameters 0±, 

(4) *(a;,t;A) solves the Lax pair (3.2,3.3) at (Q,A). 

This theorem is quite general, constructing homoclinic solutions from a wide 
class of starting solutions q{x, t). It's proof is one of direct verification [118]. 

We emphasize several qualitative features of these homoclinic orbits: (i) Q{x,t) 
is homoclinic to a torus which itself possesses rather complicated spatial and tem- 
poral structure, and is not just a fixed point, (ii) Nevertheless, the homoclinic orbit 
typically has still more complicated spatial structure than its "target torus" . (iii) 
When there are several complex double points, each with nonvanishing growth rate, 
one can iterate the Backlund-Darboux transformations to generate more compli- 
cated homoclinic orbits, (iv) The number of complex double points with nonvanish- 
ing growth rates counts the dimension of the unstable manifold of the critical torus 
in that two unstable directions are coordinatized by the complex ratio c+/c_ . Un- 
der even symmetry only one real dimension satisfies the constraint of evenness, as 
will be clearly illustrated in the following example, (v) These Backlund-Darboux 
formulas provide global expressions for the stable and unstable manifolds of the 
critical tori, which represent figure-eight structures. 

Example: As a concrete example, we take q{x,t) to be the special solution 

(3.9) qc = cexp {-i[2cH + -f]} . 
Solutions of the Lax pair (3.2,3.3) can be computed explicitly: 

(3.10) </.(±)(a;,i;A)=e±-(-+2^*)f ^ ce-'(2^^*+T)/2 \ 



where 



K = k(A) = Vc2+A2 . 

With these solutions one can construct the fundamental matrix 



(3.11) M{x;X;qc) 



cos KX + iz: sin kx 



cos KX ^ i — sin kx 



from which the Floquet discriminant can be computed: 
(3.12) A(A; qc) = 2 cos(2«;7r) . 

From A, spectral quantities can be computed: 
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(1) simple periodic points: — ±i c , 

(2) double points: K{xf^) =j/2 

(3) critical points: A^'^' = A^**' , j G 

(4) simple periodic points: Aq'^^ = . 

For this spectral data, there are 2N purely imaginary double points, 



(3.13) 
where 



{Xfr=f/A-C^ i = l,2,.. 



< < 



(iV+ 1)2/4 -c^ 



From this spectral data, the homoclinic orbits can be explicitly computed through 
Backlund-Darboux transformation. Notice that to have temporal growth (and de- 
cay) in the eigenfunctions (3.10), one needs A to be complex. Notice also that 

the Backlund-Darboux transformation is built with quadratic products in <p, thus 
choosing v ~ X^^^ will guarantee periodicity of Q in x. When iV = 1, the 

Backlund-Darboux transformation at one purely imaginary double point A^'''' yields 
Q = Q{x,t;c,r,c+/c_) [137]: 



Q 



cos2p— sinp sechr cos(x -I- ■!?) — i sin2ptanhT 



(3.14) 



1 + sinp sechr cos(a; + 1?) 



gT2»Pcg-i(2c^t+7) as p 



ce 



TOO, 



where c^jc- = exp(p -I- ifi) and p is defined by 1/2 -|- i\J(? — 1/4 = cexp(zp), 
T = at-p, and '&=p - {(3 + 7r/2). 

Several points about this homoclinic orbit need to be made: 

(1) The orbit depends only upon the ratio c+/c_, and not upon c+ and c_ 
individually. 

(2) Q is homoclinic to the plane wave orbit; however, a phase shift of — 4p 
occurs when one compares the asymptotic behavior of the orbit as t ^ 
— DO with its behavior as t ^ + oo. 

(3) For small p, the formula for Q becomes more transparent: 



Q 



(cos 2p — i sin 2p tanh r) — 2 sin p sech r cos(x + i?) 



ce 



(4) An evenness constraint on Q in a; can be enforced by restricting the phase 
6 to be one of two values 



0,7r. 



(evenness) 



In this manner, the even symmetry disconnects the level set. Each com- 
ponent constitutes one loop of the figure eight. While the target q is in- 
dependent of X, each of these loops has x dependence through the cos(a;). 
One loop has exactly this dependence and can be interpreted as a spatial 
excitation located near x = 0, while the second loop has the dependence 
cos(a; — it), which we interpret as spatial structure located near x = tt. 
In this example, the disconnected nature of the level set is clearly related 
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to distinct spatial structures on the individual loops. See Figure 1 for an 
illustration. 

(5) Direct calculation shows that the transformation matrix M (1; X^^"^ ; Q) is 
similar to a Jordan form when t G {—oo, oo), 



M(1;A(''';Q) 



-1 1 
-1 



and when t ~f ±oo, A/(l; A^''''; Q) — > — / (the negative of the 2x2 identity 
matrix). Thus, when t is finite, the algebraic multiplicity {— 2) of A = A^^'*-' 
with the potential Q is greater than the geometric multiplicity {— 1). 

In this example the dimension of the loops need not be one, but is determined 
by the number of purely imaginary double points which in turn is controlled by the 
amplitude c of the plane wave target and by the spatial period. (The dimension 
of the loops increases linearly with the spatial period.) When there are several 
complex double points, Backlund-Darboux transformations must be iterated to 
produce complete representations. Thus, Backlund-Darboux transformations give 
global representations of the figure-eight structures. 

3.1.1. Linear Instability. The above figure-eight structure corresponds to 
the following linear instability of Benjamin- Feir type. Consider the uniform solution 
to the NLS (3.1), 

Let 

and linearize equation (3.1) at qc, we have 

iqt = qxx + 2c^ [q + q] ■ 
Assume that q takes the form. 



q 



cos kn 



where kj — 2jTT, {j = 0, 1, 2, • • •), Aj and Bj are complex constants. Then, 



=±A:, J4c2-fc2 



Thus, we have instabilities when c > 1/2. 
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3.1.2. Quadratic Products of Eigenfunctions. Quadratic products of eigen- 
functions play a crucial role in characterizing the hyperbolic structures of soliton 
equations. Its importance lies in the following aspects: (i). Certain quadratic 
products of eigenfunctions solve the linearized soliton equation, (ii). Thus, they 
are the perfect candidates for building a basis to the invariant linear subbundles. 
(iii). Also, they signify the instability of the soliton equation, (iv). Most impor- 
tantly, quadratic products of eigenfunctions can serve as Melnikov vectors, e.g., for 
Davey-Stewartson equation [121]. 

Consider the linearized NLS equation at any solution q{t,x) written in the 
vector form: 

idt{5q) = iSq),, + 2[q^S^+2\q\''5q] , 
(3.15) _ _ _ 

idtiSq) = -{5q)^^ - 2[f6q + 2\q\''Sq] , 

we have the following lemma [137]. 

Lemma 3.4. Let tp^^^ — ip^j^t,x; X,q) (j — 1,2) be any two eigenfunctions 
solving the Lax pair (3.2,3.3) at an arbitrary X. Then 



Sq 



(1) (2) 
'Pi Vl 



(1) (2) 



and S 



(1) (2) 

f i Vl 



(1) (2) 
V'2 V2 



where S 



solve the same equation (3.15); thus 



( 



(1) (2) 

Vl Vl 



s 



1 (1) (2) 
V '^2 V2 



(2) \ 



(1) (2) 
\ <^2 V2 



solves the equation (3.15) and satisfies the reality condition <I>2 = $1. 
Proof: Direct calculation leads to the conclusion. Q.E.D. 

The periodicity condition $(a; + 27r) = $(a;) can be easily accomplished. For 
example, we can take (^'^^' {j = 1, 2) to be two linearly independent Bloch functions 
ipO) = e'^^^^p^^^ {j = 1,2), where <J2 = — ai and ^p'^^^ are periodic functions ip''^\x + 
27r) = V^-'Ha;)- Often we choose A to be a double point of geometric multiplicity 2, 
so that ip^^^ are already periodic or antiperiodic functions. 

3.2. Discrete Cubic Nonlinear Schrodinger Equation 

Consider the discrete focusing cubic nonlinear Schrodinger equation (DNLS) 

(3.16) iqn = -^[qn+i - 2qn + qn-i] + |g'n|^(gn+i + 9n-i) - 2a;^g„ , 
under periodic and even boundary conditions, 

qn+N = qn I q-n = qn i 

where i = V— 1, g„'s are complex variables, n € Z, w is a positive parameter, 
h = 1/N, and TV is a positive integer N > 3. The DNLS is integrable by virtue of 
the Lax pair [2] : 

(3.17) = L^fVn , 

(3.18) = Bi'^^n, 
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where 




b^n^ -izhqn + {l/z)ihqn-i \ 

-izhqn-i + {l/z)ihqn 

bl^^ = 1 - + 2i\h - h^qnqn-i + to'^h'^, 

6^ = 1/z^ - 1 + 2iXh + h^qnqn-i - 

and where z = exp{iXh). Compatibility of the over-determined system (3.17,3.18) 
gives the "Lax representation" 

of the DNLS (3.16). Let M(n) be the fundamental matrix solution to (3.17), the 
Floquet discriminant is defined as 

A = trace {M{N)} . 

Let ip~^ and 'tjj~ be any two solutions to (3.17), and let W„(V'~'~, V'~) be the Wron- 
skian 

One has 

where pn = 1 + h'^\qn 



w^„(V'+, V-) = V'i+'^Vi-'^^ - V'l+'^Vi-'^^ • 
w^„+i(^+,V-) = P„w^„(V+,^") , 



where = Ili^Lo^ Pn- Periodic and antiperiodic points z*-^-* are defined by 

A(z(f)) = ±2D . 
A critical point z^"^ is defined by the condition 

— -0 

A multiple point z^™) is a critical point which is also a periodic or antiperiodic point. 
The algebraic multiplicity of z'^"^^ is defined as the order of the zero of A(z) ± 2D. 
Usually it is 2, but it can exceed 2; when it does equal 2, we call the multiple point 
a double point, and denote it by z'''^\ The geometric multiplicity of z^™^ is defined 
as the dimension of the periodic (or antiperiodic) eigenspace of (3.17) at z^"), and 
is either 1 or 2. 

Fix a solution g„(t) of the DNLS (3.16), for which (3.17) has a double point 
Qf geometric multiplicity 2, which is not on the unit circle. We denote two 
linearly independent solutions (Bloch functions) of the discrete Lax pair (3.17,3.18) 
at z = z^"*) by ((/)+,(/)"). Thus, a general solution of the discrete Lax pair (3.17;3.18) 
at {qn{t),z^'^^) is given by 

<l)n = C+(l)+ + C~(j)~, 

where c+ and c~ are complex parameters. We use to define a transformation 
matrix r„ by 

z + (l/z)a„ bn 

Cn -1/z + zdn 



18 



3. FIGURE-EIGHT STRUCTURES 



dn 
K 



where, 

1 

|zW|4- 1 
^(d)^(<i)A„ 
1 
IP) 

Prom these formulae, we see that 



(d)|2 



A„ = — 



yral<A'ra2) 



Then we define Q„ and by 
(3.19) 
and 

(3.20) *„(t; ^) = r„(^; ^('^); .^„)V„(t; ^) 

where ■(/;„ solves the discrete Lax pair (3.17,3.18) at {qn{t), z). Formulas (3.19) and 
(3.20) are the Backlund-Darboux transformations for the potential and eigenfunc- 
tions, respectively. We have the following theorem [118]. 

Theorem 3.5. Let qn{t) denote a solution of the DNLS (3.16), for which (3.17) 

has a double point z'^'^^ of geometric multiplicity 2, which is not on the unit circle. 
We denote two linearly independent solutions of the discrete Lax pair (3.17,3.18) 
at {qn, z^^^) by {41+, cj)-). We define Q„{t) and ^„{t; z) by (3.19) and (3.20). Then 

(1) Qn{t) is also a solution of the DNLS (3.16). (The eveness of Q„ can be 
obtained by choosing the complex Bdcklund parameter c'^/c~ to lie on a 
certain curve, as shown in the example below.) 

(2) ^'„(t;z) solves the discrete Lax pair (3.17,3.18) at {Qn{t),z). 

(3) A(z; Qn) = A{z; qn), for all zgC. 

(4) Qn{t) is homoclinic to qn{t) in the sense that Qn{t) e*^* qn{t), expo- 
nentially as exp(— (7|t|) as t ^ ±00. Here 6± are the phase shifts, a is a 
nonvanishing growth rate associated to the double point z^'^\ and explicit 
formulas can be developed for this growth rate and for the phase shifts 9± . 

Example: We start with the uniform solution of (3.16) 



(3.21) 



?n = 9c , Vn; qc = aexp < - i[2(a^ - ui^)t - 7] 



We choose the amplitude a in the range 



iVtan— < o < A'^tan — 
N N 



when N >3 



(3.22) 



3 tan — < a < 00 



when A'' = 3 ; 
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SO that there is only one set of quadruplets of double points which are not on 
the unit circle, and denote one of them hy z = z^"^^ = z^^^ which corresponds to 
(3 = tt/N. The homoclinic orbit Qn is given by 



(3.23) 
where 



Qn = Qc{En+l] 



An+i - 2cos/3-\/pcos2/3 - 1B„+1 



En = ha cos /3 + \/ pcos^ (3 - 1 sech [2^it + 2p] cos[(2n - l)/3 + i?] , 
An+i = ha cos (3 + Vpcos^/? - 1 scch [2/if + 2p] cos[(2n + 3)/? + i?] , 
Bn+i = cos tp + i sin ip tanh[2/it + 2p] + sech [2/if + 2p] cos[2(n + 1)(3 + d] 
P = tt/N, p=l + h'^a'^ , n = 2h^^ y/p sin 13 (3 - 1 , 

h=l/N, c+/c_ = ie^V'' , ??e[0,27r], (-00,00), 
VPC0S/3+ \/pcos2/3- 1 , 6{t)^{a^ -uj^)t--il2 , 
sj p cos^ /? — 1 + i-ypsin/3 = ft.ae*''' , 



where ^p = sm~^[^{ha)~^ sin/3], ip e (0,7r/2). 

Next wc study the "evenness" condition: Q_„ = Qn- It turns out that the 
choices ^ = —(3 , —(3 + tt in the formula of Q„ lead to the evenness of Qn in n. In 
terms of figure eight structure of Qn, = —(3 corresponds to one ear of the figure 
eight, and -& = —(3 + tt corresponds to the other ear. The even formula for Qn is 
given by. 



(3.24) 
where 



Qn 



r/A„ - 1 



r = 1 — cos 2p — i sin 2p tanh[2/it + 2p] , 
A„ = 1 ± cos p{cos f3]~^ sech[2/xt + 2p] cos[2n^] , 

where ('+' corresponds to ■& = —(3). 

The heteroclinic orbit (3.24) represents the figure eight structure. If we denote 
by S the circle, we have the topological identification: 

(figures) 05= IJ Qn{p,l,a,uj,±,N) . 

pe(-oo,oo), 7e[0,27r] 

3.3. Davey-Stewartson II (DSII) Equations 

Consider the Davey-Stewartson II equations (DSII), 

r idtq = [dl - d^]q + [2(|g|2 - a;^) + Uy]q , 
(3.25) <^ 

[ [dl + dl]u = -4dy\q\\ 

where q and u are respectively complex-valued and real-valued functions of three 
variables {t,x,y), and a; is a positive constant. We pose periodic boundary condi- 
tions, 

q{t, x + Li,y) = q{t, x, y) = q{t, x,y + L2) , 
u{t, x + Li,y) = u{t,x,y) = u{t,x,y + L2) , 
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and the even constraint, 



q{t, -x, y) = q{t, X, y) = q{t, x, -y) , 
u{t, -X, y) = u{t, X, y) = u{t, X, -y) . 



Its Lax pair is defined as: 

(3.26) 
(3.27) 

where ip = {ijji,ip2), and 



A = i 



L = 

-dl qd^ 

qdx dl 




(3.28) 



£)+ = ady + dx , 
ri and r2 have the expressions, 

1 



D = ady — dx ■ 



1, 



a 



-1 . 



^[-w + iv] , Ti = -[w + iv], 



(3.29) n 

where u and v are real-valued functions satisfying 

(3.30) [dl + dl]w = 2[dl-dl]\q\\ 

(3.31) [dl + dl]v = Hadxdy\q\'' , 

and w = 2(|g|^ — w^) + Uy. Notice that DSII (3.25) is invariant under the transfor- 
mation cr: 

(3.32) 



(TO {q,q,ri,r2;a) = {q,q,-r2,-ri;-a) . 



Applying the transformation a (3.32) to the Lax pair (3.26, 3.27), we have a con- 
gruent Lax pair for which the compatibility condition gives the same DSII. The 
congruent Lax pair is given as: 



(3.33) 
(3.34) 

where V; = (V'l, ^^2), and 



Q 
-D- 




-r2 -{D-q) 
{D+q) -n 



-dl qdx 
qdx dl 

The compatibility condition of the Lax pair (3.26, 3.27), 

dtL = [A, L] , 

where [A, L] = AL — LA, and the compatibility condition of the congruent Lax pair 
(3.33, 3.34), 

dtL=[A,L] 
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give the same DSII (3.25). Let {q,u) be a solution to the DSII (3.25), and let Aq 
be any value of A. Let V' = {ipi,fp2) be a solution to the Lax pair (3.26, 3.27) at 
(9i 9) ''1) ^2; Ao). Define the matrix operator: 



r 

where A = ady — A, and a, b, c, d are functions defined as: 



A + a b 
c A + d 



b = 



c = 



[V'2 A2 tp2 + Ipl Ai V'l] 
^ [V'2 Ai tpi - Ipi A2 tp2] 
[tpi Ai 1p2 - tlJ2 A2 V'l] 



d = ^ [V'2 Ai V2 + t/Ji A2 V'l] , 
in which Ai = ady — Aq, A2 = ady + Aq, and 

A=-[iV'ir+iv2n . 

Define a transformation as follows: 

/ (9,n,r-2) {Q,Ri,R2), 
\ 4, ^ 

Q = q-2b, 
(3.35) Ri = ri + 2{D+a), 

R2 = T2- 2{D-d) , 
$ = Tcf,; 

where (j) is any solution to the Lax pair (3.26, 3.27) at (g, g, ri, r2; A), £>+ and D~ 
are defined in (3.28), we have the following theorem [121]. 

Theorem 3.6. The transform,ation (3.35) is a Bdcklund-Darboux transforma- 
tion. That is, the function Q defined through the transformation (3.35) is also a 
solution to the DSII (3.25). The function $ defined through the transformation 

(3.35) solves the Lax pair (3.26, 3.27) at (Q, Q, i?2; A). 

3.3.1. An Example . Instead of using Li and L2 to describe the periods of 
the periodic boundary condition, one can introduce ki and K2 as Li = ^ and 
L2 = Consider the spatially independent solution, 

(3.36) qc = rj cyLY>{—2i[rj^ ~ uj'^]t + i"f} . 
The dispersion relation for the linearized DSII at qc is 

n = ± '|^j|^ ^V-(^2^g|) , for 5q ~ qc exp{i(6a; + 6y) + ^ , 

V €1 + S2 
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where = kiKi, ^2 = ^2^2, and ki and are integers. We restrict ki and K2 as 
follows to have only two unstable modes {±ki, 0) and (0, ±K2), 

K2 < Ki < 2k2 , kI < 4?7^ < min{K;i + k\,Ak\} , 



or 



Ki < K2 < 2k;i , K2 < ^ff < min{K;f + k%,^k\} . 
The Bloch eigenfunction of the Lax pair (3.26) and (3.27) is given as, 



(3.37) 
where 



^ = c{t) 



X 



exp{i(6a^ + 6y)} , 



c{t) = Co exp {[2^i(ia^2 - A) + ir2\ t} , 
r-2-ri = 2{\qcf - w'^) , 
X = iia^2 - A) - i^i , 
(ia6 - Xf + ei = ■ 
For the iteration of the Backlund-Darboux transformations, one needs two sets of 



eigenfunctions. First, we choose ^1 = ±\k\, ^2 =0, Aq 
branch), 

(3.38) i)^ = c- 

where 



1 ^2 



(for a fixed 



-(Ic 
,± 



exp 



= Cq exp {[=Ffi;iAo + ir2] t} , 



-Ao T i-^i^i 



We apply the Backlund-Darboux transformations with ^ = + tjj , which gener- 
ates the unstable foliation associated with the {ki, 0) and (— ki, 0) linearly unstable 

modes. Then, we choose ^2 = ±5'«2, A = 0, = ~ \i^2 (for a fixed branch). 



(3.39) 
where 



-9c 



exp \ i{^lx ± -K2y) 



c± = Cj_ exp { [±iaK2S,i + ir2] t} , 
X± = ±ia^K2 - i^i = ±r]e^'^' . 

We start from these eigenfunctions (p± to generate T(p± through Backlund-Darboux 

transformations, and then iterate the Backhmd-Darboux transformations with T(j)j^ + 
T(j)- to generate the unstable foliation associated with all the linearly unstable 
modes (±ki,0) and (0, ±/t2). It turns out that the following representations are 



convenient, 
(3.40) 
(3.41) 
where 



3.3. DAVEY-STEWARTSON II (DSII) EQUATIONS 



,0 „0 f,iC{x+ir2t 



Wo 



and 



-9ce 2 



T = 2KiAot — p , i = TT'^ia; + — ,z = x— — —'&i , 



The following representations are also very useful, 
(3.42) 



t I vi 

V2 



(3.43) 
where 



+ < 



T T T T 

v\ = — gc[cosh — cos x — i sinh — sin i] , V2 = 7?[cosh — cos z — i sinh — sin z\ 



T T T T 

Wi = —gc [cosh — cos^ + i sinh — sin y] , t«2 = ?? [sinh — cos5 + i cosh — sin.2] . 

Applying the Backlund-Darboux transformations (3.35) with f/; given in (3.42), 
have the representations. 



(3.44) a = — Ao sech r sin(a; + z) sin(i — z) 

1 + sech T cos(i + z) cos(x — z) 



(3.45) h = -qcb=- 



V 



cos{x — z) — i tanhr sin(i — z) 



(3.46) 



+ sech r cos(a; + z) 
c=b, d= —a = —a . 



1 + sech r cos(i + z) cos{x — z) 



24 



3. FIGURE-EIGHT STRUCTURES 



The evenness of 6 in a; is enforced by the requirement that d — -di = ±f , and 
(3.47) = =fAo sech r cos^?i sin(Kia;) 



(3.48) 



(3.49) 



1 =F sech T sin cos(Kia;) 
Ao^c 



V 



± sech rcos(K;ia;) 



sin I?! — i tanh r cos i?! 



1 =F sech r sini?i cos(K;ia;) 



c = b, d= —a = —a. 



Notice also that is an odd function in x. Under the above Backlund-Darboux 
transformations, the eigenfunctions (p± (3.39) are transformed into 



(3.50) 



where 



A + a b 
b A - a 



and A = ady — A with A evaluated at 0. Let (p = Lp^ + p (the arbitrary constants 
Cj_ are already included in (p^), (p has the representation, 



(3.51) 



where 



VW2 



Wi = {adyWi) + awi + ribw2 , W2 = {adyW2) — aw2 + rjbwi . 

We generate the coefficients in the Backlund-Darboux transformations (3.35) with 
if (the iteration of the Backlund-Darboux transformations). 



(3.52) a(^) 

(3.53) 
(3.54) 



qc 

V 



W2{adyW2) + Wi{adyWi) 



W2{adyWi)-Wi{adyW2) 



Wlf + \W2\'^ 
\Wif + \W2f 



c(^) = bin , d(^) 
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where 



W2{adyW2) +Wi{adyWi) 



—aK2<. coshf 



+ 



aK2a + iar]{b + b) cos ■&2 
cos{y + z) sin ■i?2 + sinh f 



ar]{b — b) sin 7^2 



\WiY + \W2 

coshf 



+ +^^1^1^ +iaK2r}^(]} + b)cos'&2 



1 



V'\b\' 



sin(y + z) sin-(?2 + sinhf 



W2iadyWi)-Wi{adyW2) 



= -a«;2'| coshf 



an2'rjb + i{—a'^ + —K2 + Tf^b^) cos 7^2 



+ sinhf (a^ - -K2 + f] b ) sini?2| • 

The new solution to the focusing Davey-Stewartson II equation (3.25) is given by 

(3.55) Q = qc-2b- 2b^^^ . 

The evenness in y is enforced by the requirement that 1? — '!?2 = ±f • In fact, 

we have 

Lemma 3.7. Under the requirements that "& — "di = and 1? — i?2 = ±f , 

(3.56) b{-x) = b{x), b^'\-x,y) = b^'\x,y) = b^'\x,-y), 
and Q = qc — 2b— 2b^^^ is even in both x and y. 

Proof: It is a direct verification by noticing that under the requirements, a is 
an odd function in x. Q.E.D. 

The asymptotic behavior of Q can be computed directly. In fact, we have the 
asymptotic phase shift lemma. 

Lemma 3.8 (Asymptotic Phase Shift Lemma). For Aq > 0, > 0, andia = 1, 
as t ^ ±00, 



(3.57) 



Q = qc-2b- 26^^) ^ g,e''^e=F'2(''i-^2) 



In comparison, the asymptotic phase shift of the first application of the Bdcklund- 
Darboux transformations is given by 

3.4. Other Soliton Equations 

In general, one can classify soliton equations into two categories. Category I 
consists of those equations possessing instabilities, under periodic boundary condi- 
tion. In their phase space, figure-eight structures (i.e. separatrices) exist. Category 
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II consists of those equations possessing no instability, under periodic boundary con- 
dition. In their phase space, no figure-eight structure (i.e. separatrix) exists. Typ- 
ical Category I soliton equations are, for example, focusing nonlinear Schrodinger 
equation, sine-Gordon equation [127], modified KdV equation. Typical Category II 
soliton equations are, for example, KdV equation, defocusing nonlinear Schrodinger 
equation, sinh-Gordon equation, Toda lattice. In principle, figure-eight structures 
for Category I soliton equations can be constructed through Backlund-Darboux 
transformations, as illustrated in previous sections. It should be remarked that 
Backlund-Darboux transformations still exist for Category II soliton equations, but 
do not produce any figure-eight structure. A good reference on Backlund-Darboux 
transformations is [152]. 



CHAPTER 4 



Melnikov Vectors 

4.1. ID Cubic Nonlinear Schrodinger Equation 

We select the NLS (3.1) as our first example to show how to establish Melnikov 
vectors. We continue from Section 3.1. 

Definition 4.1. Define the sequence of functionals Fj as follows, 
(4.1) F,{q) = A{X'^{q),q), 

where A^'s arc the critical points, q — (g, —q). 

We have the lemma [137]: 

Lemma 4.2. IfXj{q) is a simple critical point o/A(A) [i.e., A"(Ap ^ OJ, Fj is 
analytic in a neighborhood of q, with first derivative given by 



(4.2) 

where 

(4.3) 



6F, 


SA 




OA 


5X^ _ 


SA 


6q 


" 5q 






A=Ac Sq 


" Sq 



S - 4 



V'2" {x; X)7p2 {x; A) 
tpt{x;X)7p^ (a;; A) 



and the Bloch eigenfunctions ip^ have the property that 
(4.4) {x + 27r; A) = p± {x; A) , 

for some p, also the Wronskian is given by 

w^[v+,^-] = vfV'2"-V'2+V'r • 

In addition, A' is given by 
.^ dA . VA 

^^■'^ Ix = w 



2 _ 4 f^TT 



Proof: To prove this lemma, one calculates using variation of parameters: 

5M{x;X)=M{x) [ M-\x')SQ{x')M{x')dx' , 
Jo 



where 



6Q = i 



Sq 
Sq 

27 
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Thus; one obtains the formula 



(5A(A; q) = trace 



M(2 



(.27r 

tt) / M-^{x' 
Jo 



) 5Q{x') M{x')dx' 



which gives 



(4.6) 



(5A(A) 
5q{x) 

(5A(A) 
(5g(a;) 



i trace 



i trace 



M-i(x)( U J )M(x)M(27r) 



M-i(a;)( J [J ^M{x)M{2tt) 



Next, we use the Bloch eigenfunctions {V'^} to form the matrix 



N{x; A) = 



_( i,+ i,- 



i^t ^-2 



Clearly, 

N{x;\) = M(a;; A)iV(0; A) ; 

or equivalently, 

(4.7) M{x; A) = N{x; \)[N{<d- X)]'^ . 

Since V'^ are Bloch eigenfunctions, one also has 

P 



iV(a; + 27r; A) = Ar(a;; A) 







which implies 

that is, 
(4.8) 



N{2'k; A) = M(27r; A)Ar(0; A) = N{Q; A) 



P 

p-i 



M(27r;A)=7V(0;A)( J ) [iV(0; A)]"! . 



For any 2x2 matrix cr, equations (4.7) and (4.8) imply 

trace {[M(.t)]-i a M{x) M(27r)} = trace {[N{x)\-'^ a N{x) diag{p, p-^)}, 

which, through an explicit evaluation of (4.6), proves formula (4.3). Formula (4.5) 
is established similarly. These formulas, together with the fact that A'^(g) is dif- 
ferentiable because it is a simple zero of A', provide the representation of 
Q.E.D. " 

Remark 4.3. Formula (4.2) for the ^ is actually valid even if A"(A|(g); q) = 
0. Consider a function at which 

A"(A,^(g;);g;) = o, 

and thus at which Xj{q) fails to be analytic. For q near one has 

A'(A5(j);s} = 0j 

AA'(A;(,1;,-) = A"iA; + AA' = 0; 
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that is, 



Thus, 



Since ^ 



Sq ^ A" Sq ' 



5/' 



Sq ^ Sq "' 



■5 A ^' ^ A/ , 



0, as g ^ Q*, one still has formula (4.2) at q = q^: 
S S 



(4.9) 
where 



Sq ' Sq 
The NLS (3.1) is a Hamiltonian system: 

5H 



nt = 



Sq 



H = r{-\q.? + \q\^] dx . 
Jo 



Corollary 4.4. For any fixed A e C, A(A, q) is a constant of motion of the 
NLS (3.1). In fact, 

{A{\,q),H{ci)} = Q , {A(A,g-),A(A',g)} = , VA,A'eC, 

where for any two functionals E and F, their Poisson bracket is defined as 



{E,F} 



27T 



SE SF SE SF' 
Sq Sq Sq Sq 



dx 



Proof: The corollary follows from a direction calculation from the spatial part 
(3.2) of the Lax pair and the representation (4.3). Q.E.D. 

For each fixed q, A is an entire function of A; therefore, can be determined by 
its values at a countable number of values of A. The invariance of A characterizes 
the isospectral nature of the NLS equation. 

Corollary 4.5. The functionals Fj are constants of motion of the NLS (3.1). 
Their gradients provide Melnikov vectors: 



(4.10) 



grad Fj (q) = i 



W['ip+,^l;- 



i;+{x;X-j)i;,{x;X-j) 
i,+ {x;\'^)i,^{x;X'^) 



The distribution of the critical points A^ are described by the following counting 
lemma [137], 

Lemma 4.6 (Counting Lemma for Critical Points). For q G , set N = 
N{\\q\\,)€Z+ by 



jo cosh ( I 



alio 



3||g||isinh ||g||o 



where [x] = first integer greater than x. Consider 



A'{X;q) = -A{X;q). 



Then 
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(1) A'{X;q) has exactly 2N + 1 zeros (counted according to multiplicity) in 
the interior of the disc D = {X € C : |A| < {2N + l)f }; 

(2) Vfc G Z, \k\ > N,A'{X,q) has exactly one zero in each disc 
{AG C: |A-fc7r| < f}. 

(3) A'(A;g) has no other zeros. 

(4) For \X\ > {2N + l)f , the zeros of A',{X''j,\j\ > N}, are all real, simple, 
and satisfy the asymptotics 

Xj=jn + o{l) as\j\^oo. 



4.1.1. Melnikov Integrals. When studying perturbed integrable systems, 

the figure-eight structures often lead to chaotic dynamics through homoclinic bifur- 
cations. An extremely powerful tool for detecting homoclinic orbits is the so-called 
Melnikov integral method [158], which uses "Melnikov integrals" to provide esti- 
mates of the distance between the center-unstable manifold and the center-stable 
manifold of a normally hyperbolic invariant manifold. The Melnikov integrals are 
often integrals in time of the inner products of certain Melnikov vectors with the 
perturbations in the perturbed integrable systems. This implies that the Melnikov 
vectors play a key role in the Melnikov integral method. First, we consider the 
case of one unstable mode associated with a complex double point v, for which the 
homoclinic orbit is given by Backlund-Darboux formula (3.7), 

Q{x,t) = q{x,t) + 2{y-U) ^ /"^^^ - , 



where q lies in a normally hyperbolic invariant manifold and <j) denotes a general 
solution to the Lax pair (3.2, 3.3) at (g, v), (j) = c+0"'" + C-(f)~ , and (f)^ are Bloch 
eigenf unctions. Next, we consider the perturbed NLS, 

iqt = Ixx + 2|gpg + ief{q, q) , 

where e is the perturbation parameter. The Melnikov integral can be defined using 
the constant of motion Fj, where Xj = i' [137]: 

(4.11) M,.y_^ I {^f+^f},.^dxdt, 

where the integrand is evaluated along the unperturbed homoclinic orbit q = Q, and 

the Melnikov vector has been given in the last section, which can be expressed 
rather explicitly using the Backlund-Darboux transformation [137] . We begin with 
the expression (4.10), 

5Fj _ , \^A23I / 



^^'^^^ 6q V[$+,$-] V ^1^1 

where $^ are Bloch eigenfunctions at {Q,y), which can be obtained from Backlund- 
Darboux formula (3.7): 

^^{x,t-,v) = G{v]v\(t)) 4)^{x,t\v) , 

with the transformation matrix G given by 

G = G{X;v;cI>)=n( ^ " X-v)^~'^ 
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Ni 



These Backlund-Darboux formulas are rather easy to manipulate to obtain explicit 
information. For example, the transformation matrix G(A, v) has a simple limit as 



(4.13) 

where is defined by 



lim G(A, u) 

A— 



<j)2<t>l 



-0102 
0101 



= 0101 + 0202- 



With formula (4.13) one quickly calculates 

= ±Czp iy[0+,0-] 



from which one sees that and $ are linearly dependent at {Qjv), 

$+ = -— 

c+ 

Remark 4.7. For Q on the figure-eight, the two Bloch eigenfunctions are 
linearly dependent. Thus, the geometric multiplicity of v is only one, even though 
its algebraic multiplicity is two or higher. 



(4.14) 



Using L'Hospital's rule, one gets 

VA^ -4 _ ^A{u)A"{u) 



With formulas (4.12, 4.13, 4.14), one obtains the explicit representation of the 
grad Fj [137]: 



(4.15) 



5F^ 
Sq 



where the constant Ci, is given by 



With these ingredients, one obtains the following beautiful representation of the 
Melnikov function associated to the general complex double point u [137]: 



/+00 /'27r 
/ W^[0+,0-] 



(0i)/(Q,Q) + (0f)/(Q,Q) 



dxdt. 



In the case of several complex double points, each associated with an instability, 
one can iterate the Backlund-Darboux tranformations and use those functionals 
Fj which are associated with each complex double point to obtain representations 
Melnikov Vectors. In general, the relation between -j4- and double points can be 
summarized in the following lemma [137] , 
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Lemma 4.8. Except for the trivial case q = 0, 

5F SF- 
(a). ^=o^'Ll=o^ M(27r, A^; g,) = ±1. 
oq oq ■' 

5F- 

{b). U = Q A'(A,^(g;); q.) = 0,^ \Fj{q^)\ = 2, 

=> A^(g*) is a multiple point, 
where I is the 2x2 identity matrix. 

The Backlund-Darboux transformation theorem indicates that the figure-eight 

structure is attached to a complex double point. The above lemma shows that at 
the origin of the figure-eight, the gradient of Fj vanishes. Together they indicate 
that the the gradient of Fj along the figure-eight is a perfect Melnikov vector. 

Example: When ^ < c < 1 in (3.9), and choosing = tt in (3.14), one can 
get the Melnikov vector field along the homoclinic orbit (3.14), 

SFi „ . o , [(— sinp + zcosptanhr) cosx + sechrl ,a 

(4.17) = 27r sin^ p sechr r-^ r—^ ^ c e'^ , 

oq [1 — smp sech Tcosa;]^ 

SFi _'SF\ 
6q 6q 

4.2. Discrete Cubic Nonlinear Schrodinger Equation 

The discrete cubic nonlinear Schrodinger equation (3.16) can be written in the 
Hamiltonian form [2] [118]: 

(4.18) iqn = PndH/dqn , 

where p„ = 1 -|- /i^|g„|^, and 

1 f_ 2 1 

H = T2^ \ Qn{qn+1 + qn-l) - -r^il + Uj'^h^)ln Pn > . 

'^n=o ^Qn{qn+i+qn-i)^ itsclf is also a constant of motion. This invariant, together 
with H, implies that J2n=o ^ constant of motion too. Therefore, 

N-l 

(4.19) D^=l[pn 

n=0 

is a constant of motion. We continue from Section 3.2. Using D, one can define a 
normalized Floquet discriminant A as 

A = A/D . 

Definition 4.9. The sequence of invariants Fj is defined as: 

(4.20) F,{q} = A{zf{q);q), 

where g = {q, -q), q = (go,gi, • • ■,qN-i)- 

These invariants Fj^s are perfect candidate for building Melnikov functions. 
The Melnikov vectors are given by the gradients of these invariants. 
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Lemma 4.10. Let z^^\q) be a simple critical point; then 



(4.21) 



(4.22) 



(5A 



w/iere V'n = ('^i'n^'"^\ V'n^'^^)^ '^'^e *w SZoc/i functions of the Lax pair (3.2,3.3), 
such that 

where •ip^ are periodic in n with period N, Wn = det (V'^ , V'n ) • 

For Zj''^ = z^'^\ the Melnikov vector field located on the heteroclinic orbit (3.19) 
is given by 



(4.23) 



5Qn 



En An 



+ 1 



where Qn = {Qn,-Qn), 

Wn = \ \ , 

i5„ = |^(l)|2 + |,(^)|2|^(2)|2^ 



K ^ — 



-|zW|4(|2(d)|4_ i)[^(d)j-lWA(^(d);q-)A"(0('');5) , 



where 



A"(z('');g) 



a^A(z('^);g) 
dz^ 



ExEunple: The Melnikov vector evaluated on the heteroclinic orbit (3.23) is 
given by 



(4.24) 
where 



5Fi _ 



K 



EnAn+l 



sech [2ijt + 2p] 



(1) 



(2) 



cos 13 sech [2^lt + 2p] + cos[(2n + l)/3 + 1? + 



-itaMh.[2iJLt + 2p] sin[(2n + 1)/? + + (p] 



i26(t) 



COS /? sech [2/ii + 2p] + cos[(2n + 1)/? + 1? - 



-i tanh[2/ii + 2p] sin[(2n +l)(3 + 'd-ip\ 



-i2e(t) 



k = -2Nh''a{l - ^*)[8p^/2^2]-Vpcos2/3-l . 
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Under the even constraint, the Melnikov vector evaluated on the heterochnic orbit 
(3.24) is given by 



(4.25) 
where 



SFi 



5Q„ 



(l,e) 



even 



= sech[2/xt + 2p] [n„]-i ^"(2, 



e) 



k^^^ = -2N{1 - ^4)[8ap=^/202]-i^pcos2/J-l , 
n„, = cos /? ± COS (f sech[2/zt + 2p] cos[2(n — 1)/?] 
cos /3 ± cos sech[2/it + 2p] cos[2(n + l)/3] 



cos /3 sech [2/if + 2p] ± (cos if 



—i sin tanh[2/it + 2p]) cos[2n/?] 



cos (3 sech [2;uf + 2p] ± (cos if 
+i sin (f ta.nh[2idt + 2p]) cos[2n/3] 



-i2e{t) 



4.3. Davey-Stewartson II Equations 

The DSII (3.25) can be written in the Hamiltonian form, 



(4.26) 
where 



SH/Sq , 
-6H/Sq 



f"^" f^" 2 2 i 2 



We have the lemma [121]. 

Lemma 4.11. The inner product of the vector 

'>p2'4'2 



u 



V'lV'i 



where tjj = (V'i,V'2) is an eigenfunction solving the Lax pair (3.26, 3.27), and 
■ip = ('^i'l, V'2) is an eigenfunction solving the corresponding congruent Lax pair (3.33, 
1 



3.34), and S 



1 



of (4.26) vanishes. 



, with the vector field JVH given by the right hand side 



{U , JVH) = , 



where J ■ 



1 
-1 
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If we only consider even functions, i.e., q and u = r2 — ri are even functions in 
both X and y, then we can spht U into its even and odd parts. 



{e,y) (o,y) {e,v) (o,y) ' 



where 



j,{e,x) 


1 

4 




-x,y) 


+ U{x, 


-y)+u{- 


-X, 


-V) 




1 
4 


U{x,y)+U{- 


-x,y) 


-U{x, 


-y)-U{- 


-X, 


-y) 




1 
4 


U{x,y)-U{- 


-x,y) 


+ U{x, 


-y)-U{- 


-X, 


-y) 


(o,V) 


1 

4 


U{x,y)-U{- 


-x,y) 


-U{x, 


-y)+u{- 


-X, 


-y) 



Then wc have the lemma [121]. 



Lemma 4.12. When q and u = r2 — ri are even functions in both x and y, we 



have 



(w£g,JVif) = o. 



(4.27) 



4.3.1. Melnikov Integrals. Consider the perturbed DSII equation, 
r idtq = [dl - d^]q + [2{\qf - tu^) + Uy]q + eif , 



-Ady \q\ 



where q and u arc respectively complex-valued and rcal-valucd functions of three 
variables [t, x, y), and G = (/, /) are the perturbation terms which can depend on 
q and q and their derivatives and t, x and y. The Melnikov integral is given by 
[121], 



M 



/OO 
{U,G) 
-OO 



dt 



OO 

OO r2-K /'27r 



(4.28) 



OO ^0 ^0 



Re[[i)2i)2)f + {ipii)i)f} dxdydt, 



where the integrand is evaluated on an unperturbed homoclinic orbit in certain 
center-unstable (= center-stable) manifold, and such orbit can be obtained through 
the Backlund-Darboux transformations given in Theorem 3.6. A concrete example 
is given in section 3.3.1. When we only consider even functions, i.e., q and u are 
even functions in both x and y, the corresponding Melnikov function is given by 
[121], 



/OO 
{U, G) dt , 
-OO 

which is the same as expression (4.28). 
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4.3.2. An Example . We continue from the example in section 3.3.1. We 
generate the following eigenfunctions corresponding to the potential Q given in 
(3.55) through the iterated Backlund-Darboux transformations, 



(4.30) 

(4.31) 
where 



>'i I 



r = 



A + a 



at A = Ao = y ry^ 

at A = , 



b A - a 
where A = ady — A for general A. 

Lemma 4.13 (see [121]). The eigenfunctions "^/^ and $± defined in (4-30) and 
(4-31) have the representations, 



(4.32) 

(4.33) 
where 



±2AoVK(V+,V'") 



(-Ao + a(^))^2-^^'Vi 
600^2 + (Ao+^)V^i 



^1 

^2 



-miCo Cg gcexp{z2r2t} , 



A = - 



I |2 , I |2 



A(^) = 

= ai A = , 

Si = ^i((^+v.f) + ip+{(f+ip^) + ip^{if^ip+) , 

If we take r2 to be real (in the Melnikov vectors, r2 appears in the form r2 — ri = 

2(|q'c|^ — <^'^)), then 

(4.34) ^-^^ ^ , $d 



, as t ^ ±00 . 



Next we generate eigenfunctions solving the corresponding congruent Lax pair 
(3.33, 3.34) with the potential Q, through the iterated Backlund-Darboux trans- 
formations and the symmetry transformation (3.32) [121]. 



Lemma 4.14. Under the replacements 
a — > — a (■i?2 — >7r — ■!?2), ^- 
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the coefficients in the iterated Bdcklund-Darboux transformations are transformed 
as follows, 

Lemma 4.15 (see [121]). Under the replacements 

a<-^-a (i?2 — >7r-i?2), ri<-^—r2, 

(4.35) 

r2 -ri , I? — > 1? + TT - 2'&2 , P — > -P , 
the potentials are transformed as follows, 

Q , 

{R = Q)^iR = Q), 
Ri -R2 , 
R2 — > —Ri ■ 



The cigcnfunctions and $± given in (4.32) and (4.33) depend on the vari- 
ables in the replacement (4.35): 

*± = ^^{a,ri,r2,S,p), 

$± = $±(a,ri,r2,i9,/5) . 

Under replacement (4.35), ^'^ and <1>± arc transformed into 

(4.36) $± = «'±(-a,-r2,-ri,t? + 7r-2i92,-p), 

(4.37) $± = $±(-a,-r2,-ri,'!? + 7r- 2i?2,-/5) . 

Corollary 4.16 (see [121]). and ^± solve the congruent Lax pair (3. S3, 
3.34) at {Q,Q,Ri,R2;Xo) and {Q,Q,Ri,R2;0), respectively. 

Notice that as a function of 77, has two (plus and minus) branches. In 
order to construct Melnikov vectors, we need to study the effect of the replacement 

Lemma 4.17 (sec [121]). Under the replacements 

(4.38) {^2 ^-^2), + IT -2^2, p~^-p, 

the coefficients in the iterated Bdcklund-Darboux transformations are invariant, 

(c(^) = m)^ (c(^) = W)) , (d(^) = -^) ^ (d(^) = -^) ; 
thus the potentials are also invariant, 

Q^Q, {R = Q)^iR^Q), 
Ri — > Ri , i?2 — > R2 ■ 
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The eigenfunction <i?± given in (4.33) depends on the variables in the replace- 
ment (4.38): 

Under the replacement (4.38), $± is transformed into 

(4.39) $± = $±(-ei",'(9 + 7r-2?92,-p) . 

Corollary 4.18. $± solves the Lax pair (3.26,3.27) at (Q, Q, i?2 ; 0). 

In the construction of the Melnikov vectors, we need to replace $± by $± to 
guarantee the periodicity in x of period ii = 

The Melnikov vectors for the Davey-Stewartson II equations are given by, 



(4.40) U 



(4.41) U± = 




where S = [ \^ 1 . The corresponding Melnikov functions (4.28) are given by. 



1 

= I (U^ , G) dt 



/OO 
-OO 

/CO /'27r />27T f 
/ / i?e [4'^*^]/(Q,g) 
OO Jo Jo L 

(4.42) +[*f$f]7(0,g)|dxdydt, 

/OO 
{U± , G) dt 
-OO 

/OO /'27r /'27r 
-ooJO Jo l 

(4.43) +[$2^$i'^]7(Q,Q)|dxdydi, 

where Q is given in (3.55), is given in (4.32), <l>± is given in (4.33) and (4.39), 
is given in (4.32) and (4.36), and $± is given in (4.33) and (4.37). As given in 
(4.29), the above formulas also apply when we consider even function Q in both x 
and y. 



CHAPTER 5 



Invariant Manifolds 

Invariant manifolds have attracted intensive studies which led to two main 
approaches: Hadamard's method [78] [64] and Perron's method [176] [42]. For 
example, for a partial differential equation of the form 

dfU = Lu + N{u) , 

where L is a linear operator and N{u) is the nonlinear term, if the following two 
ingredients 

(1) the gaps separating the unstable, center, and stable spectra of L are large 
enough, 

(2) the nonlinear term N{u) is Lipschitzian in u with small Lipschitz constant, 

are available, then establishing the existence of unstable, center, and stable mani- 
folds is rather straightforward. Building invariant manifolds when any of the above 
conditions fails, is a very challenging and interesting problem [124]. 

There has been a vast literature on invariant manifolds. A good starting point 
of reading can be from the references [103] [64] [42]. Depending upon the emphasis 
on the specific problem, one may establish invariant manifolds for a specific flow, 
or investigate the persistence of existing invariant manifolds under perturbations 
to the flow. 

In specific applications, most of the problems deal with manifolds with bound- 
aries. In this context, the relevant concepts arc overflowing invariant, inflowing 
invariant, and locally invariant submanifolds, defined in the Chapter on General 
Setup and Concepts. 

5.1. Nonlinear Schrodinger Equation Under Regular Perturbations 

Persistence of invariant manifolds depends upon the nature of the perturbation. 

Under the so-called regular perturbations, i.e., the perturbed evolution operator is 
close to the unperturbed one, for any fixed time; invariant manifolds persist 

"nicely". Under other singular perturbations, this may not be the case. 

Consider the regularly perturbed nonlinear Schrodinger (NLS) equation [136] 

[140], 

(5.1) iqt = + 2[|(?|2 - u;^]q + ie[dlq - aq + P] , 

where q = q{t, x) is a complex-valued function of the two real variables t and x, 
t represents time, and x represents space. q{t,x) is subject to periodic boundary 
condition of period 27r, and even constraint, i.e., 

q{t, x + 2ir) = q{t, x) , q{t, -x) = q{t, x) . 
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w is a positive constant, a > and /3 > are constants, is a bounded Fourier 
multiplier, 

AT 

^x1 = -'^ k'^^kQk COS kx , 
fe=i 

^fc = 1 when k < N, S^k = 8k~^ when k > N, for some fixed large A?", and e > is 
the perturbation parameter. 

One can prove the following theorems [136] [140]. 

Theorem 5.1 (Persistence Theorem). For any integers k and n (1 < k,n < 

oo ), there exist a positive constant eq o.'f^d a neighborhood U of the circle = 
{<? I dxQ = 0, |g| = u, 1/2 < uj < 1} in the Sobolev space , such that inside U, 
for any e G (— eo,eo), there exist C" locally invariant submanifolds and W^^ 
of codimension 1, and (= W^" D W^''^ ) of codimension 2 under the evolution 
operator F^ given by (5.1). When e = 0, Wq", Wq" , and Wq are tangent to the 
center-unstable, center- stable, and center subspaces of the circle of fixed points So,, 
respectively. W™, W^^, and are C" smooth in e for e G (— eo, eo). 

VKf", ly^*, and W^r arc called persistent center-unstable, center-stable, and 
center submanifolds near imdcr the evolution operator _F* given by (5.1). 

Theorem 5.2 (Fiber Theorem). Inside the persistent center-unstable subman- 
ifold W™ near S^^, there exists a family of 1- dimensional C" smooth submanifolds 
(curves) {.^"^"''^^(g) : q G W^}, called unstable fibers: 

• W^'^ can be represented as a union of these fibers, 

W™= IJ J-("'^)(g). 

• J^^'^''^\q) depends C"~^ smoothly on both e and q for e G (— eo,eo) and 
q e W^, in the sense that W defined by 



|(gi,g,e) gi e^("'^)(g), qGW^, eG(-eo,eo)| 



is a C"~^ smooth submanifold of H'' x H'' x (— eo,eo). 

Each fiber J^^^''^\q) intersects transversally at q, two fibers T^'^'^^qi) 

anrf ^^"'"^^(52) o.re either disjoint or identical. 

The family of unstable fibers {.7-"^"'^^ (g) : q G W^} is negatively invariant, 
in the sense that the family of fibers commutes with the evolution operator 
Fl in the following way: 

F*(.F("'^)(g))c.F("^^)(F*(g)) 

for all q G and allt <0 such that \Jre[t,o] Plil) ^ 

There are positive constants k and C which are independent of e such that 

ifqG Wl and qi G J^^"''^?); then 



Fl{q{)-Fl{q) 



< Ce 



fit 
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for all t < such that Urelf.o] ^eil) C W^, where 
of periodic functions of period 2-k. 



denotes 
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For any q,p & W^, q ^ p, any qi e jF("'')(g) and any pi G J^^"''^\p); if 
F'^iq),F'^{p)eW^, ViG(-oo,0], 



and 
then 

Similarly for 



e5«* ^0, ast^ -oo. 



\F!{q,)~F^{q)\\ \ 
\F^ip,)-Fi{q)\\j 



When e = 0, certain low-dimensional invariant submanifolds of the invariant 
manifolds, have explicit representations through Darboux transformations. Specif- 
ically, the periodic orbit (3.9) where 1/2 < c < 1, has two-dimensional stable and 
unstable manifolds given by (3.14). Unstable and stable fibers with bases along the 
periodic orbit also have expressions given by (3.14). 

5.2. Nonlinear Schrodinger Equation Under Singular Perturbations 

Consider the singularly perturbed nonlinear Schrodinger equation [124], 
(5.2) iqt = qxx + 2[|g|^ - uJ^]q + iefg^x -aq + 13] , 

where q = q{t,x) is a complex-valued function of the two real variables t and x, 

t represents time, and x represents space. q{t,x) is subject to periodic boundary 
condition of period 27r, and even constraint, i.e., 

q{t, X + 27r) = q{t, x) , q(t, -x) = q(t, x) . 

CO e (1/2, 1) is a positive constant, a > and /3 > are constants, and e > is the 
perturbation parameter. 

Here the perturbation term ed"^ generates the semigroup 6*^*^=, the regularity 
of the invariant mainfolds with respect to e at e = will be closely related to the 
regularity of the semigroup e^*^« with respect to e at e = 0. Also the singular 
perturbation term ed'^q breaks the spectral gap separating the center spectrum and 
the stable spectrum. Therefore, standard invariant manifold results do not apply. 
Invariant manifolds do not persist "nicely" . Nevertheless, certain persistence results 
do hold. 

One can prove the following unstable fiber theorem and center-stable manifold 
theorem [124]. 

Theorem 5.3 (Unstable Fiber Theorem). Let A be the annulus: A = {q\ dxq = 
0, 1/2 < \q\ < 1}, for any p £ A, there is an unstable fiber which is a curve. 

has the following properties: 

(1) !Fp is a smooth in H^-norm, k > 1. 

(2) ^+ is also smooth in e, a, /3, w, and p in H'^-norm, k > 1, e G [0, eo) 

for some eg > 0. 

(3) has the exponential decay property: Vpi S 

\\^y^-^;ph<ce^^, v*<o, 

\\pl-p\\k 

where F* is the evolution operator, fi > 0. 
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(4) {^p}peA forms an invariant family of unstable fibers, 
F*J-+C^+p, yt€[-T,0], 
and \/T>0 (T can be +oo), such that F*p eA,\/te [-T, 0]. 

Theorem 5.4 (Center-Stable Manifold Theorem). There exists a smooth 
codimension 1 locally invariant center-stable m,anifold W^^ in a neighborhood of the 
annulus A (Theorem 5.3) in for any k > 1. 

(1) At points in the subset Wj:^^ ofW^^, W^^ is smooth in e in H^-norm 
for e S [0, eo) and some eg > 0. 

(2) W^" is smooth in (a,^,u;). 

Remark 5.5. regularity in e is crucial in locating a homoclinic orbit. As 
can be seen later, one has detailed information on certain unperturbed (i.e. e = 
0) homoclinic orbit, which will be used in tracking candidates for a perturbed 
homoclinic orbit. In particular, Mclnikov measurement will be needed. Melnikov 
measurement measures zeros of 0(e) signed distances, thus, the perturbed orbit 
needs to be 0{e) close to the unperturbed orbit in order to perform Melnikov 
measurement. 

5.3. Proof of the Unstable Fiber Theorem 5.3 

Here we give the proof of the unstable fiber theorem 5.3, proofs of other fiber 
theorems in this chapter are easier. 

5.3.1. The Setup of Equations. First, write q as 

(5.3) q{t,x) = [p{t) + f{t,x)]e''^'\ 

where / has zero spatial mean. We use the notation (•) to denote spatial mean, 

(5.4) ^'^^'^2^1 

Since the L^-norm is an action variable when e = 0, it is more convenient to replace 
pby: 

(5.5) 7=(|g|2)=p2^(|/|2). 

The final pick is 



(5.6) J = I 



.2 



In terms of the new variables (J, 9, /), Equation (5.2) can be rewritten as 



(5.7) j = e -2a{J + lu^) + 20^.7 + uj^ cos 



sm( 



(5.8) e = -2j-e^^f^ + nl 

(5.9) ft = LJ + VJ - iN2 - i^3, 
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where 
(5.10) 

(5.11) 
(5.12) 

(5.13) 

(5.14) 

(5.15) 



LJ = + e{-a.f + /,,) - i2u;2(/ + /), 

VJ = -t2J{f + f)+tePf ^'""^ 



^2 =-2(|/.|') + 2/3cOS^ V/J + W2-(|/|2)- VJ + w2 



-{{f+m--{\f\'(f+f)) 



1 



1 



■ w 



AA2 = 2p[2(|/|2-(|/|2)) + (f -(f))], 

= -{f + P + 6|/p)/ + 2(|/|V - (|/|2/)) 

-^(l/P(/ + /))/-2(|/H/ 



el3 sine 



/• 



Remark 5.6. The singular perturbation term "eS^fj" can be seen at two loca- 
tions, Lc and 7^2 (5.10,5.12). The singular perturbation term {\fx\'^) in T^i does 
not create any difficulty. Since is a Banach algebra [6], this term is still of 
quadratic order, (|/,|2) ^ 0{\\f\\l). 

Lemma 5.7. The nonlinear terms have the orders: 

\T^i\-oi\\m, \n',\^o{\\.m), 

i!A/'2iu-o(ii/ii^), imis-oiwm, is>i). 

Proof: The proof is an easy direct verification. Q.E.D. 

5.3.2. The Spectrum of L^. The spectrum of consists of only point spec- 
trum. The eigenvalues of ie are: 



(5.16) = -eia + k^)±kVAuj'^-k'^, (fc = l,2,...); 

where u! G (5? l)j only ^if are real, and /i^ arc complex for fc > 1. 

The main difficulty introduced by the singular perturbation ed^f is the breaking 
of the spectral gap condition. Figure 1 shows the distributions of the eigenvalues 
when e = and e ^ 0. It clearly shows the breaking of the stable spectral gap 
condition. As a result, center and center-unstable manifolds do not necessarily 
persist. On the other hand, the unstable spectral gap condition is not broken. 
This gives the hope for the persistence of center-stable manifold. Another case of 
persistence can be described as follows: Notice that the plane 11, 



(5.17) 



n = {<z| 9.g = 0}, 



is invariant under the flow (5.2). When e = 0, there is an unstable flbration with 
base points in a neighborhood of the circle S^, of fixed points, 



(5.18) 



S^ = {qeU\ \q\=uj}, 
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e = 



e > 



Figure 1. The point spectra of the Unear operator L^. 



in n, as an invariant sub-fibration of the unstable Fenichel fibration with base points 
in the center manifold. When e > 0, the center manifold may not persist, but 11 
persists, moreover, the unstable spectral gap condition is not broken, therefore, 
the unstable sub-fibration with base points in 11 may persist. Since the semiflow 
generated by (5.2) is not a perturbation of that generated by the unperturbed 
NLS due to the singular perturbation e9|, standard results on persistence can not 
be appHed. 

The eigenfunctions corresponding to the real eigenvalues are: 
(5.19) ef =e±^''^cosx, ,±^s. ^ItW^EI. 

Notice that they are independent of e. The eigenspaces corresponding to the com- 
plex conjugate pairs of eigenvalues are given by: 

El = spandcoscc}. 

and have real dimension 2. 

5.3.3. The Re-Setup of Equations. For the goal of this subsection, we need 
to single out the eigen-directions (5.19). Let 



where are real variables, and 

{h) — {hcosx) = 0. 
In terms of the coordinates {^f, J, 0, h), (5.7)-(5.9) can be rewritten as: 



(5.20) 


it 


= f^tit + v{'^t+^ft, 


(5.21) 


j 


= e 


-2a{J + u)'^) + 2(3^ J + u)^ cos 


(5.22) 







^ sin^ 


(5.23) 


ht 


= L^h + V.h+M, 


(5.24) 


4" 





5.3. PROOF OF THE UNSTABLE FIBER THEOREM 5.3 45 

where Hi are given in (5.16), A/f*^ and J\f are projections of —iAf2 — iJ^s to the 
corresponding directions, and 

sin Q 

Ff = -2ci + ^r) + e/3^==(cr^i+ - c+^r). 

1 _,_ _ 2a;2 - 1 _ 2up- 



V4a;2 - 1 ' ^ \/4a;2 - 1 ' ^ \/4a;2 - 1 ' 

5.3.4. A Modification. 

Definition 5.8. For any 5 > 0, we define the annular neighborhood of the 
circle (5.18) as 

A{6) = {{J,d) \ |J|<<5}. 

To apply the Lyapunov-Perron's method, it is standard and necessary to modify 
the J equation so that ^(4^) is overflowing invariant. Let r] G C'^{R,R) be a 
"bump" function: 

'O, in (-2,2) U (-00,-6) U (6,oo), 
1, in (3,5), 
^-1, in (-5,-3), 

and \r]'\ <2, \r]"\ < C. Wc modify the J equation (5.21) as follows: 

(5.25) j = ebn{J/5) + e \-2a{J + w^) + 2/^^7+7^ cos e\ + eTZ(, 

where b > 2(2q:cj^ + 2/3uj). Then A{45) is overflowing invariant. There are two 
main points in adopting the bump function: 

(1) One needs A{4:d) to be overflowing invariant so that a Lyapunov-Perron 
type integral equation can be set up along orbits in A{4:5) for t g (—00, 0). 

(2) One needs the vector field inside A{2S) to be unchanged so that results 
for the modified system can be claimed for the original system in A{S). 

Remark 5.9. Due to the singular perturbation, the real part of /x^ approaches 
—00 as fc — » 00. Thus the h equation (5.23) can not be modified to give overflowing 
flow. This rules out the construction of unstable fibers with base points having 

general h coordinates. 

5.3.5. Existence of Unstable Fibers. For any (Jo, 6*0) S A{4:S), let 

(5.26) j = j^(t), e^e^{t), te(-oo,o], 

be the backward orbit of the modified system (5.25) and (5.22) with the initial 
point {Jo,0o). If 

i^tit), Mt) + Jit),e4t) + m,h{t), ^r(i)) 

is a solution of the modified full system, then one has 

(5.27) it = ^^t^t + Ft, 

(5.28) ut = Au + F, 
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where 



u = 







(' 








0^ 




/Fj\ 


9 















, F = 


Fe 


h 















Fh 














Mr/ 







F+ = V+^t+^fi^, 

Fj = eb[v{J/S)-v{Mt)/S)]+e 



- 2P y/jjt)+^ cos e^t) 



2a J + 20^/7+1^ cose 



Fe = -e/3 



sinf 



+ e/3- 



sin0^ (t) 



;2 ^jJf)TZ? 
Fh = V,h + N, 

Fr = v,-^^+Mr, 

j = j^{t) + j, e = e4t) + e. 

System (5.27)-(5.28) can be written in the equivalent integral equation form: 



(5.29) 
(5.30) 



J to 

u{t) = e^(*-*<')u(io) + f e^^*-^^F{T)dT. 

Jto 



By virtue of the gap between /i^ and the real parts of the eigenvalues of A, one can 
introduce the following space: For a e and n > 1, let 

Ga,n = {g{t) = {it{t),u{t)) 



t £ (—00,0], g{t) is continuous 

in t in iJ" norm , \\g\\a,n = supe"'^*[|^j'"(t)| + ||u(t)||„] < oo I . 

t<o J 

Gf^^n is a Banach space under the norm || • ||(7,ra- Let B^^n{r) denote the ball in Go-,n 
centered at the origin with radius r. Since A only has point spectrum, the spectral 
mapping theorem is valid. It is obvious that for i > 0, 

\\e^'u\U<C{l+t)\\u\\n, 

for some constant C. Thus, if g{t) S Bcr^nir), r < oo is a solution of (5.29)-(5.30), 
by letting to —oo in (5.30) and setting to = in (5.29), one has 



(5.31) 
(5.32) 



^+(t) =^+(0)e'^^*+ / e'^^(*-^)i^+(r)rfr, 

u{t)= j e^^*-^^F{T)dT. 

<y —CO 



For g{t) e Ba,n{r), let r{g) be the map defined by the right hand side of (5.31)- 
(5.32). Then a solution of (5.31)-(5.32) is a fixed point of F. For any n > 1 and 
e < d^, and 5 and r are small enough, Fi+ and F arc Lipschitz in g with small 
Lipschitz constants. Standard arguments of the Lyapunov-Perron's method readily 
imply the existence of a fixed point g^ of F in B,T,n{r)- 
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5.3.6. Regularity of the Unstable Fibers in e. The difficulties lie in the 
investigation on the regularity of with respect to (e, a, /3, w, Jg, ^o, C]^(0))- The 
most difRcult one is the regularity with respect to e due to the singular perturbation. 
Formally differentiating in (5.3f )-(5.32) with respect to e, one gets 



(5.33) 

(5.34) 
where 

(5.35) 



(5.36) 

(5.37) 
(5.38) 
(5.39) 



Jo 

Ue{t) = j 



A{t-T) 



duF+ . + d^+F+ . ^+ ] {T)dT + 7^+(^), 

duF -u. + d^+F- Ci+J {T)dT + n{t), 



^0 

+ / e^^ [d,F+ + du, F+ • {T)dT, 
Jo 

TZ{t)= f (i - T)A,e^(*-^)F(T)rfT 

J — oo 

+ r e^^'-^\deF + du,F-u,,e]{r)dT, 

J — oo 



\ 





-(a + 1)/ 



-(a + 1), 

/O 

-a + dl 
\0 

where T = transpose, and (J*, 6'*) are given in (5.26). The troublesome terms are 
the ones containing Ae or in (5.35)-(5.36). 

(5.40) \\AF\\^ < C\\F\\n+2 < d{\\u\\n+2 + |^i+|), 

where c is small when ( • ) on the right hand side is small. 

dj.Fj ■ J,,, = ^b[7j'{J/d) - 7j\J,/5)]J,,, 
cos6> cos6>« 

+ easier terms. 

\dj,Fj ■ J,,,| < sup |r?"([7J* + (1 - l)J]/S)\ \J\ | J*,e| 

" 0<7<1 

+ e/3C(| J| + l^iDI J*,£| + easier terms 
< Ci(| J| + |0|)| J«,e| + easier terms, 
supe-'^*|aj.Fj • J,,,| < Ci sup[e-('^+^)*(| J| + |^|)] sup[e^*i J*,e|] 

t<0 t<0 t<0 

+ easier terms, 

where sup«o e"* | J*,e | can be bounded when e is sufficiently small for any fixed > 0, 
through a routine estimate on Equations (5.25) and (5.22) for (J*(f), ^*(t)). Other 



48 



5. INVARIANT MANIFOLDS 



terms involving m*^c can bo estimated similarly. Thus, the || \\(y,n norm of terms 
involving has to be bounded by || ||cr+p,n norms. This leads to the standard 
rate condition for the regularity of invariant manifolds. That is, the regularity is 
controlled by the spectral gap. The || \\a^n norm of the term involving A^^ has to 
be bounded by || ||o-,n+2 norms. This is a new phenomenon caused by the singular 
perturbation. This problem is resolved by virtue of a special property of the fixed 
point 5* of r. Notice that if 1T2 > ci, n2 > ni, then 0^2,112 C Go-i,ni- Thus by 
the uniqueness of the fixed point, if is the fixed point of F in Gcr2,n2i 9* is also 
the fixed point of F in Gai,ni- Since exists in Ga;n for an fixed n > 1 and 

a G ij^, ^ lOP) where D is small enough, 

\\nt\\a,n<Ci+C2\\g4a+C,n, 

\\T^\\cr,n < C'3||.g*||cr,„+2 + C4 || 5* || ct+S,™ + C5, 

where Cj (1 < j < 5) depend upon ||g*(0)||„ and ||.9*(0)||„+2- Let 

M = 2{\\n+\\,,n + \\na,n), 

and F' denote the linear map defined by the right hand sides of (5.33) and (5.34). 
Since the terms duF^ , d^+F^, i9„F, and d^+F all have small || ||„ norms, F' is 
a contraction map on B{M) C L{TZ,G„^n), where B{M) is the ball of radius M. 
Thus F' has a unique fixed point g^^. Next one needs to show that ^ is indeed 
the partial derivative of with respect to e. That is, one needs to show 

(5 41) lim llg4£ + Ae)-ff.(e)-g.,eAe||^,„ ^ ^ 

^ ■ ' Ae^O Ae 

This has to be accomplished directly from Equations (5.31)-(5.32), (5.33)-(5.34) 
satisfied by g^, and g^ ^. The most troublesome estimate is still the one involving 
A^. First, notice the fact that is holomorphic in e when e > 0, and not 
differentiable at e = 0. Then, notice that g^ e G^.n for any n > 1, thus, e^^^g^ is 
differentiable, up to certain order m, in e at e = from the right, i.e. 

(rf+/rfe)"e^^^5*|.=o 

exists in H". Let 

z(i, Ae) = e(^+^^)*^'3* - e^*^'^* - iAe)tdy^'g» 
= e^*^'w;(Ae), 
where f > 0, Ae > 0, and 

w(Ae)=e(^^)*^^g,-5,-(Ae)<a25.. 
Since w{0) = 0, by the Mean Value Theorem, one has 

imAe)||„ = |mAe)-w(0)||„ < sup || -^(AAe)||„|Ae|, 

o<A<i "^e 



where at A = 0, 77^7 = 77x7 j and 

-.t[e^^^^*''dlg.-dlg4. 



dAe dAe 

dw 
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Since ^^(0) = 0, by the Mean Value Theorem again, one has 

11 dw . ,,1 ,1 dw . , dw ,, d^w 

ll^(AAe)|U = ll^(AAe) - — (0)|U < jup J|^(. 

where 



d'w ^i2fg(A.)taJg4^^] 



dA. 

Therefore, one has the estimate 

(5.42) \\z{t,Ae)\\n<\M''t^\\94n+4 ■ 

This estimate is sufficient for handling the estimate involving A^. The estimate 
involving u^,^ can be handled in a similar manner. For instance, let 

z{t, Ae) = F{u4t, e + Ae)) - e)) - Aedu,F ■ u*,^ , 

then 

\\z{t,Ae)\\a,n< |Ae|2 sup • S^^F • w*,, + • M*,«](AAe)||a,„ • 

0<A<1 

Prom the expression of F (5.28), one has 

< C'i||s'*||a+25,n[(supe''*|u,,e|)^ +supe^''*|u*,„|], 

t<0 t<0 

and the term [ ] on the right hand side can be easily shown to be bounded. In 
conclusion, let 

h = g^{e + Ae) - g^{e) - S(*,eAe, 

one has the estimate 

\\h\\,,n<ii\\h\\^,n + \Ae\^C{\\g4 <7,n+4:', ||5*||(T+2j>,n)) 
where k is small, thus 

||/i|U,„ < 2|Ae|2C(||5*|U,„+4; ||5*IU+2^,„). 

This implies that 

llml!ifei = 0, 

Ae^O |Ae| 

which is (5.4f). 

= 0. 



Let g*{t) = {^i{t),u{t)). First, let me comment on 4" 



Prom (5.32), one has 



C+(0)=0,e=0 



du dg. 



"9^+(0)'"^'"-'''"9^+(0)""'"' 
where by letting ^]^(0) ^ and e ^ 0+, ki ^ 0. Thus 

=0. 

4+(0)=0,e=0 



ae+(o) 

I shall also comment on "exponential decay" property. Since + <r, 



\\9*mn<re—\ yt<0. 
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Definition 5.10. Let = {t) , u{t)) , where 




J — oo 



depends upon ^^'(0). Thus 



u°:ei+(0)-«(0), 



defines a curve, which we call an unstable fiber denoted by where p = ( Jq, ^o) 
is the base point, £,i{0) G [— r, r] x [— r, r]. 



Let denote the evolution operator of (5.27)-(5.28), then 



That is, {^p}peA(4S) is an invariant family of unstable fibers. The proof of the 
Unstable Fiber Theorem is finished. Q.E.D. 

Remark 5.11. If one replaces the base orbit {J^{t),6*{t)) by a general orbit 

for which only j| ||„ norm is boimdcd, then the estimate (5.40) will not be possible. 
The II ||<T,n+2 norm of the fixed point 5* will not be bounded either. In such case, 
5, may not be smooth in e due to the singular perturbation. 

Remark 5.12. Smoothness of 5* in e at e = is a key point in locating 

homoclinic orbits as discussed in later sections. From intcgrable theory, information 
is known at e = 0. This key point will link "e = 0" information to "e ^ 0" studies. 
Only continuity in e at e = is not enough for the study. The beauty of the entire 
theory is reflected by the fact that although e/'^^ is not holomorphic at e = 0, e'^=: 
can be smooth at e = from the right, up to certain order depending upon the 
regularity of . This is the beauty of the singular perturbation. 



Here we give the proof of the center-stable manifold theorem 5.4, proofs of 
other invariant manifold theorems in this chapter are easier. 



5.4.1. Existence of the Center-Stable Manifold. We start with Equa- 
tions (5.20)-(5.24), let 



En{r) is of codimension 1 in the entire phase space coordinatized by (^^, J, 6, h,^^ ). 
Let X e C°°{R,R) be a "cut-ofF' function: 



5.4. Proof of the Center-Stable Manifold Theorem 5.4 



(5.43) 






We apply the cut-off 
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to Equations (5.20)-(5.24), so that the equations in a tubular neighborhood of the 
circle S^j (5.18) arc unchanged, and linear outside a bigger tubular neighborhood. 
The modified equations take the form: 



(5.45) 
(5.46) 

where A is given in (5.28), 



it = 

Vt=Av + F, 



F^{Fj,F0,Fh,F^f, T = transpose, 
Fj =XS e ^ 

sin t 



-2a{j + w^) + 2^^7+7^ cos e + n. 



Fe =XS 



7^« 



- w 



Fh = xs[Veh + Af], 

Equations (5.45)-(5.46) can be written in the equivalent integral equation form: 



J to 

v{t) = e^(*-*°)w(to) + /* e^(*-^)F(T)dT. 

J to 

We introduce the following space: For a G (^j^, ^) , and n > 1, let 



(5.47) 
(5.48) 



G, 



in iJ" norm, \\g\\a,n = 



t e [0, oo), g{t) is continuous in t 



supe-nknOI + II^^WIIn] <oo 

t>0 



Ga,n is a Banach space under the norm 
tubular neighborhood of S^; (5.18): 



Let Aa,n{'>') denote the closed 



supe--*[|^+(t)| + ||i;(t)|U]<r 

t>o 



where v is defined in (5.43). If g{t) G Aa,n{'''), r < oo, is a solution of (5.47)- (5.48), 
by letting to +oo in (5.47) and setting to = in (5.48), one has 

ft 

(5.49) ^+{t) = / e''^(*-^)F+(T)dr, 

ft 



(5.50) 



v{t) = e^'v{Q) + [ e^^'-^W{T)dT. 
Jo 



For any g{t) G Acr^nir), let T{g) be the map defined by the right hand side of (5.49)- 
(5.50). In contrast to the map T defined in (5.31)-(5.32), F contains constant terms 
of order 0(e), e.g. Fj and Fg both contain such terms. Also, Aa.nir) is a tubular 
neighborhood of the circle S^^ (5.18) instead of the ball Ba,n{r) for F. Fortunately, 
these facts will not create any difficulty in showing F is a contraction on Aa,n{'f)- 
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For any n > 1 and e < 6^, and 6 and r are small enough, and F are Lipschitz in 
g with small Lipschitz constants. F has a unique fixed point 5* in A„^n{r), following 
from standard arguments. 

5.4.2. Regularity of the Center-Stable Manifold in e. For the regularity 

of with respect to (e, a, /?, a;, w(0)), the most difficult one is of course with respect 
to e due to the singular perturbation. Formally differentiating 5* in (5.49)-(5.50) 
with respect to e, one gets 

(5.51) ^+eit) = f e""^*""^ htFt ■ ^te + dvFt ■ V,] {T)dT + Rtit), 

J+00 '- -' 

(5.52) v{t) = f e^(*-^) k+F • + dyP ■ v,] {T)dT + R{t), 

Jo L ^1 ' J 

where 

(5.53) R+{t)= f iJ,l^{t-T)e''^^*-''^F+{T)dT+ f e''^^*-''^deF+{T)dT, 

J +00 J +00 

(5.54) R{t) = tA,e^*v{Q) + f {t- T)A,e^^*-^^ F{T)dT + f e'^^'-^^d,F{T)dT, 

Jo Jo 

and /Li^j and are given in (5.37)-(5.39). The troublesome terms are the ones 

containing in (5.54). These terms can be handled in the same way as in the 
Proof of the Unstable Fiber Theorem. The crucial fact utilized is that if 11(0) G H^^ , 
then is the unique fixed point of F in both Gcr,ni and Ga^n2 for any n2 < ni. 

Remark 5.13. In the Proof of the Unstable Fiber Theorem, the arbitrary initial 
datum in (5.31)-(5.32) is ^+(0) which is a scalar. Here the arbitrary initial datum 
in (5.49)-(5.50) is v{0) which is a function of x. If v{0) E i/"^ but not H"^ for some 
ni > n2, then ^ Ga.nn in contrast to the case of (5.31)-(5.32) where g^^ G Ga^n 
for any fixed n > 1. The center-stable manifold W^" stated in the Center-Stable 
Manifold Theorem will be defined through w(0). This already illustrates why W^^ 
has the regularity in e as stated in the theorem. 

We have 

\\R\\a,n<C2rgAU,n+2 + Cz, 

for e Aa.n+2{T), whcic Gj (j = 1,2,3) are constants depending in particular 
upon the cut-off in F^ and F. Let P' denote the linear map defined by the right 
hand sides of (5.51)-(5.52). If v{0) G and e Aa,n+2{r), standard argument 

shows that P' is a contraction map on a closed ball in L{R,G„^n)- Thus P' has 
a unique fixed point g^,^^. Furthermore, if w(0) S and 9* G Aa.n+Aii^), one 

has that g^^^ is indeed the derivative of ^* in e, following the same argument as 
in the Proof of the Unstable Fiber Theorem. Here one may be able to replace 
the requirement v{0) E jj^l+a .-^-^^]^ ^ _4^^^^(7.) by just ?;(0) G and 5* G 

„-l_2(r). But we are not interested in sharper results, and the current result is 
sufficient for our purpose. 
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Definition 5.14. For any v{0) e En{r) where r is sufficiently small and En{r) 
is defined in (5.44), let g*{t) = (Ci^(i)) ^^(i)) be the fixed point off in Gcr,n, where 
one has 

^+(0) = /" e''^(*--)F+(r)dr, 

J+oo 

which depend upon v{0). Thus 

e:«(0)^^i+(0), 

defines a codimension 1 surface, which we call center-stable manifold denoted by 

^'^ n ■ 

The regularity of the fixed point 5, immediately imphes the regularity of W°*. 
We have sketched the proof of the most difficult regularity, i.e. with respect to e. 
Uniform boimdcdncss of in v{0) E En-^-4{r) and e G [0, eo), is obvious. Other 
parts of the detailed proof is completely standard. We have that is a locally 
invariant submanifold which is in (a, /3, w). W^* is in e at point in the subset 

Remark 5.15. Let S* denote the evolution operator of the perturbed nonlinear 
Schrodinger equation (5.2). The proofs of the Unstable Fiber Theorem and the 
Center-Stable Manifold Theorem also imply the following: S* is a map on H"- 
for any fixed t > 0. n > 1. 5* is also in (a,/3, w). 5* is in e as a map from 
to H" for any fixed n > 1, e € [0, eo), eo > 0. 

5.5. Perturbed Davey-Stewartson II Equations 

Invariant manifold results in Sections 5.2 and 5.1 also hold for perturbed Davey- 
Stewartson 11 equations [132], 



iqt = Tq + 



2{\q\' -u')+u^ 



Au = -4dy\q\^ , 

where g is a complex- valued function of the three variables {t,x,y), u is a real- 
valued function of the three variables (t,x,y), T = dxx — dyy, A = dxx + dyy, 
w > is a constant, and / is the perturbation. We also consider periodic boundary 
conditions. 

Under singular perturbation 

f ^ Aq-aq + (3 , 

where a > 0, /3 > are constants, Theorems 5.3 and 5.4 hold for the perturbed 
Davey-Stewartson II equations [132]. 

When the singular perturbation A is mollified into a bounded Fourier multiplier 

Aq = — Ph\k\^Qk cos kix cos k2y , 

in the case of periods {2n, 27r), 

Pk = 1, \k\ <N, pk = \k\-^, \k\ > N, 

for some large A'', jfcj^ = kf + fc|. Theorems 5.1 and 5.2 hold for the perturbed 
Davey-Stewartson II equations [132]. 
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5.6. General Overview 

For discrete systems, i.e., the flow is given by a map, it is more convenient 
to use Hadamard's method to prove invariant manifold and fiber theorems [84]. 
Even for continuous systems, Hadamard's method was often utilized [64]. On the 
other hand. Perron's method provides shorter proofs. It involves manipulation of 
integral equations. This method should be a favorite of analysts. Hadamard's 
method deals with graph transform. The proof is often much longer, with clear 
geometric intuitions. It should be a favorite of geometers. For finite-dimensional 
continuous systems, N. Fenichel proved persistence of invariant manifolds under 

perturbations of flow in a very general setting [64] . He then went on to prove 
the fiber theorems in [65] [66] also in this general setting. Finally, he applied 
this general machinery to a general system of ordinary differential equations [67]. 
As a result, Theorems 5.1 and 5.2 hold for the following perturbed discrete cubic 
nonlinear Schrodinger equations [138], 

1 



(5.55) +ie 



Qn+l — '^Qn + Qn-l 



Qnl {qn+l +qn-l) - 2w g„ 



where i = V— 1, ^Zn's are complex variables, 

qn+N = q-n, (periodic condition); and q-n = qn, (even condition); 
h = jf, and 

iVtan^ < w < TVtan— , for TV > 3, 
N N 

TT 

3 tan — < Lu < oo, for = 3. 

e £ [0, ei), a (> 0), f3 (> 0) are constants. 
This is a 2(M + 1) dimensional system, where 

M = N/2, (iVeven); and M = (A^ - l)/2, (iV odd). 

This system is a finite-difference discretization of the perturbed NLS (5.2). 

For a general system of ordinary differential equations, Kellcy [103] used the 
Perron's method to give a very short proof of the classical unstable, stable, and cen- 
ter manifold theorem. This paper is a good starting point of reading upon Perron's 
method. In the book [84], Hadamard's method is mainly employed. This book is 
an excellent starting point for a comprehensive reading on invariant manifolds. 

There have been more and more invariant manifold results for infinite dimen- 
sional systems [140]. For the employment of Perron's method, we refer the readers 
to [42]. For the employment of Hadamard's method, we refer the readers to [18] 
[19] [20] which are terribly long papers. 



CHAPTER 6 



Homoclinic Orbits 

In terms of proving the existence of a homoclinic orbit, the most common tool 
is the so-called Melnikov integral method [158] [12]. This method was subsequently 
developed by Holmes and Marsden [76], and most recently by Wiggins [201]. For 
partial differential equations, this method was mainly developed by Li et al. [136] 
[124] [138] [128]. 

There are two derivations for the Melnikov integrals. One is the so-called geo- 
metric argument [76] [201] [138] [136] [124]. The other is the so-called Liapunov- 
Schmitt argument [41] [40]. The Liapunov-Schmitt argument is a fixed-point type 
argument which directly leads to the existence of a homoclinic orbit. The condition 
for the existence of a fixed point is the Melnikov integral. The geometric argument 
is a signed distance argument which applies to more general situations than the 
Liapunov-Schmitt argument. It turns out that the geometric argument is a much 
more powerful machinary than the Liapunov-Schmitt argument. In particular, the 
geometric argument can handle geometric singular perturbation problems. I shall 
also mention an interesting derivation in [12]. 

In establishing the existence of homoclinic orbits in high dimensions, one often 
needs other tools besides the Melnikov analysis. For example, when studying or- 
bits homoclinic to fixed points created through resonances in (n > 4)-dimensional 
near-integrable systems, one often needs tools like Fenichel fibers, as presented in 
previous chapter, to set up geometric measurements for locating such homoclinci 
orbits. Such homoclinic orbits often have a geometric singular perturbation nature. 
In such cases, the Liapunov-Schmitt argument can not be applied. For such works 
on finite dimensional systems, see for example [108] [109] [138]. For such works 
on infinite dimensional systems, see for example [136] [124]. 

6.1. Silnikov Homoclinic Orbits in NLS Under Regular Perturbations 

We continue from section 5.1 and consider the regularly perturbed nonlinear 
Schrodinger (NLS) equation (5.1). The following theorem was proved in [136]. 

Theorem 6.1. There exists a cq > 0, such that for any e G (0, eg), there exists 
a codimension 1 surface in the external parameter space {a, (3, lo) € M+ x R+ x M+ 
where uj G (5- 1); and auo < [3. For any (a, /3, w) on the codimension 1 surface, the 
regularly perturbed nonlinear Schrodinger equation (5.1) possesses a symmetric pair 
of Silnikov homoclinic orbits asymptotic to a saddle Qe ■ The codimension 1 surface 
has the approximate representation given by a = 1/k{(jj), where k{co) is plotted in 
Figure 1. 

The proof of this theorem is easier than that given in later sections. To prove 
the theorem, one starts from the invariant plane 
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Figure 1. The graph of ^(a;). 

On n, there is a saddle = \fle}^ to which the symmetric pair of Silnikov homo- 
clinic orbits will be asymptotic to, where 

(6.1) i = j^-el-^^r^M^^... ^ cos^ = ^, ^e(o,|). 

Its eigenvalues are 

(6.2) A± = -e[a + i^n^\ ± 2^(^ + a;^ - /)(3/ - - !^) , 

where n = 0, 1, 2, • • • , w G (|, 1), ^„ = 1 when n < N, = 8n~^ when n > N, 
for some fixed large N, and / is given in (6.1). The crucial points to notice are: 
(1). only Aq and have positive real parts, Re{AQ } < Re{A^}; (2). all the other 
eigenvalues have negative real parts among which the absolute value of Re{Aj} — 
Re{A2 } is the smallest; (3). |Re{A2 }| < RejAg }. Actually, items (2) and (3) are 
the main characteristics of Silnikov homoclinic orbits. 

The unstable manifold VF"(Qe) of has a fiber representation given by Theo- 
rem 5.2. The Melnikov measurement measures the signed distance between W^"((5e) 
and the center-stable manifold W^^ proved in Thoerem 5.1. By virtue of the Fiber 
Theorem 5.2, one can show that, to the leading order in e, the signed distance is 
given by the Melnikov integral 

/+00 /'27r 
/ [dgF,{qomdlqo{t) - aqo{t) + /3) 
-oo Jo 

+dgFi{qo{t)){d^'^ - + p)]dxdt, 

where qo{t) is given in section 3.1, equation (3.14); and dqFi and dgFi are given in 
section 4.1, equation (4.15). The zero of the signed distance implies the existence 
of an orbit in W"{Qe) n W^''. The stable manifold W'^{Qe) of Qe is a codimension 
1 submanifold in W^^. To locate a homoclinic orbit, one needs to set up a second 
measurement measuring the signed distance between the orbit in W"(Qe) H W^^* 
and ^{Qe) inside W^'^. To set up this signed distance, first one can rather easily 
track the (perturbed) orbit by an unperturbed orbit to an 0{e) neighborhood of 
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n, then one needs to prove the size of W^{Qe) to be 0{e'^) [v < 1) with normal 
form transform. To the leading order in e, the zero of the second signed distance is 
given by 

aLo{A'y) 
2sm 

where A7 = —Ado and i^q is given in (3.14). To the leading order in e, the common 
zero of the two second signed distances satisfies a = 1/k{u)), where k{u)) is plotted 
in Figure 1. Then the claim of the theorem is proved by virtue of the implicit 
function theorem. For rigorous details, see [136]. In the singular perturbation case 
as discussed in next section, the rigorous details are given in later sections. 

6.2. Silnikov Homoclinic Orbits in NLS Under Singular Perturbations 

We continue from section 5.2 and consider the singularly perturbed nonlinear 
Schrodinger (NLS) equation (5.2). The following theorem was proved in [124]. 

Theorem 6.2. There exists a€(j> 0, such that for any e S (0,69)7 there exists 
a codimension 1 surface in the external parameter space (a, /?, lu) G M+ x M+ x 
M"*" where uj € {^,1)/S, S is a finite subset, and aoj < j3. For any (a,/?,^) on 
the codim,ension 1 surface, the singularly perturbed nonlinear Schrodinger equation 

(5.2) possesses a symmetric pair of Silnikov homoclinic orbits asymptotic to a saddle 
Qt- The codimension 1 surface has the approximate representation given by a = 
1/k{u), where k{uj) is plotted in Figure 1. 

In this singular perturbation case, the persistence and fiber theorems are given 
in section 5.2, Theorems 5.4 and 5.3. The normal form transform for proving the 
size estimate of the stable manifold W^{Qe) is still achievable. The proof of the 
theorem is also completed through two measurements: the Melnikov measurement 
and the second measurement. 

6.3. The Melnikov Measurement 
6.3.1. Dynamics on an Invariant Plane. The 2D subspace H, 

(6.3) n = {g| d^q = 0}, 

is an invariant plane under the flow (5.2). The governing equation in 11 is 

(6.4) iq = 2[\q\'^ -LO^]q + i€[-aq + (3], 

where • = ^ • Dynamics of this equation is shown in Figure 2. Interesting dynamics 
is created through resonance in the neighborhood of the circle S^j'- 

(6.5) S^ = {q€n\ \q\=io}. 

When e = Q, Sui consists of fixed points. To explore the dynamics in this neighbor- 
hood better, one can make a series of changes of coordinates. Let q = \//e*^, then 
(6.4) can be rewritten as 

(6.6) 7 = e(-2a/+2/3V7cos6') , 

(6.7) e = -2{I - oj'') - ed^'"'^ 



There are three fixed points: 
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(1) The focus Oe in the neighborhood of the origin, 



(6.8) 



cose = ^, ^e(o,f). 



Its eigenvalues are 



sin 9 — ea, 

where I and 6 are given in (6.8). 
(2) The focus in the neighborhood of S^j (6.5), 



(6.10) 

Its eigenvalues are 



cos9=^, ^e(-f,o) 



(6.11) Mi,2 = ±«A/^W-4V7/3sin0 + e(^^^^ 



2 

— ea. 



where / and 9 are given in (6.10). 
(3) The saddle in the neighborhood of (6.5), 



(6.12) 



lcos^ = ^, ^e(o,f). 

Its eigenvalues are 



2 

— ea. 



(6.13) Mi,2 =±V^y4V7/3sin0-e(^^^^ 

where / and 6 arc given in (6.12). 
Now focus our attention to order ^ neighborhood of S^i (6.5) and let 
J = I — ui^, J = y/ej, T = y/et, 

we have 



(6.14) / = 2 



' cost 



(6.15) 9' = -2j - 



smf 



where ' = ^ . To leading order, wc get 

(6.16) / = 2[-aLo^ + ^Lo cos 9] , 

(6.17) 9' = -2j . 

There are two fixed points which are the counterparts of and Qe (6.10) and 
(6.12): 

(1) The center P,, 

(6.18) j = 0, cos9=^, ^e(-|,o). 

Its eigenvalues are 

(6.19) /ii,2 = ±^2^^;(/3^ - aV)i 
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Figure 2. Dynamics on the invariant plane H. 



(2) The saddle Q,, 
(6.20) j = 0, cos( 

Its eigenvalues are 



(^) 



(6.21) 



Mi,2 = ±2V^J(/32 - a 



'2 "V)i 



In fact, (6.16) and (6.17) form a Hamiltonian system with the Hamiltonian 
(6.22) f + 2uj{-auje + (3sme). 

Connecting to is a fish-like singular level set of 7i, which intersects the axis 
j = at and Q = (0,^), 



(6.23) 



at^(6'-6l,) =/3(sin6l-sin6l,), ( — ^,0), 



where 6** is given in (6.20). See Figure 3 for an illustration of the dynamics of 
(6.14)-(6.17). For later use, we define a piece of each of the stable and unstable 
manifolds of Q.^, 

for some small S > Q, and 



(6.24) 



-^y2uj[auj{0 - 61,) - P{sin0- sin 61*) 
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Figure 3. The fish-like dynamics in the neighborhood of the res- 
onant circle S,^. 



<j)^{9) and cj)l{0) perturb smoothly in 9 and -^e into and 0^ for (6.14) and 
(6.15). 

The homoclinic orbit to be located will take off from along its unstable 
curve, flies away from and returns to 11, lands near the stable curve of and 
approaches Qe spirally. 

6.3.2. A Signed Distance. Let p be any point on 0^ (6.24) which is the 
unstable curve of in 11 (6.3). Let (7e(0) and (7o(0) be any two points on the 
unstable fibers |e and with the same coordinate. See Figure 4 for 

an illustration. By the Unstable Fiber Theorem 5.3, J^^ is in e for e € [0,eo), 
eo > 0, thus 

|lge(0)-go(0)||„+8<Ce. 
The key point here is that C if* for any fixed s > 1. By Remark 5.15, the 
evolution operator of the perturbed NLS equation (5.2) 5* is in e as a map from 
to i?" for any fixed 7i > 1, e e [0, eo), eo > 0. Also is a map on i?" for 
any fixed i > 0, n > 1. Thus 

\MT) - (7o(r)||n+4 = \\s'^{qM) - 5^(go(0))|l„+4 < C,e, 
where T > is large enough so that 

<zo(r) e w^l^ |,.o . 

Our goal is to determine when q^iT) € W^" through Melnikov measurement. Let 
qt{T) and qo{T) have the coordinate expressions 

(6.25) q,(T)^{it'',Ve), qo{T)^i^t'",v,). 

Let qe(T) be the unique point on W.^'^^, which has the same u-coordinate as qe{T), 

q,{T) = {i+-\v,)eW:U- 
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Figure 4. The Melnikov measurement. 

By the Center-Stable Manifold Theorem, at points in the subset W^''_^^, W^* is 
smooth in e for e € [0, eo), eo > 0, thus 

(6.26) ||g,(r)-g,(T)||„ < Cse. 

Also our goal now is to determine when the signed distance 

is zero through Melnikov measurement. Equivalently, one can define the signed 
distance 

di-(VFi(go(r)),<7,(T)-g,(r)) 

= a,Fi(go(T))fe(T) - UT)) + 9,-Fi(<zo(T))fe(T) - 9^(7))", 

where VFi is given in (4.17), and qo{t) is the homoclinic orbit given in (3.14). In 
fact, qe{t), qe{t), qo{t) € H^, for any fixed n > 1. The rest of the derivation for 
Melnikov integrals is completely standard. For details, see [136] [138]. 

(6.27) di = eMi +o(e), 
where 

/ + 00 p2-K 
/ [dgF,{qf,{t)){dlqo{t) ^ aqoit) + P) 
-oo Jo 

+ dqFi{qo{t)){dl'^ - + f3)]dxdt, 

where again qo{t) is given in (3.14), and dqFi, and dqFi are given in (4.17). 

6.4. The Second Measurement 

6.4.1. The Size of the Stable Manifold of Q,. Assume that the Melnikov 
measurement is successful, that is, the orbit q^it) is in the intersection of the unsta- 
ble manifold of and the center-stable manifold W'^{Qe) n W^^. The goal of the 
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second measurement is to determine when qt{t) is also in the codimension 2 stable 
manifold of Qe, W^{Qe), in The existence of W^{Qe) follows from standard 
stable manifold theorem. W^{Qe) can be visualized as a codimension 1 wall in 
with base curve 0^ (6.24) in 11. 

Theorem 6.3. [124] The size of W^{Q^) off U is of order 0{^) for to G 
(^, l) /S, where S is a finite subset. 

Starting from the system (5.7)-(5.9), one can only get the size of W^{Qf) off 
n to be 0{e) from the standard stable manifold theorem. This is not enough for 
the second measurement. An estimate of order 0(6"), k < 1 can be achieved 
if the quadratic term Af2 (5.14) in (5.9) can be removed through a normal form 
transformation. Such a normal form transformation has been achieved [124], see 
also the later section on Normal Form Transforms. 

6.4.2. An Estimate. From the explicit expression of qo{t) (3.14), we know 

that qo{t) approaches 11 at the rate ©(e"^^""""). Thus 

(6.28) distance |go(r + ^|lne|),n| < Ce 

Lemma 6.4. For all t e 



|2 



r,T+i|lne| 

(6.29) \\qe{t) - qomn < C\e\lne\ 
where Ci = C'i(T). 

Proof: We start with the system (5.49)-(5.50). Let 

q.it) = {^t''it),J%t),e^it),h%t),cr'it)), 

Qo{t) = iit'°{t),J%t),9%t),h%t),i^''{t)). 
Let Ti(> r) be a time such that 

(6.30) \\q4t)-qnit)\\n<C2e\ln€f, 

for all t e [T,Ti], where C2 = C2{T) is independent of e. From (6.26), such a 
Ti exists. The proof will be completed through a continuation argument. For 

ie [r,Ti], 

(ki+'°(i)l + ICr'"(i)l) + Wh^'iMn < C,re-i^(''^\ \J^it)\ < C^Vc 

(6.31) I J^(t)| < |j"(t)| + \.rit) - J°(t)| < |J"(0| + C2e\ lne|2 < C,V~e, 
{\^t'\t)\ + ler-'WI) + W{t)\U < Csre-i^^'-^^+C2e\\ne\^ 

where r is small. Since actually qdt), qo{t) G for any fixed n > 1, by Theo- 
rem 5.4, 

(6.32) \^t'%t) - it'°{t)\ < Ce\\ve{t) - vo{t)\\n + Cre, 

whenever Ve{t), vo{t) G En+4{r), where Ve{T) = and vo{T) = vq are defined in 
(6.25). Thus we only need to estimate \\ve{t) — vo{t)\\n- Prom (5.50), we have for 
t e [T, Ti] that 

(6.33) v{t) = e^^*-^\{T) + j\^^*-^^F{T)dT. 
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Let Av{t) = Ve{t) - VQ{t). Then 

Av{t) = [e^(*-^) - e^l-''(*-^)]^;o(T) + e^^'-'^^ Av{T) 
(6.34) +^*e^(*--)[F(r)-F(r)|e=o]rfT 



By the condition (6.31), we have for t G [T,Ti] that 

(6.35) \\F{t) - F{t)U=o\\n < [C78x/i+C9re-5''(*-^)]e|lne|2. 
Then 

(6.36) ||Ai;(i)||„ < Cioe(t-T) + Ciire|lne|2 + Ci2V^(t-r)2e|lne|2. 

Thus by the continuation argument, for t £ [T,T + ^|lne|], there is a constant 
Ci = Ci(T), 

(6.37) \\Av{t)\\n<Cie\\ne\''. 
Q.E.D. 

By Lemma 6.4 and estimate (6.28), 

(6.38) distance I (r + ^1 In e|),n| < Celine]^. 

By Theorem 6.3, the height of the wall W^{Qe) off n is larger than the distance 
between ^^(T + ^| lne|) and 11. Thus, ii qf:{T + \ lne|) can move from one side of 
the wall W^{Qf_) to the other side, then by continuity qe{T + ^| lne|) has to be on 
the wall W^(Qe) at some values of the parameters. 

6.4.3. Another Signed Distance. Recall the fish-like singular level set given 
by H (6.22), the width of the fish is of order 0{y/e), and the length of the fish is of 
order 0(1). Notice also that qo(t) has a phase shift 

(6.39) e°i=e°{T+-\lne\)-e^{0). 

M 

For fixed (3, changing a can induce 0{1) change in the length of the fish, 0{^) 
change in 6*?, and £»(!) change in 0°{O). See Fi gurc 5 for an illustration. The 
leading order signed distance from ge(T + i| lne|) to W^{Qc) can be defined as 

d = nuo, - nuo, o'{o) + ^?) 

^ ' ' =2w [awOl + /J[sin^O(0) - sin(^0(0) + e\)]] , 

where H is given in (6.22). The common zero of Mi (6.27) and d and the implicit 

function theorem imply the existence of a homoclinic orbit asymptotic to Qf . The 
common roots to Mi and d are given by a = 1/K(a;), where is plotted in 

Figure 1. 
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/ 



WXQ) 





7 




0(e) 



Figure 5. The second measurement. 



6.5. Silnikov Homoclinic Orbits in Vector NLS Under Perturbations 

In recent years, novel results have been obtained on the solutions of the vector 
nonlinear Schrodinger equations [3] [4] [206] . Abundant ordinary integrable results 
have been carried through [203] [70], including linear stability calculations [69]. 
Specifically, the vector nonlinear Schrodinger equations can be written as 



where p and q are complex valued functions of the two real variables t and x, and x is 
a positive constant. These equations describe the evolution of two orthogonal pulse 
envelopes in birefringent optical fibers [159] [160], with industrial applications 
in fiber communication systems [81] and all-optical switching devices [93]. For 
linearly birefringent fibers [159], x = 2/3. For elliptically birefringent fibers, x 
can take other positive values [160]. When x = 1; these equations are first shown 
to be integrable by S. Manakov [150], and thus called Manakov equations. When 
X is not 1 or 0, these equations are non-integrable. Propelled by the industrial 
applications, extensive mathematical studies on the vector nonlinear Schrodinger 
equations have been conducted. Like the scalar nonlinear Schrodinger equation, 
the vector nonlinear Schrodinger equations also possess figure eight structures in 
their phase space. Consider the singularly perturbed vector nonlinear Schrodinger 
equations. 



where x) and q{t, x) are subject to periodic boundary condition of period 27r, 
and are even in x, i.e. 



(6.42) 



(6.41) 




p(t, X + 27r) — p{t, x) , p{t, —x) — p{t, x) , 
q{t, x + 2tt) ^ q{t, x) , q{t, -x) = q{t, x) , 
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we (1, 2), a > and j3 are real constants, and e > is the perturbation parameter. 
We have 

Theorem 6.5 ([135]). There exists a eo > 0, suc:h that for any e € (0, eo), 

there exists a codimension 1 surface in the space of {a, (i, oj) G K"*" x R"*" x K"*" 
where ui G [1,2)/ S, S is a finite subset, and auj < ^/2(3. For any {a,P,Lo) on 
the codimension 1 surface, the singularly perturbed vector nonlinear Schrddinger 
equations (6.41)-(6.42) possesses a homoclinic orbit asymptotic to a saddle Q^. 
This orbit is also the homoclinic orbit for the singularly perturbed scalar nonlinear 
Schrddinger equation studied in last section, and is the only one asymptotic to 
the saddle Qf for the singularly perturbed vector nonlinear Schrddinger equations 
(6.41)-(6.42). The codimension 1 surface has the approximate representation given 
by a = 1/k{w), where k{w) is plotted in Figure 1. 

6.6. Silnikov Homoclinic Orbits in Discrete NLS Under Perturbations 

We continue from section 5.6 and consider the perturbed discrete nonlinear 
Schrodinger equation (5.55). The following theorem was proved in [138]. 
Denote by Sjv {N >7) the external parameter space, 



OJ e (ATtan— ,A/^tan— ), 



Sjv = |(a;,a,/3) 

a€ (0,ao),/3G (0,/3o); 

where ao and Po are any fixed positive numbers. 

Theorem 6.6. For any N (7 < N < oo), there exists a positive number €q, 
such that for any e G (0, eg), there exists a codimension 1 surface in Ejv/ for any 
external parameters (u),a,(3) on E^, there exists a homoclinic orbit asymptotic to a 
saddle Qe ■ The codimension 1 surface E^ has the approximate expression a = 1 /k, 
where k = K(a;; N) is shown in Fig. 6. 

In the cases (3 < < 6), k is always negative. For N > 7, k can be positive 
as shown in Fig. 6. When N is even and > 7, there is in fact a symmetric pair of 
homoclinic orbits asymptotic to a fixed point at the same values of the external 
parameters; since for even N, we have the symmetry: If g„ = f{n,t) solves (5.55), 
then g„ = /(n + N/2, t) also solves (5.55). When N is odd and > 7, the study can 
not guarantee that two homoclinic orbits exist at the same value of the external 
parameters. 



6.7. Comments on DSII Under Perturbations 

Wc continue from section 5.5. Under both regular and singular perturbations, 
the rigorous Melnikov measurement can be established [132]. It turns out that only 
local well-posedness is necessary for rigorously setting up the Melnikov measure- 
ment, thanks to the fact that the unperturbed homoclinic orbits given in section 
3.3.1 are classical solutions. Thus the Melnikov integrals given in section 4.3.2 
indeed rigorously measure signed distances. For details, see [132]. 

The obstacle toward proving the existence of homoclinic orbits comes from a 
technical difficulty in solving a linear system to get the normal form for proving the 
size estimate of the stable manifold of the saddle. For details, see [132]. 
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Figure 6. The graph of k{uj;N). 



(6.43) 



6.8. Normal Form Transforms 

Consider the singularly perturbed Davey-Stewartson II equation in section 5.5, 

iqt = Tq+ [2(|g|2 - u;2) + Uy]q + ie[Aq - aq + P] , 
Au=-4:dy\q\^ . 
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Let 

q{t, X, y) = [p{t) + fit, X, y)]e'^W, (/) = 0, 
where ( , ) denotes spatial mean. Let 

I={\q\')=p' + {\ff), J = I-w\ 
In terms of the new variables (J, 6, /), Equation (6.43) can be rewritten as 



(6.44) 


J 


(6.45) 


e 


(6.46) 


ft 


where 





2a( J + w^) + 2/3VV + cos 6 



„ ^ ^ sin 9 
-2J-eP , =+-R2> 

ft = L,f + V,f-m2-i^fs, 



+ eR. 



-iTf + e(A - a)f - 2iL0^A-^T{f + /), 



Afo = 2uj 



A-iT|/|2 + /A-iT(/ + /) - (/A-iT(/ + /)) 



and Ri, R^, VJ, and are higher order terms. It is the quadratic term that 
blocks the size estimate of the stable manifold of the saddle [132]. Thus, our goal is 
to find a normal form transform g = f + K{f, f) where K is & bilinear form, that 
transforms the equation 

ft = LJ - iAf2, 
into an equation with a cubic nonlinearity 

gt = L,g + 0{\\g\\l), (s > 2), 
where is given in (6.46). In terms of Fourier transforms. 



where k = {ki,k2) S Z^, ^ = {kix, K2y)- The terms in Af2 can be written as 



f{k)f{-i)+fii)f{-k) 

/A-iT/ - (/A-^T/) = 1 [«(^) + a(^)]/>)/We^('=+^)-« , 
/A-iT/-(/A-iT/) = i Y 



where 



a{e)f{k)f{-t)+a{k)f{t)f{-k) 



1.2 2 _ 1.2 2 



i(k+eyi 



Q/(^k^ — ; 2 2 / 2 2 



We will search for a normal form transform of the general form, 

g = f + K{f,f), 
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where 



K{fJ) = [Ki{k,i)f{k)f{^) + K2{Kl)f{k)f{-i) 



+K2{e, k)f{-k)m + Ksik, e)f{-k)f{-e) 



where Kj{k,£), [j = 1,2,3) are the unknown coefficients to be determined, and 
Kj{k,£) = Kj{i,k), {j = 1,3). To eliminate the quadratic terms, we first need to 
set 

iL,K{f, f) - iK{L,f, /) - iK{f, LJ) = A/-2, 
which takes the explicit form: 

(ai + ia)ki{k,l) + B{£)k2{k, I) + B{k)k2{l, k) 



1 



(6.47) +B{k + l)Ks{k,£) = —[B{k)+B{i)], 

2uj 



-B{e)Ki{k, £) + ((72 + ia)K2{k, i) + B{k + e)K2{£, k) 
(6.48) +B{k)ki{k,£) = ^[B{k + £) + B{£)\, 



-B{k)Ki{k, £) + B{k + £)K2{k, I) + {a^ + i(T)K2{i, k) 
(6.49) +B{t)k^{kJ.) = ^[B{k + C) + B{k)], 



B{k + e)Ki{k, £) - B{k)K2{k,£) - B{£)K2{£, k) 
(6.50) +{(Ti + ia)k:i{k,£)=Q, 

where B{k) = 2uj'^a{k), and 

a - 2{ki£iK\ + k2£2Kl) 

CTi = 2{k2£2Kl - ki£iKl) + B{k + £)- B{k) - B{£) , 

(72 = 2[{k2 + £2)£24 - {ki + £i)eiKl] + B{k + £)- B{k) + B(£) , 

C73 = 2[(fc2 + £2)k2>^l - (fci + ^i)fciK?] + B{k + £) + B{k) - B{£) , 

CT4 = 2[{kl + fc2^2 + £l)4 - {kj + ki£i + £1)kI] + B{k + £) + B{k) + B{£) . 

Since these coefficients are even in (/c,^), we will search for even solutions, i.e. 

kj{-k,-£) = kj{k,£), j = 1,2,3. 

The technical difficulty in the normal form transform comes from not being able 
to answer the following two questions in solving the linear system (6.47)-(6.50): 

(1) Is it true that for all k,£ G Z^/{0}, there is a solution ? 

(2) What is the asymptotic behavior of the solution as k and/or £ — > oo ? In 
particular, is the asymptotic behavior like and/or £~"^ (m > 0) ? 

Setting dy = in (6.43), the singularly perturbed Davey-Stewartson II equation 
reduces to the singularly perturbed nonlinear Schrodinger (NLS) equation (5.2) in 
section 5.2. Then the above two questions can be answered [124]. For the regularly 
perturbed nonlinear Schrodinger (NLS) equation (5.1) in section 5.1, the answer to 
the above two questions is even simplier [136] [124]. 
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xdl^ll) = 



6.9. Transversal Homoclinic Orbits in a Periodically Perturbed SG 

Transversal homoclinic orbits in continuous systems often appear in two types 
of systems: (1). periodic systems where the Poincarc period map has a transversal 
homoclinic orbit, (2). autonomous systems where the homoclinic orbit is asymptotic 
to a hyperbolic limit cycle. 

Consider the periodically perturbed sine-Gordon (SG) equation, 

(6.51) utt = (?Uxx + siiiM + e[-aut + u^xdl^ll) cost], 
where 

1, ||u||<M, 
0, ||u||>2M, 

for M < \\u\\ < 2M, x(ll^ll) is a smooth bump function, under odd periodic bound- 
ary condition, 

u{x + 27r, t) = u{x, t), u{x, t) = —u{x, t), 
^ < < 1, a > 0, e is a small perturbation parameter. 

Theorem 6.7 ([136], [183]). There exists an interval / C R+ sueh that for 
any a € I, there exists a transversal homoclinic orbit u = ^{x,t) asymptotic to in 

6.10. Transversal Homoclinic Orbits in a Derivative NLS 

Consider the derivative nonlinear Schrodinger equation, 

(6.52) «gt = fea; + 2|g|^5-|-ie 



(^ - \q?)q + iJ-\dxq\^q 



where 9 is a complex-valued function of two real variables t and x, e > is the 
perturbation parameter, /i is a real constant, and dx is a bounded Fourier multiplier, 

K 00 

dxq = — ^ kqk sin kx , for g = qk cos kx , 

k=l k=0 

and some fixed K. Periodic boundary condition and even constraint are imposed, 
q{t, x + 2ir) = q{t, x) , q{t, -x) = q{t, x) . 

Theorem 6.8 ([125]). There exists a eo > 0, such that for any e e (0,eo), and 
\^\ > 5.8, there exist two transversal homoclinic orbits asymptotic to the limit cycle 
(?e = |exp{-i[|t + 7]}. 



CHAPTER 7 



Existence of Chaos 

The importance of homoclinic orbits with respect to chaotic dynamics was first 
realized by Poincare [177]. In 1961, Smale constructed the well-known horseshoe 
in the neighborhood of a transversal homoclinic orbit [190] [191] [192]. In partic- 
ular, Smale's theorem implies Birkhoff 's theorem on the existence of a sequence of 
structurely stable periodic orbits in the neighborhood of a transversal homoclinic 
orbit [27]. In 1984 and 1988 [172] [173], Palmer gave a beautiful proof of Smale's 
theorem using a shadowing lemma. Later, this proof was generalized to infinite 
dimensions by Steinlein and Walther [193] [194] and Henry [83]. In 1967, Silnikov 
proved Smale's theorem for autonomous systems in finite dimensions using a fixed 
point argument [186]. In 1996, Palmer proved Smale's theorem for autonomous 
systems in finite dimensions using shadowing lemma [50] [174]. In 2002, Li proved 
Smale's theorem for autonomous systems in infinite dimensions using shadowing 
lemma [125]. 

For nontransversal homoclinic orbits, the most well-known type which leads to 
the existence of Smale horseshoes is the so-called Silnikov homoclinic orbit [184] 
[185] [188] [55] [56]. Existence of Silnikov homochnic orbits and new constructions 
of Smale horseshoes for concrete nonlinear wave systems have been established in 
finite dimensions [138] [139] and in infinite dimensions [136] [120] [124]. 

7.1. Horseshoes and Chaos 

7.1.1. Equivariant Smooth Linearization. Linearization has been a pop- 
ular topics in dynamical systems. It asks the question whether or not one can 
transform a nonlinear system into its linearized system at a fixed point, in a small 
neighborhood of the fixed point. For a sample of references, see [80] [182]. Here we 
specifically consider the singularly perturbed nonlinear Schrodinger (NLS) equation 
(5.2) in sections 5.2 and 6.2. The symmetric pair of Silnikov homoclinic orbits is 
asymptotic to the saddle = \fle^^ , where 

(7.1) / = ,^2_^^y^T3^+... , cos0=^, ^e(0,|). 
Its eigenvalues are 

(7.2) A± = -e[a + ± 2^(^ + a;^ - 7)(3/ - a;^ - !^) , 

where n = 0, 1, 2, • • • , a; e (|, 1), and / is given in (7.1). In conducting linearization, 
the crucial factor is the so-calUed nonresonance conditions. 

Lemma 7.1 ([120] [129]). For any fixed e € (0, eq); let he the codimension 
1 surface in the external parameter space, on which the symmetric pair of Silnikov 
homoclinic orbits are supported (cf: Theorem 6.2). For almost every (a, /3, w) e 
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the eigenvalues (6-2) satisfy the nonresonance condition of Siegel type: There 
exists a natural number s such that for any integer n>2, 



3=1 



for all r = 2, 3, • • • ,n and all hjh,- " Jr € Z, where A„ = A+ for n > 0, 
An = Al„_i for n< 0. 

Thus, in a neighborhood of the singularly perturbed NLS (5.2) is analyti- 
cally equivalent to its linearization at [166]. In terms of eigenvector basis, (5.2) 
can be rewritten as 

X = —ax — by + Nx{Q), 
y = bx-ay+Afy(Q), 

(7.3) ii = 7i^i+A4,(g), 

Z2 = 722^2 +A42(Q), 

Q = LQ+Mq{Q); 

where a = — Re{A2'}, b = ImjAg }, 7i = Aq", 72 = A]*"; AA's vanish identically in 
a neighborhood Q of Q = 0, Q = {x,y,Zi,Z2,Q), Q is associated with the rest of 
eigenvalues, L is given as 

LQ = -iQ^^ - 2i[{2\Q,f - u^)Q + QIQ] + e[-aQ + Qfc] , 

and is given in (7.1). 

7.1.2. Conley-Moser Conditions. We continue from last section. Denote 

by /ife (fc = 1,2) the symmetric pair of Silnikov homoclinic orbits. The symmetry 
a of half spatial period shifting has the new representation in terms of the new 
coordinates 

(7.4) cr o (x, y, Zi,Z2,Q) = (x, y, zi, -Z2, aoQ). 
Definition 7.2. The Poincarc section Eo is defined by the constraints: 

2/ = 0, rycxp{— 27ra/6} < x < rj; 
< Zl < ?7, -Tj<Z2<V, WQW < m 
where 77 is a small parameter. 

The horseshoes are going to be constructed on this Poincarc section. 

Definition 7.3. The auxiliary Poincare section Si is defined by the con- 
straints: 

Zl =ri, -Tj < Z2 < rj, 



y/x^ +y^ <ri, |1Q|1<?7. 

Maps from a Poincare section to a Poincare section are induced by the flow. 
Let and be the coordinates on Sq and Si respectively, then the map Pg 
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from So to Si has the simple expression: 



y = 



x° cos 



■ 

X sm 



h 7] 
— In-Q 

71 2^1 

b 



In- 
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Let (3° and be the intersection points of the homocHnic orbit hi with Sq and 

El respectively, and let = — Q^, and = — Q\. One can define slabs 
in Eq. 

Definition 7.4. For sufficiently large natural number I, we define slab Si in 
So as follows: 



Si 



QeSo 



r?exp{-7i(io,2(;+i) - ^)} < 



z\{Q) < ?7exp{-7i(to,2i - 2^)}, 



\5r{Q)\ < ?7exp{--a to,2;}, 
1^2(^0 (<3))| <»?exp{-ia io,2;}, 
||g'(Po'(Q))ll <^exp{-lato,2a|, 



where 



'-[Iw - (fii] + o{l) 



as Z ^ +00, are the time that label the fixed points of the Poincare map P from 
So to So [120] [129], and the notations x°{Q), z^{Po{Q)), etc. denote the x° 
coordinate of the point Q, the Z2 coordinate of the point Pq{Q), etc.. 

Si is defined so that it includes two fixed points and pj' of P. Let Si,^ = croSi 
where the symmetry a is defined in (7.4). We need to define a larger slab Si such 
that SiUSi^a- C Si. 

Definition 7.5. The larger slab §1 is defined as 



Si = 



€ So 



7?exp{-7i(io,2(;+i) - ^)} < 



ZiiQ) < ??exp{-7i(to,2; - ^)}, 
|a;°(Q) - a;°| < ?7exp{-^a to,2i}, 
14(^0 (Q))l < 14*1 +^exp{-ia to,2j}, 
||Q'(Po(Q))|| <r?exp{-iaio,2/}|, 
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where z\ ^ is the z|-coordinatc of Q\. 

Definition 7.6. In the coordinate system Q^}, the stable boundary 

of Si, denoted by dgSi, is defined to be the boundary of Si along [Sp, (5'^)-directions, 
and the unstable boundary of Si, denoted by duSi, is defined to be the boimdary of 
Si along (5°, Z2)-directions. A stable slice V in Si is a subset of Si, defined as the 
region swept out through homeomorphically moving and deforming dgSi in such a 
way that the part 

dsSi n duSi 

of daSi only moves and deforms inside duSi. The new boundary obtained through 
such moving and deforming of dsSi is called the stable boundary of V, which is 

denoted by dgV . The rest of the boundary of V is called its unstable boundary, 
which is denoted by duV . An unstable slice of Si, denoted by H, is defined similarly. 

As shown in [120] and [129], under certain generic assumptions, when I is suf- 
ficiently large, P{Si) and P{Si^„) intersect Si into four disjoint stable slices {Fi, V2} 
and {V-\,V-2) in §1. Vj^s {j = 1, 2,-1, —2) do not intersect dgSi; moreover, 

(7.5) dsVi C P{dsSi),{i = 1,2); dsVi C P{dsSi,,),{i = -1,-2). 
Let 

(7.6) H,=P-\V,). (.7 = 1,2,-1,-2), 

where and for the rest of this article, P^^ denotes prcimagc^ of P. Then Hj (j = 
1, 2, —1, —2) are unstable slices. More importantly, the Conk^y-AIost^r condit ions arc 



satisfied. Specifically, Conley-Moser conditions are: Conley-Moser condition (i): 



V, = pm, 

dsVj = P(dsHj), (j = 1,2,-1,-2) 
duV, = P{duH,). 



Conley-Moser condition (ii): There exists a constant < v < 1, such that for any 
stable slice V cVj {j = 1, 2, —1, —2), the diameter decay relation 

d{V) < vd{V) 

holds, where d{-) denotes the diameter [120], and V = P{VnHk), {k = 1, 2, -1, -2); 
for any unstable slice H c Hj (j = 1, 2, — 1, — 2), the diameter decay relation 

diH) < vd{H) 

holds, where H = p-^{H n Vfc), {k = 1, 2, -1, -2). 

The Conley-Moser conditions are sufficient conditions for establishing the topo- 
logical conjugacy between the Poincare map P restricted to a Cantor set in Eg, and 
the shift automorphism on symbols. It also display the horseshoe nature of the in- 
tersection of Si and Si^„ with their images P{Si) and P{Si^cr)- 

7.1.3. Shift Automorphism. Let W be a set which consists of elements of 
the doubly infinite sequence form: 

a = (• • •a_2a_iao,aia2 • ••), 

where Ofe G {1,2,-1,-2}; k G Z. We introduce a topology in W by taking as 
neighborhood basis of 

a* = (. . ■a*_2a*_iao,ala2 ■ ••), 
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the set 



Wj = Uew 



Uk = al (|fc| < j) 



for j = 1, 2, • • •. This makes W a topological space. The shift automorphism x is 
defined on W by 

Va e W, x(<i) = ^) where bk = afe+i- 

The shift automorphism x exhibits sensitive dependence on initial conditions, which 
is a hallmark of chaos. 
Let 

a = (• • •a_2a_iao, aia2 • ••), 
be any element of W. Define inductively for fc > 2 the stable slices 

KoO_i = P{Ha_i) n 7?ao7 
Kioo_i...a_fc = -P(^_i...o_/i) I"! -f^oo- 

By Conley-Moser condition (ii), 

diVaoa-i.-.a-k) < J^l<Koa_i...a_(;,_i) ) < ••• < J^i "^(^(Koa-i ) ■ 

Then, 

OO 

^(«) = n ^aoa_,...a_. 
fe=l 

defines a 2 dimensional continuous surface in Sq; moreover, 

(7.7) 9y(a) C duSi. 
Similarly, define inductively for fc > 1 the unstable slices 

Haoai...ak ~ P (-^^oi...afc ^ ^aa)- 

By Conley-Moser condition (ii), 

d{Haoai...ak) < '^2d{Haoai...ak-i) < ■•• < t^2'l{Hao)- 

Then, 

OO 

-f^('^) = n Haoai...ak 
fe=0 

defines a codimension 2 continuous surface in Sq; moreover, 

(7.8) aif(a) C dsSi. 
By (7.7;7.8) and dimension count, 

V{a) n H{a) ^ 

consists of points. Let 

p G V{a) n ff (a) 
be any point in the intersection set. Now we define the mapping 

(t){a) = p. 
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By the above construction, 

P{p) = Hxia)). 

That is, 

Po(t) = (f) ox- 
Let 

A = 

then A is a compact Cantor subset of Si , and invariant under the Poincare map P. 
Moreover, with the topology inherited from §1 for A, ^ is a homeomorphism from 
Wto A. 

7.2. Nonlinear Schrodinger Equation Under Singular Perturbations 

Finally, one has the theorem on the existence of chaos in the singularly per- 
turbed nonlinear Schrodinger system (5.2). 

Theorem 7.7 (Chaos Theorem, [129]). Under certain generic assumptions, 
for the perturbed nonlinear Schrodinger system (5.2), there exists a compact Cantor 
subset A of Si, A consists of points, and is invariant under P. P restricted to A, 
is topologically conjugate to the shift automorphism x on four symbols 1, 2, —1, —2. 
That is, there exists a homeomorphism 

such that the following diagram commutes: 




W ^ A 



Although the symmetric pair of Silnikov homoclinic orbits is not structurally 
stable, the Smale horseshoes are structurally stable. Thus, the Cantor sets and the 
conjugacy to the shift automorphism are structurally stable. 

7.3. Nonlinear Schrodinger Equation Under Regular Perturbations 

We continue from section 6.1 and consider the regularly perturbed nonlinear 
Schrodinger (NLS) equation (5.1). Starting from the Homoclinic Orbit Theorem 
6.1, the arguments in previous sections equally apply. In fact, under the regular 
perturbation, the argument can be simplified due to the fact that the evolution 
operator is a group. Thus, the Chaos Theorem 7.7 holds for the regularly perturbed 
nonlinear Schrodinger (NLS) equation (5.1). 

7.4. Discrete Nonlinear Schrodinger Equation Under Perturbations 

We continue from section 6.6. Starting from the Homoclinic Orbit Theorem 6.6, 
the arguments in previous sections equally apply and are easier. Thus, the Chaos 
Theorem 7.7 holds for the perturbed discrete nonlinear Schrodinger equation (5.55). 
In the case N >7 and is odd, we only define one slab Si . Consequently, the Poincare 
map is topologically conjugate to the shift automorphism on two symbols. 
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7.5. Numerical Simulation of Chaos 

The finite-difference discretization of both the regularly and the singularly per- 
turbed nonlinear Schrodinger equations (5.1) and (5.2) leads to the same discrete 
perturbed nonlinear Schrodinger equation (5.55). 

In the chaotic regime, typical numerical output is shown in Figure 1. Notice 
that there arc two typical profiles at a fixed time: (1). a breather type profile with 
its hump located at the center of the spatial period interval, (2). a breather type 
profile with its hump located at the boundary (wing) of the spatial period interval. 
These two types of profiles are half spatial period translate of each other. If we 
label the profiles, with their humps at the center of the spatial period interval, by 
"C" ; those profiles, with their humps at the wing of the spatial period interval, by 
"W", then 

(7.9) 'W'=CTo"C", 

where a is the symmetry group element representing half spatial period translate. 
The time series of the output in Figure 1 is a chaotic jumping between "C" and 
"W" , which we call "chaotic center-wing jumping" . By the relation (7.9) and the 
study in previous sections, we interpret the chaotic center-wing jumping as the 
numerical realization of the shift automorphism x on four symbols 1,2,-1,-2. 
Wc can make this more precise in terms of the phase space geometry. The figure- 
eight structure (3.14) projected onto the plane of the Fourier component cos a; is 
illustrated in Figure 2, and labeled by Lc and Lw From the symmetry, we know 
that 

Lc has the spatial-temporal profile realization as in Figure 3 with the hump located 
at the center, i^y corresponds to the half spatial period translate of the spatial- 
temporal profile realization as in Figure 3, with the hump located at the boundary 
(wing). An orbit inside Lc-, Lcin has a spatial-temporal profile realization as in 
Figure 4. The half period translate of Lcin-, Lwin is inside Lw- An orbit outside 
Lc and Lw, Lout bas the spatial-temporal profile realization as in Figure 5. The 
two slabs (unstable slices of Si) Si and Si^a projected onto the plane of "cosx", are 
also illustrated in Figure 2. From these figures, we can see clearly that the chaotic 
center-wing jumping (Figure 1) is the realization of the shift automorphism on four 
symbols 1, 2, -1, -2 in Si U Si^^- 

7.6. Shadowing Lemma and Chaos in Finite-D Periodic Systems 

Since its invention [9], shadowing lemma has been a useful tool for solving 
many dynamical system problems. Here we only focus upon its use in proving the 
existence of chaos in a neighborhood of a transversal homoclinic orbit. According 
to the type of the system, the level of diSiculty in proving the existence of chaos 
with a shadowing lemma is different. 

A finite-dimensional periodic system can be written in the general form 

X = F{x,t) , 

where x G M", and F{x,t) is periodic in t- Let / be the Poincare period map. 
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Figure 2 . Figure-eight structure and its correspondence with the 
chaotic center-wing jumping. 
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Figure 3. Spatial- temporal profile realization of Lc in Figure 2 
(coordinates are the same as in Figure 1). 




Figure 4. Spatial-temporal profile realization of Lcin in Figure 
2 (coordinates are the same as in Figure 1). 

Definition 7.8. A doubly infinite sequence {y^} in R" is a 5 pseudo-orbit of 
a diffeomorphism / : R" i— > M" if for all integers j 

An orbit {/■'(a;)} is said to e-shadow the 5 pseudo-orbit {yj} if for all integers j 

\r{x)-y^\<e. 

Definition 7.9. A compact invariant set S is hyperbolic if there are positive 
constants K, a and a projection matrix valued function P{x), x € 5, of constant 
rank such that for all a; in 5 

P{f{x))Df{x) = Df{x)P{x) , 
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Figure 5. Spatial- temporal profile realization of Lout in Figure 2 
(coordinates are the same as in Figure 1). 

\Dp{x)P{x)\<Ke-"^ , (j>0), 
\DP{x)iI - P{x))\ < Ke"^ , (j < 0) . 

Theorem 7.10 (Shadowing Lemma [173]). Let S be a compact hyperbolic set 
for the diffeomorphism / : R" i— > M". Then given e > sufficiently small there 
exists S > such that every S pseudo-orbit in S has a unique e-shadowing orbit. 

The proof of this theorem by Palmer [173] is overall a fixed point argument 
with the help of Green functions for linear maps. 

Let yo be a transversal homoclinic point asymptotic to a saddle xo of a 
diffeomorphism / : K" i-^ K". Then the set 

S = {xo} U {f{yo) :j€Z} 

is hyperbolic. Denote by Aq and Ai the two orbit segments of length 2m + 1 

Ao = {xo,xo,--- ,xo} , Ai = {f-^{yo),f-"'+'{yo),--- ,r-Hyo),riyo)} ■ 
Let 

a = {■■■ ,a_i,ao,ai,- • • ) , 

where aj G {0,1}, be any doubly infinite binary sequence. Let A be the doubly 
infinite sequence of points in S, associated with a 

A = {• • • , , A^Q , Aai , • • • } • 

When m is sufficiently large, A is a S pseiido-orbit in S. By the shadowing lemma 
(Theorem 7.10), there is a unique e-shadowing orbit that shadows A. In this man- 
ner, Palmer [173] gave a beautiful proof of Smale's horseshoe theorem. 

Definition 7.11. Denote by S the set of doubly infinite binary sequences 

a= {■■■ ,a_i,ao,ai, ■■■) , 
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where aj G {0, 1}. We give the set {0, 1} the discrete topology and S the product 
topology. The Bernoulli shift x is defined by 

[x{a)]'j = aj+i . 

Theorem 7.12. Let yo be a transversal homoclinic point asymptotic to a saddle 
Xo of a diffeomorphism f : R" i-^ R" . Then there is a homeomorphism cj) ofE 
onto a com,pact subset of R" which is invariant under f and such that when m is 
sufficiently large 

that is, the action of f^™~^^ on is topologically conjugate to the action of x on 
E. 

Here one can define (l){a) to be the point on the shadowing orbit that shadows 
the midpoint of the orbit segment Aag, which is either xo or yo- The topological 
conjugacy can be easily verified. For details, see [173]. Other references can be 
found in [172] [213]. 

7.7. Shadowing Lemma and Chaos in Infinite-D Periodic Systems 

An infinite-dimensional periodic system defined in a Banach space X can be 
written in the general form 

X = F{x,t) , 

where x G X, and F{x,t) is periodic in t. Let / be the Poincare period map. 
When / is a map which needs not to be invertible, shadowing lemma and 
symbolic dynamics around a transversal homoclinic orbit can both be established 
[193] [194] [83]. Other references can bo found in [79] [43] [214] [30]. There exists 
also a work on horseshoe construction without shadowing lemma for sinusoidally 
forced vibrations of buckled beam [86] . 

7.8. Periodically Perturbed Sine-Gordon (SG) Equation 

We continue from section 6.9, and use the notations in section 7.6. For the 
periodically perturbed sine-Gordon equation (6.51), the Poincare period map is a 

diffeomorphism in H^. As a corollary of the result in last section, we have the 
theorem on the existence of chaos. 

Theorem 7.13. There is an integer m and a homeomorphism <p of T, onto a 
compact Cantor subset A of . A is invariant under the Poincare period-2'K map 
P of the periodically perturbed sine-Gordon equation (6.51). The action of P^"^^^ 
on A is topologically conjugate to the action of x on T, : p2m+i ^ ^ _ ^ ^ That 
is, the following diagram commutes: 
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7.9. Shadowing Lemma and Chaos in Finite-D Autonomous Systems 

A finite-dimensional autonomous system can be written in the general form 

X = F{x) , 

where x e M". In this case, a transversal homoclinic orbit can be an orbit asymp- 
totic to a normally hyperbolic limit cycle. That is, it is an orbit in the intersection 
of the stable and unstable manifolds of a normally hyperbolic limit cycle. Instead 
of the Poincare period map as for periodic system, one may want to introduce the 
so-called Poincare return map which is a map induced by the flow on a codimcn- 
sion 1 section which is transversal to the limit cycle. Unfortunately, such a map 
is not even well-defined in the neighborhood of the homoclinic orbit. This poses a 
challenging difficulty in extending the arguments as in the case of a Poincare period 
map. In 1996, Palmer [174] completed a proof of a shadowing lemma and existence 
of chaos using Newton's method. It will be difficult to extend this method to infi- 
nite dimensions, since it used heavily differentiations in time. Other references can 
be found in [50] [51] [186] [71]. 

7.10. Shadowing Lemma and Chaos in Infinite-D Autonomous Systems 

An infinite-dimensional autonomous system defined in a Banach space X can 
be written in the general form 

X = F{x) , 

where x G X. In 2002, the author [125] completed a proof of a shadowing lemma 
and existence of chaos using Fenichel's persistence of normally hyperbolic invariant 
manifold idea. The setup is as follows, 

• Assumption (Al): There exist a hyperbolic limit cycle S and a transver- 
sal homoclinic orbit ^ asymptotic to S. As curves, S and ^ are C^. 

• Assumption (A2): The Fenichel fiber theorem is valid at S. That is, 
there exist a family of unstable Fenichel fibers : q E S} and a 
family of stable Fenichel fibers {J^^{q) : q G S}. For each fixed q & S, 
J='"{q) and J^^iq) are submanifolds. and are in q,\fq e 
S. The unions Uges-^"('Z) Uqes-^^C?) ^re the unstable and stable 
manifolds of S. Both families arc invariant, i.e. 

F\j^-{q))(lT-{FHq)),yt<0,qeS, 

F\j^'{q))CJ'%F\q)),yt>0,qGS, 

where is the evolution operator. There are positive constants k and C 
such that Vq G 5, Vg" G and Vg+ G T'iq), 

\\F\q~)-F\q)\\<Ce-'\\q- -q\\yt<0 , 

\\F'iq+) - F'iq)\\ < de--'\\q+ - q\\,y t >0 . 

• Assumption (A3): F*(q) is C° in t, for t G (—00,00), q e X. For any 
fixed t G (—00,00), F*{q) is a diffeomorphism on X. 

Remark 7.14. Notice that we do not assume that as functions of time, S and 
^ are , and we only assume that as curves, S and ^ are C^. 

Under the above setup, a shadowing lemma and existence of chaos can be 
proved [125]. Another crucial element in the argument is the establishment of a 
A-lemma (also called inclination lemma) which will be discussed in a later section. 
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7.11. A Derivative Nonlinear Schrodinger Equation 

We continue from section 6.10, and consider tiie derivative nonlinear Schrodinger 
equation (6.52). The transversal homoclinic orbit given in Theorem 6.8 is a classical 
solution. Thus, Assumption (Al) is valid. Assumption (A2) follows from the stan- 
dard arguments in [139] [136] [124]. Since the perturbation in (6.52) is bounded. 
Assumption (A3) follows from standard arguments. Thus there exists chaos in the 
derivative nonlinear Schrodinger equation (6.52) [125]. 

7.12. A-Lemma 

A-Lcmma (inclination lemma) has been utilized in proving many significant 
theorems in dynamical system [171] [170] [197]. Here is another example, in [125], 
it is shown that A-Lemma is crucial for proving a shadowing lemma in an infinite 
dimensional autonomous system. Below we give a brief introduction to A-Lemma 
[171]. 

Let / be a C" (r > 1) diffeomorphism in with as a saddle. Let and 

i?" be the stable and unstable subspaces. Through changing coordinates, one can 
obtain that for the local stable and unstable manifolds, 

W^loc^^^ ^locCi^"- 
Let c Wf^^ and S" c W^^^ be balls, and V = x B". Let q e Wf^J{0}, 
and be a disc of the same dimension as B", which is transversal to W^^^ at q. 

Lemma 7.15 (The A-Lemma). Let be the connected component of f"'{D")n 
V to which f^{q) belongs. Given e > 0, there exists no such that if n > no, then 

Dl is e -close to S". 

The proof of the lemma is not complicated, yet nontrivial [171]. In different 
situations, the claims of the A-lemma needed may be different [171] [125]. Never- 
theless, the above simple version illustrated the spirit of A-lemmas. 

7.13. Homoclinic Tubes and Chaos Cascades 

When studying high-dimensional systems, instead of homoclinic orbits one is 
more interested in the so-called homoclinic tubes [126] [127] [119] [128] [130]. The 
concept of a homoclinic tube was introduced by Silnikov [187] in a study on the 
structure of the neighborhood of a homoclinic tube asymptotic to an invariant torus 
a under a diffeomorphism F in a finite dimensional phase space. The asymptotic 
torus is of saddle type. The homoclinic tube consists of a doubly infinite sequence 
of tori {(Tj. i = 0, ±1,±2,- • •} in the transversal intersection of the stable and 
unstable manifolds of a, such that cTj+i = F oaj for any j. It is a generalization of 
the concept of a transversal homoclinic orbit when the points are replaced by tori. 

We are interested in homoclinic tubes for several reasons [126] [127] [119] 
[128] [130]: 1. Especially in high dimensions, dynamics inside each invariant 
tubes in the neighborhoods of homoclinic tubes are often chaotic too. We call 
such chaotic dynamics chaos in the sm,aW , and the symbolic dynamics of the in- 
variant tubes chaos in the large". Such cascade structures are more important 
than the structures in a neighborhood of a homoclinic orbit, when high or infi- 
nite dimensional dynamical systems are studied. 2. Symbolic dynamics structures 
in the neighborhoods of homoclinic tubes are more observable than in the neigh- 
borhoods of homoclinic orbits in numerical and physical experiments. 3. When 
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studying high or infinite dimensional Hamiltonian system (for example, the cubic 
nonlinear Schrodinger equation under Hamiltonian perturbations), each invariant 
tube contains both KAM tori and stochastic layers (chaos in the small). Thus, not 
only dynamics inside each stochastic layer is chaotic, all these stochastic layers also 
move chaotically under Poincare maps. 

Using the shadowing lemma technique developed in [125], we obtained in [126] 
a theorem on the symbolic dynamics of submanifolds in a neighborhood of a ho- 
moclinic tube under a C^-diffeomorphism defined in a Banach space. Such a proof 
removed an uncheckable assumption in [187]. The result of [187] gives a sym- 
bolic labeling of all the invariant tubes around a homoclinic tube. Such a symbolic 
labeling docs not imply the symbolic dynamics of a single map proved in [126]. 

Then in [127], as an example, the following sine-Gordon equation under chaotic 
perturbation is studied, 

(7.10) utt = c^Uxx + smu + e[-au + f{t,62,d3,0i){smu - u)] , 
which is subject to periodic boundary condition and odd constraint 

(7.11) u{t,x + 2n) =u{t,x) , u{t,—x) = —u{t,x), 

where u is a real-valued function of two real variables t > and x, c is a parameter, 
i<c<l, a>Oisa parameter, e > is a small parameter, / is periodic in 
t, 02, 03, 04, and (6*2, 6*3, 04) e T*^, by introducing the extra variable 0i = LOit + 0°, f 
takes the form 

4 

/(*) = '^anCOs[0n{t)] , 

n=l 

and Un's are parameters. Let 0n — i^nt + + e^^^nj n = 2, 3, 4, > 1, and i?„'s are 
given by the ABC flow [57] which is verified numerically to be chaotic, 

'&2 = j4sini?4 -I- Ccosi?3 , 

■!?3 = B sin ^2 + A cos 1^4 , 

^4 = C sin + _B cos 1^2 , 

where A, B, and C are real parameters. Existence of a homoclinic tube is proved. 
As a corollary of the theorem proved in [126], Symbolic dynamics of tori around 
the homoclinic tube is established, which is the "chaos in the large". The chaotic 
dynamics of {02, 03, 6*4) is the "chaos in the small". 

Here we see the embedding of smaller scale chaos in larger scale chaos. By 
introducing more variables, such embedding can be continued with even smaller 
scale chaos. This leads to a chain of embeddings. We call this chain of embeddings 
of smaller scale chaos in larger scale chaos, a "chaos cascade" . We hope that such 
"chaos cascade" will be proved important. 
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Unlike soliton equations under periodic boundary conditions, phase space struc- 
tures, especially hyperbolic structures, of soliton equations under decay boundary 
conditions are not well understood. The application of Backlund-Darboux trans- 
formations for generating hyperbolic foliations is not known in decay boundary 
condition case. Nevertheless, we believe that Backlund-Darboux transformations 
may still have great potentials for understanding phase space structures. This can 
be a great area for interested readers to work in. The common use of Backlund- 
Darboux transformations is to generate multi-soliton solutions from single soliton 
solutions. The relations between different soliton solutions in phase spaces (for 
example certain Sobolev spaces) are not clear yet. There are some studies on the 
linear and nonlinear stabilities of traveling- wave (soliton) solutions [175] [148]. Re- 
sults concerning soliton equations arc only stability results. Instability results are 
obtained for non-soliton equations. On the other hand, non-soliton equations are 
not integrable, have no Lax pair structures, and their phase space structures are 
much more difficult to be understood. The challenging problem here is to identify 
which soliton equations possess traveling-wave solutions which are unstable. The 
reason why we emphasize instabilities is that they are the sources of chaos. Below 
we are going to review two types of studies on the phase space structures of soliton 
equations under decay boundary conditions. 



If we only consider traveling-wave solutions to the soliton equations or gen- 
eralized ( or perturbed ) soliton equations, the resulting equations are ordinary 
differential equations with parameters. These ODEs often offer physically mean- 
ingful and mathematically pleasant problems to be studied. One of the interesting 
features is that the traveling- wave solutions correspond to homoclinic orbits to such 
ODEs. We also naturally expect that such ODEs can have chaos. Consider the 
KdV equation 

(8.1) Ut + QuUx + Uxxx = ^ , 

and let u{x, t) = U{^), ^ = x — ct {c is a real parameter), we have 



where Ci is a real integration constant. This equation can be rewritten as a system, 



8.1. Traveling Wave Reduction 



U'" + 6UU' -cU' = . 




U" -cU + ci=0 , 



(8.3) 
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For example, when ci = and c > we have the sohton, 

(8.4) = I sech'[±y^C] , 

which is a homochnic orbit asymptotic to 0. 

We may have an equation of the following form describing a real physical phe- 
nomenon, 

Ut + [f{u)]x + Uxxx + egiu) = , 

where e is a small parameter, and traveling waves are physically important. Then 

we have 

U"' + [f(U)]'-cU' + eg{U)^Q, 
which can be rewritten as a system, 

U' = V , 

v = w , 

W' = -fV + cV- eg{U) . 

This system may even have chaos. Such low dimensional systems are very popular 
in current dynamical system studies and are mathematically very pleasant [94] [95]. 
Nevertheless, they represent a very narrow class of solutions to the original PDE, 
and in no way they can describe the entire phase space structures of the original 
PDE. 

8.2. Stabilities of the Traveling- Wave Solutions 

The first step toward understanding the phase space structures of soliton equa- 
tions under decay boundary conditions is to study the linear or nonlinear stabilities 
of traveling- wave solutions [175] [148] [75] [98]. Consider the KdV equation (8.1), 
and study the linear stability of the traveling wave (8.4). First we change the vari- 
ables from (x, t) into (^, t) where ^ — x — ct; then in terms of the new variables 
(^,i), the traveling wave (8.4) is a fixed point of the KdV equation (8.1). Let 
u{x,t) = Us{£,) +f(^,i) and linearize the KdV equation (8.1) around C/s(^), we 
have 

dtv - cd^v + 6[Usd^v + vd^U^] + dfv = . 

We seek solutions to this equation in the form v = e'^*Y{^) and we have the eigen- 
value problem; 

(8.5) XY - cd^Y + 6[Usd^Y + Yd^U^] + d^Y = . 

One method for studying such eigenvalue problem is carried in two steps [175]: 
First, one studies the |^] — > oo limit system 

XY - cd^Y + dlY = , 

and finds that it has solutions Y{^) = e'*^^ for j = 1, 2, 3 where the /ij satisfy 

Re{/ii} < < Re{w} for 1 = 2,3. 

Thus the equation (8.5) has a one-dimensional subspace of solutions which decay 

as X — > oo, and a two-dimensional subspace of solutions which decay as a; ^ — oo. 
A will be an eigenvalue when these subspaces meet nontrivially. In step 2, one 
measures the angle between these subspaces by a Wronskian-like analytic function 
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D{X), named Evans's function. One interpretation of this function is that it is hke 
a transmission coefficient, in the sense that for the solution of (8.5) satisfying 

y(C) ~ e''^« as C ^ oo , 

we have 

F(0 ~ i?(A)e^i^ as C ^ -cx) . 

In equation (8.5), for Re{A} > 0, if D{X) vanishes, then A is an eigenvalue, and 
conversely. The conclusion for the equation (8.5) is that there is no eigenvalue with 
Rc{A} > 0. In fact, for the KdV equation (8.1), Us is F^-orbitally stable [26] [116] 
[199] [117]. 

For generalized soliton equations, linear instability results have been established 
[175]. Let f{u) = uP+'^/{p + 1), and 

f/s(0 = asech2/P(60 

for appropriate constants a and b. For the generalized KdV equation [175] 

dtu + dj{u) + dlu = , 

if p > 4, then Us is linearly unstable for all c > 0. If 1 < p < 4, then Us is 
iJ^-orbitally stable. For the generalized Benjamin-Bona-Mahoney equation [175] 

dtu + dxu + dxf{u) - dtdlu = , 

if p > 4, there exists a positive number Co(p), such that Us with 1 < c < co(p) are 
linearly unstable. For the generalized regularized Boussinesq equation [175] 

d^,u-dlu-dlf{u)-d^,dlu = Q, 

if p > 4, Us is linearly unstable if 1 < < Co(p), where 

cl{p) = ^p/{A + 2p) . 

8.3. Breathers 

Breathers are periodic solutions of soliton equations under decay boundary 
conditions. The importance of breathers with respect to the phase space structures 
of soliton equations is not clear yet. Backlund-Darboux transformations may play a 
role for this. For example, starting from a breather with one frequency, Backlund- 
Darboux transformations can generate a new breather with two frequencies, i.e. 
the now breather is a quasipcriodic solution. The stability of breathers will be a 
very important subject to study. 

First we study the reciprocal relation between homoclinic solutions generated 
through Backlund-Darboux transformations for soliton equations under periodic 
boundary conditions and breathers for soliton equations under decay boundary 
conditions. Consider the sine-Gordon equation 

(8.6) utt - Uxx + sinu = , 
under the periodic boundary condition 

(8.7) u{x + L,t) = u{x,t) , 
and the decay boundary condition 

(8.8) u{x,t) — > as |a;| — > oo . 
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Starting from the trivial solution u = ir, one can generate the homoclinic solution 
to the Cauchy problem (8.6) and (8.7) through Backlund-Darboux transformations 
[59], 

tanj/cos[(cos i')x] 



(8.9) 



UH{x,t) = IT + 4tan 



cosh[(sin i^)t] 

where i = 27r/ cosu. Applying the sine-Gordon symmetry 

(a;, t, u) — > {t, x,u — n) , 

one can immediately generate the breather solution to the Cauchy problem (8.6) 
and (8.8) [59], 

tan:^cos[(cosi^)t] 



(8.10) 



UB{x,t) = 4 tan 



cosh[(sin iy)x] 

The period of this breather is L = 27r/ cos v. The linear stability can be investigated 
through studying the Floquct theory for the linear partial differential equation 
which is the linearized equation of (8.6) at ub (8.10): 

(8.11) utt - Uxx + (cosub)u = , 
under the decay boundary condition 

(8.12) u{x,t) — > as |a;| — > oo . 

Such studies will be very important in terms of understanding the phase space 

structures in the neighborhood of the breather and developing infinite dimensional 
Floquet theory. In fact, McLaughlin and Scott [156] had studied the equation 
(8.11) through infinitessimal Backlund transformations, which lead to the so-called 
radiation of breather phenomenon. In physical variables only, the Backlund trans- 
formation for the sine-Gordon equation (8.6) is given as, 



d_ 

d_ 



2 

u+ — U- 



— - sm 



= iC, sin 



u+ 



2 

u+ + U- 



where ^ = ^{x + t), r] = ^{x — t), and is a complex Backlund parameter. Through 



the nonlinear superposition principle [179] [8], 

C1+C2 



.13) 



tan 



tan 



Ui - U2 



C1-C2 

built upon the Bianchi diagram [179] [8]. The breather ub (8.10) is produced by 
starting with U- = 0, and choosing ^1 = cosf + isinv and C2 = — cosv + isinu. 
Let U- be a solution to the linearized sine-Gordon equation (8.11) with ub replaced 
by u- = 0, 

u- = e'(''t+^^^^) . 

Then variation of the nonlinear superposition principle leads to the solution to the 
linearized sine-Gordon equation (8.11): 

-1 



u+ 



U- 



{1 + Ay 



-I- ^cos2(mb/4) 
sin^ Pt 



ai + a2 



cos [(Wi - W2)/4] 



^ as (cosh^ ax)' ^ 



cosh ax 2 cosh ax 



(ui - U2) 

ai (sin^ pty 
2 cosh ax 
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where 



a = sin 1/ , 13 = cos v 



A 



a sin j3t 



P cosh ax 




|2 ' 



which represents the radiation of the breather. 

Next we briefly survey some interesting studies on breathers. The rigidity 
of sine-Gordon breathers was studied by Birnir, Mckean and Weinstein [29] [28]. 

Under the scahng 

X ^ x' = {l + eay^/^x , t ^t' = + , u ^ u' = {1 + €b)u 

the sine-Gordon equation (8.6) is transformed into 



The breather (8.10) is transformed into a breather for the new equation. The 
remarkable fact is the rigidity of the breather (8.10) as stated in the following 
theorem. 

Theorem 8.1 (Birnir, Mckean and Weinstein [29]). If f{u) vanishes at u = 
and is holomorphic in the open strip |i?e{w}| < tt, then uu — u^x + sinu + ef{u) + 
O(e^) = has breathers u = ub + eui + 0{e^), one to each s\nv from an open 
subinterval of (0, l/v^] issuing smoothly from u = ub at e = 0, with ui vanishing 
at X = ±oo, if and only if f{u) is a linear combination of sin u, ucosu, and/or 
1 + ScosM — 4cos(u/2) + 4cos'ulncos(u/4). 

In fact, it is proved that sine-Gordon equation (8.6) is the only nonlinear wave 
equation possessing small analytic breathers [104]. These results lend color to the 
conjecture that, for the general nonlinear wave equation, breathing is an "arithmeti- 
cal" phenomenon in the sense that it takes place only for isolated nonlinear ities, 
scaling and other trivialities aside. 



Utt — Uxx + sinu -|- e(asinM -|- 6ucosu) -|- O(e^) = . 



CHAPTER 9 



Lax Pairs of Euler Equations of Inviscid Fluids 

The governing equations for the incompressible viscous fluid flow are the Navier- 
Stokes equations. Turbulence occurs in the regime of high Reynolds number. By 
formally setting the Reynolds number equal to infinity, the Navier-Stokes equations 
reduce to the Euler equations of incompressible inviscid fluid flow. One may view 
the Navier-Stokes equations with large Reynolds number as a singular perturbation 
of the Euler equations. 

Results of T. Kato show that 2D Navier-Stokes equations are globally well- 
posed in C"([0,oo);i?''(i?2)), s > 2, and for any < T < oo, the mild solu- 
tions of the 2D Navier-Stokes equations approach those of the 2D Euler equations 
in C°{[0,T];H%R'^)) [101]. 3D Navier-Stokes equations are locally well-posed in 
C°{[0,T];H''iR^)), s > 5/2, and the mild solutions of the 3D Navier-Stokes equa- 
tions approach those of the 3D Euler equations in C'^ {[0, t]; (R^)) , where r de- 
pends on the norms of the initial data and the external force [99] [100]. Extensive 
studies on the inviscid limit have been carried by J. Wu et al. [204] [49] [205] [31]. 
There is no doubt that mathematical study on Navier-Stokes (Euler) equations is 
one of the most important mathematical problems. In fact, Clay Mathematics In- 
stitute has posted the global well-posedness of 3D Navier-Stokes equations as one 
of the one million dollars problems. 

V. Arnold [15] realized that 2D Euler equations are a Hamiltonian system. 
Extensive studies on the symplectic structures of 2D Euler equations have been 
carried by J. Marsdcn, T. Ratiu et al. [151]. 

Recently, the author found Lax pair structures for Euler equations [123] [141] 
[131] [128]. Understanding the structures of solutions to Euler equations is of 
fundamental interest. Of particular interest is the question on the global well- 
posedness of 3D Navier-Stokes and Euler equations. Our number one hope is that 
the Lax pair structures can be useful in investigating the global well-posedness. 
Our secondary hope is that the Darboux transformation [141] associated with the 
Lax pair can generate explicit representation of homoclinic structures [121]. 

The philosophical significance of the existence of Lax pairs for Euler equations is 
even more important. If one defines intcgrability of an equation by the existence of 
a Lax pair, then both 2D and 3D Euler equations arc integrable. More importantly, 
both 2D and 3D Navier-Stokes equations at high Reynolds numbers are singularly 
perturbed integrable systems. Such a point of view changes our old ideology on 
Euler and Navier-Stokes equations. 

9.1. A Lax Pair for 2D Euler Equation 

The 2D Euler equation can be written in the vorticity form, 
(9.1) dtfl + {^,n} = , 
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where the bracket { , } is defined as 

{f,9} = id.f){dyg)-{dyf)id,g) , 
fl is the vorticity, and is the stream function given by, 



u = -dy^ , V = 



and the relation between vorticity fl and stream function ^ is, 

Sl = dxV- dyU = A* . 
Theorem 9.1 (Li, [123]). The Lax pair of the 2D Euler equation (9.1) is given 



Lip = \if , 
dt^fi + Aip = Q 



(9.2) 
where 

Lip^{n,<p}, Aip = {^,if,}, 

and X is a complex constant, and (p is a complex-valued function. 

9.2. A Darboux Transformation for 2D Euler Equation 

Consider the Lax pair (9.2) at A = 0, i.e. 

(9.3) {n,p} = o, 

(9.4) dtp + {^,p} = 0, 
where we replaced the notation ip by p. 

Theorem 9.2 ([141]). Let f = f{t,x,y) be any fixed solution to the system 
(9.3, 9.4), we define the Gauge transform Gf: 

(9.5) p = G/p = ^ [p, - {d. In f)p] , 

and the transforms of the potentials SI and ^: 

(9.6) * = * + F, n = n + AF, 
where F is subject to the constraints 

(9.7) {n, AF} = , {AF, F} = . 

Then p solves the system (9.3, 9.4) at (fi,^'). Thus (9.5) and (9.6) form the 
Darboux transformation for the 2D Euler equation (9.1) and its Lax pair (9.3, 9.4). 

Remark 9.3. For KdV equation and many other soliton equations, the Gauge 
transform is of the form [152], 

P = Px- {da;lnf)p . 

In general. Gauge transform does not involve potentials. For 2D Euler equation, a 
potential factor is needed. Prom (9.3), one has 

^ Pv_ 

Six Sly 

The Gauge transform (9.5) can be rewritten as 
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The Lax pair (9.3, 9.4) has a symmetry, i.e. it is invariant under the transform 
{t,x,y) {—t,y,x). The form of the Gauge transform (9.5) resulted from the 
inclusion of the potential factor jp, is consistent with this symmetry. 

Our hope is to use the Darboux transformation to generate homoclinic struc- 
tures for 2D Euler equation [121]. 

9.3. A Lax Pair for Rossby Wave Equation 

The Rossby wave equation is 

dtn + {^,fl} + (3d^'i' = , 

where Cl = n{t,x,y) is the vorticity, = ^x^y - '^y^x, and * = A~^n 

is the stream function. Its Lax pair can be obtained by formally conducting the 
transformation, f2 = i7 + f3y, to the 2D Euler equation [123], 

where ip is a complex-valued function, and A is a complex parameter. 

9.4. Lax Pairs for 3D Euler Equation 
The 3D Euler equation can be written in vorticity form, 

(9.8) dtfi + {u ■ -{n-V)u = , 

where u = (^1,^2,^3) is the velocity, O = (fii, 1^2, ila) is the vorticity, V = 
{dx,dy,dz), = V X u, and V • u = 0. u can be represented by O for exam- 
ple through Biot-Savart law. 

Theorem 9.4. The Lax pair of the 3D Euler equation (9.8) is given as 
(9 9) r = , 

where 

Lcl)= (n- V)(j) , A(p = (u ■ V)^ , 
X is a complex constant, and (p is a complex scalar-valued Junction. 

Theorem 9.5 ([39]). Another Lax pair of the 3D Euler equation (9.8) is given 

as 

(9-10) I 5^1^/' n 

where 

Lip={Q- V)ip - {ip ■ V)0 , Aip={u- V)ip - {if ■ V)u , 
X is a complex constant, and = (<pi, (^2, ^3) is a complex 3-vector valued function. 

Our hope is that the infinitely many conservation laws generated by A G C can 
provide a priori estimates for the global well-posedness of 3D Navier-Stokes equa- 
tions, or better understanding on the global well-posedness. For more informations 
on the topics, see [141]. 



CHAPTER 10 



Linearized 2D Euler Equation at a Fixed Point 



To begin an infinite dimensional dynamical system study, we investigate the 
linearized 2D Euler equation at a fixed point [122]. We consider the 2D Euler 
equation (9.1) under periodic boundary condition in both x and y directions with 
period 2-k. Expanding ^ into Fourier series, 



feez2/{o} 



ik-X 



where ui-k = Uk , k = {ki, ^2) , and X — {x,y) . The 2D Euler equation can be 
rewritten as 



(10.1) 

where A{p, q) is given by, 
(10.2) 



k=p+q 



k\ ^ - \p\ ^](pi92 -P2gi) 



\p\- 



Pi 91 
P2 q2 



where = Qi+Q^ fo^' Q — (Sii 92)^, similarly for p. The 2D Euler equation (10.1) 
has huge dimensional equilibrium manifolds. 

Proposition 10.1. For any k e Z^/{0}, the infinite dimensional space 



Wk'} 

and the finite dimensional space 



ujk' =0, ifk'j^rk, Vr e i? 



El = !^{ujk'} ojk'=0, if\k'\j^\k\^ , 

entirely consist of fixed points of the system (10.1). 

Fig.l shows an example on the locations of the modes (fc' = rk) and (|fc'| = |fc|) 
in the definitions of E^^ and E'^ (Proposition 10.1). 

10.1. Hamiltonian Structure of 2D Euler Equation 

For any two functional F\ and F2 of {uJk}^ define their Lie-Poisson bracket: 



(10.3) 



{F,,F,}= 

fe+P+g=0 



?l P\ 
q2 P2 



<^k 



dFi dF2 
dUJ^ dUJ^ 
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Then the 2D Euler equation (10.1) is a Hamiltonian system [15], 

(10.4) Uk^ioJk.H}, 
where the Hamiltonian H is the kinetic energy, 

(10.5) H=\ \k\-^W?. 

feezvio} 

Following are Casimirs (i.e. invariants that Poisson commute with any functional) 
of the Hamiltonian system (10.4): 

(10.6) J„ = ^ UJki---UJkn- 

fciH hfc„=0 

10.2. Linearized 2D Euler Equation at a Unimodal Fixed Point 

Denote {oJk}kez^/{o} by We consider the simple fixed point u* [122]: 

(10.7) w;=r, u;*=0, iffc^por ~ p, 

of the 2D Euler equation (10.1), where F is an arbitrary complex constant. The 
linearized two-dimensional Euler equation at oj* is given by, 

(10.8) Wk = Mp^ k-p)T uik-p + A{-p, k+p)f Wk+p ■ 

Definition 10.2 (Classes). For any k e Z^/{0}, we define the class to be 
the subset of ZV{0}: 

= |fc + np e ZV{0} neZ, p is specified in (10.7) I . 
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See Fig. 2 for an illustration of the classes. According to the classification defined 
in Definition 10.2, the linearized two-dimensional Euler equation (10.8) decouples 
into infinitely many invariant subsystems: 

(10-9) + A{-p, k + {n+ l)p) f • 




10.2.1. Linear Hamiltonian Systems. Each invariant subsystem (10.9) can 
be rewritten as a linear Hamiltonian system as shown below. 

Definition 10.3 (The Quadratic Hamiltonian). The quadratic Hamiltonian 
is defined as: 



(10.10) 







Pi 


k 


P2 


k 


2 _ 





'^k+(n-l)p '^k+np 



where Pn = [\k + np\ — \p\ ], " Im " denotes " imaginary part ". 
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Then the invariant subsystem (10.9) can be rewritten as a hnear Hamiltonian 
system [122], 

(10-11) ^^k+nr, = P^' 



k+np 



For a finite dimensional Hnear Hamiltonian system, it is well-known that the eigen- 
values are of four types: real pairs (c, — c), purely imaginary pairs {id, —id), quadru- 
ples {±c±id), and zero eigenvalues [177] [142] [17]. There is also a complete 
theorem on the normal forms of such Hamiltonians [17]. For the above infinite 
dimensional system (10.9), the classical proofs can not be applied. Nevertheless, 
the conclusion is still true with different proof [122] [134]. Let £^ be the linear 
operator defined by the right hand side of (10.9), and be the Sobolev space 
where s > is an integer and H'^ ~ £2- 

Theorem 10.4 ([122] [134]). The eigenvalues of the linear operator Cj, in 
are of four types: real pairs (c,—c), purely imaginary pairs (id,— id), quadruples 
C±c± id), and zero eigenvalues. 

10.2.2. Liapunov Stability. 

Definition 10.5 (An Important Functional). For each invariant subsystem 
(10.9), we define the functional /j, 



h = 



■'^(restricted to Sj) 



(10.12) 



np\ 



\P\-'} 



Lemma 10.6. If. is a constant of motion for the system (10.11). 
Definition 10.7 (The Disk). The disk of radius \p\ in 1? /{Q}, denoted by 



£>|j,|, is defined as 



D 



k G ZV{0} 



See Fig. 2 for an illustration. 



Theorem 10.8 (Unstable Disk Theorem, [122]). // n D|p| = 0, then the 
invariant subsystem (10.9) is Liapunov stable for all t G R, in fact, 



where 



nez 



a = 



10 



k+np 



it) 



< 



CO 



k-\-np 



(0) 



nez 



max{-p„} 

nez 



min{-p„} 

nez 



yt€R, 



< (7 < 00. 



10.2.3. Spectral Theorems. Again denote by £^ the linear operator defined 

by the right hand side of (10.9). 

Theorem 10.9 (The Spectral Theorem, [122] [134]). We have the following 
claims on the spectrum of the linear operator Cf, : 
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(1) IfSf. n -D|p| = 0, then the entire spectrum of the linear operator Cf, is 



its continuous spectrum. See Figure 3, where b = — ^|r| 



Pi 
Pi 



ki 
k 



That is, both the residual and the point spectra of £^ are empty. 
(2) // n -D|p| ^ 0, then the entire essential spectrum of the linear 
operator Cj^ is its continuous spectrum. That is, the residual spectrum of 
Lf. is empty. The point spectrum of Lf, is symmetric with respect to both 
real and imaginary axes. See Figure 4- 



9{A} 





i2\b\ 




* 




-i2\b\ 



3?{A} 



Figure 3. The spectrum oi C; in case (1). 



For a detailed proof of this theorem, sec [122] [134] . Denote by L the right hand 
side of (10.8), i.e. the whole linearized 2D Euler operator, the spectral mapping 
theorem holds. 



Theorem 10.10 ([112]). 



a(e*^) = e*'^(^),i7^0. 



10.2.4. A Continued Fraction Calculation of Eigenvalues. Since the 

introduction of continued fractions for calculating the eigenvalues of steady fluid 
flow, by Meshalkin and Sinai [161], this topics had been extensively explored [208] 
[143] [144] [145] [146] [147] [22] [122]. Rigorous justification on the continued 
fraction calculation was given in [122] [147]. 
Rewrite the equation (10.9) as follows. 



(10.13) 



Pn 



'71+ 1 
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3{A} 



• 

• • 


i2|6| 

• 

• • 




* 






• • 


• • 


• 






-i2\b\ 



U{X} 



Figure 4. The spectrum of £^ in case (2). 



where Zn = P„e'"(«+'^/2)a;^^„p, ^ + 7 = n/2, T = \T\e'^, a = ^\T\ 
Pn = \k + np|~^ — Let Zn = e^^Zn, where X G C; then Zn satisfies 

(10.14) UnZn + Zn-1 - Zn+1 = , 

where a„ = X{apn)~^- Let w„ = Zn/zn-i [161]; then w„ satisfies 



Pi ki 
P2 h 



an-\ = Wn+l 

Wr,. 



(10.15) 

Iteration of (10.15) leads to the continued fraction solution [161], 
(10.16) 



an-2 + 



an-3 + 



Rewrite (10.15) as follows, 
(10.17) 



Wr, 



-an + Wn+l 

Iteration of (10.17) leads to the continued fraction solution [161], 



10.2. LINEARIZED 2D EULER EQUATION AT A UNIMODAL FIXED POINT 



101 



The eigenvalues are given by the condition u'[^'' = i.e. 
(10.19) / = ao + { ^ - ) + { \ ) = , 



where / = /(A, fc,p), A = A/a. 

As an example, wo take p = (1, 1)-^. When F 0, the fixed point has 4 
eigenvalues which form a quadruple. These four eigenvalues appear in the invariant 
linear subsystem labeled by = (—3, —2)^. One of them is [122]: 

(10.20) A = 2A/|r| = 0.24822302478255 + i 0.35172076526520 . 

See Figure 5 for an illustration. The essential spectrum (= continuous spectrum) 
of Cj^ with k = (—3, —2)^ is the segment on the imaginary axis shown in Figure 
5, where b = — ^T. The essential spectrum (= continuous spectrum) of the linear 
2D Euler operator at this fixed point is the entire imaginary axis. Denote by L the 



S>{A} 



i2\b\ 



3?{A} 



-i2\b\ 



Figure 5. The spectrum of £^ with k = (-3, -2)^, when p = (1, 1)^. 

right hand side of (10.8), i.e. the whole linearized 2D Euler operator. Let C denote 
the number of points q £ I? /{Q} that belong to the open disk of radius \p\, and 
such that q is not parallel to p. 

Theorem 10.11 ([112]). The number of nonimaginary eigenvalues of L (count- 
ing the multiplicities) does not exceed 2(. 

Another interesting discussion upon the discrete spectrum can be found in [61]. 



102 10. LINEARIZED 2D EULER EQUATION AT A FIXED POINT 

A rather well-known open problem is proving the existence of unstable, stable, 
and center manifolds. The main difficulty comes from the fact that the nonlinear 
term is non-Lipschitzian. 



CHAPTER 11 



Arnold's Liapunov Stability Theory 



11.1. A Brief Summary 

Let D he a region on tlie (a;,y)-plane bounded by the curves Fj {i = 1,2), an 
ideal fluid flow in D is governed by the 2D Euler equation written in the stream- 
function form: 

(11.1) ^AV' = [VV,VAV'] , 

where 

VAi/)] dAip dip dA-ip 

dx dy dy dx ' 

with the boundary conditions, 

d f dtp 

For every function f{z), the functional 

(11.2) J In ^^^"^^ '^'"^'^ 

is a constant of motion (a Casimir) for (11.1). The conditional extremum of the 
kinetic energy 



W • VV' dxdy 



for fixed F is given by the Lagrange's formula [14] , 

(11.4) dH = 5{E + \F)=Q , ^ Vo = A/'(AV'o) • 

where A is the Lagrange multiplier. Thus, ipo is the stream function of a stationary 
flow, which satisfies 

(11.5) Vo = $(AV'o) , 
where $ = A/'. The second variation is given by [14], 

(11.6) ^2^=1/ / \v(p-V(p + 'i>'{AtPo) {A(Pf\dxdy . 



D 



Let %p = ipQ + ip he a. solution to the 2D Euler equation (11.1), Arnold proved the 
estimates [16]: (a), when c < <I>'(A'(/'o) <C, 0<c<C<oo, 

j j i^^{t)-V^{t)+c[A^{t)f'^dxdy< j j |v^(0)-V^(0)+C[A^(0)]2|da;dy, 
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for all t e (-00, +00), (b). when c < -^'{Atpo) <C,0<c<C<oo, 

J J i^c[Aip{t)]^-V>fi{t)-V>fi{t)^ dxdy< J J i^C[Aip{0)f-V>fi{0}V>fi{0)^dxdy, 

for all t e (—00, +00). Therefore, when the second variation (11.6) is positive 
definite, or when 



V(t>-V(l)+ [max$'(Af/'o)] (A^)^ \ dxdy 



is negative definite, the stationary flow (11.5) is nonlinearly stable (Liapunov sta- 
ble). 

Arnold's Liapunov stability theory had been extensively explored, see e.g. [85] 
[151] [202]. 

11.2. Miscellaneous Remarks 

Establishing Liapunov instability along the line of the above, has not been suc- 
cessful. Some rather technical, with no clear physical meaning as above, argument 
showing nonlinear instability starting from linear instability, has been established 
by Guo et al. [77] [7] [72]. 

Yudovich [209] had been promoting the importance of the so-called "slow col- 
lapse" , that is, not finite time blowup, rather growing to infinity in time. The main 
thought is that if derivatives growing to infinity in time, the function itself should 
gain randomness. Yudovich [208] had been studying bifurcations of fluid flows. 

There are also interests [74] in studying structures of divergent free 2D vector 
flelds on 2-tori. 



CHAPTER 12 



Miscellaneous Topics 

This chapter serves as a guide to other interesting topics. Some topics are 
already well-developed. Others are poorly developed in terms of partial differential 
equations. 

12.1. KAM Theory 

KAM (Kolmogorov-Arnold-Moser) theory in finite dimensions has been a well- 
known topic [11]. It dealt with the persistence of Liouville tori in integrable Hamil- 

tonian systems under Hamiltonian perturbations. A natural idea of constructing 
such tori in perturbed systems is conducting canonical transformations which lead 
to small divisor problem. To overcome such difficulties, Kolmogorov introduced the 
Newton's method to speed up the rate of convergence of the canonical transforma- 
tion series. Under certain non-resonance condition and non-degeneracy condition 
of certain Hessian, Arnold completed the proof of a rather general theorem [11]. 
Arnold proved the case that the Hamiltonian is an analytic function. Moser was 
able to prove the theorem for the case that the Hamiltonian is 333-times differ- 
entiable [162] [163], with the help of Nash implicit function theorem. Another 
related topic is the Arnold theorem on circle map [10]. It answers the question 
when a circle map is equivalent to a rotation. Yoccoz [207] was able to prove an if 
and only if condition for such equivalence, using Brjuno number. 

KAM theory for partial differential equations has also been studied [198] [110] 
[111] [34] [36]. For partial differential equations, KAM theory is studied on a 
case by case base. There is no general theorem. So far the common types of 
equations studied arc nonlinear wave equations as perturbations of certain linear 
wave equations, and solition equations under Hamiltonian perturbations. Often 
the persistent Liouville tori are limited to finite dimensional, sometimes, even one 
dimensional, i.e. periodic solutions [54]. 

12.2. Gibbs Measure 

Gibbs measure is one of the important concepts in thermodynamic and statis- 
tical mechanics. In an effort to understand the statistical mechanics of nonlinear 
wave equations, Gibbs measure was introduced [155] [33] [35], which is built upon 
the Hamiltonians of such systems. The calculation of a Gibbs measure is similar to 
that in quantum field theory. In terms of classical analysis, the Gibbs measure is 
not well-defined. One of the central questions is which space such Gibbs measure 
is supported upon. Bourgain [33] [35] was able to give a brilliant answer. For ex- 
ample, for periodic nonlinear Schrodinger equation, it is supported on H~^^^ [33]. 
For the thermodynamic formalism of infinite dynamical systems, Gibbs measure 
should be an important concept in the future. 
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12.3. Inertial Manifolds and Global Attractors 

Global attractor is a concept for dissipative systems, and inertial manifold is 
a concept for strongly dissipative systems. A global attractor is a set in the phase 
space, that attracts all the big balls to it as time approaches infinity. An inertial 
manifold is an invariant manifold that attracts its neighborhood exponentially. A 
global attractor can be just a point, and often it is just a set. There is no manifold 
structure with it. Based upon the idea of reducing the complex infinite dimensional 
flows, like Navier-Stokes flow, to finite dimensional flows, the concept of inertial 
manifold is introduced. Often inertial manifolds are finite dimensional, have man- 
ifold structures, and most importantly attract their neighborhoods exponentially. 
Therefore, one hopes that the complex infinite dimensional dynamics is slaved by 
the finite dimensional dynamics on the inertial manifolds. More ambitiously, one 
hopes that a finite system of ordinary differential equations can be derived to gov- 
ern the dynamics on the inertial manifold. Under certain spectral gap conditions, 
inertial manifolds can be obtained for many evolution equations [46] [48] . Usually, 
global attractors can be established rather easily. Unfortunately, inertial manifolds 
for either 2D or 3D Navier-Stokes equations have not been established. 

12.4. Zero-Dispersion Limit 

Take the KdV equation as an example 

ut - 6uux + e^Uxxx = , 

Lax [115] asked the question: what happens to the dynamics as e — > 0, i.e. what 
is the zero-dispersion limit ? One can view the KdV equation as a singular pertur- 
bation of the corresponding inviscid Burgers equation, do the solutions of the KdV 
equation converge strongly, or weakly, or not at all to those of the Burgers equation 
? In a series of three papers [115], Lax and Levermore investigated these questions. 
It turns out that in the zero-dispersion limit, fast oscillations are generated instead 
of shocks or multi-valuedness. Certain weak convergences can also established. 

In comparison with the singular perturbation studies of soliton equations in 
previous sections, our interests are focused upon dynamical systems objects like 
invariant manifolds and homoclinic orbits. In the zero-parameter limit, regularity 
of invariant manifolds changes [124]. 

12.5. Zero- Viscosity Limit 

Take the viscous Burgers equation as an example 

Ut + UUx = fiUxx , 

through the Colc-Hopf transformation [88] [44], this equation can be transformed 
into the heat equation which leads to an explicit expression of the solution to the 
Burgers equation. Hopf [88], Cole [44], and Whitham [200] asked the question 
on the zero-viscosity limit. It turns out that the limit of a solution to the viscous 
Burgers equation can form shocks instead of the multi-valuedness of the solution 
to the inviscid Burgers equation 

Ut + UUx = . 



12.7. SLOW COLLAPSE 



107 



That is, strong convergence does not happen. Nevertheless, the location of the 
shock is determined by the multi-valued portion of the solution to the inviscid 
Burgers equation. 

One can also ask the zero-viscosity limit question for Navier-Stokes equations. 
In fact, this question is the core of the studies on fully developed turbulence. In 
most of Kato's papers on fluids [99] [100] [101], he studied the zero- viscosity limits 
of the solutions to Navier-Stokes equations in finite (or small) time interval. It turns 
out that strong convergence can be established for 2D in finite time interval [101], 
and 3D in small time interval [100]. Order yjv rate of convergence was obtained 
[99]. More recently, Constantin and Wu [49] had investigated the zero- viscosity 
limit problem for vortex patches. 

12.6. Finite Time Blowup 

Finite time blowup and its general negative Hamiltonian criterion for nonlin- 
ear Schrodinger equations have been well-known [195]. The criterion was obtained 
from a variance relation found by Zakharov [195]. Extension of such criterion 
to nonlinear Schrodinger equations under periodic boundary condition, is also ob- 
tained [102]. Extension of such criterion to Davey-Stewartson type equations is 
also obtained [73]. 

An explicit finite time blowup solution to the integrable Davey-Stcwartson II 
equation, was obtained by Ozawa [169] by inventing an extra conservation law due 
to a symmetry. This shows that integrability and finite time blowup are compatible. 
The explicit solution q{f) has the property that 

qit)eH%R^) , se(0,i) , te[0,T] , 

for some T > 0. When t T, 

MmH' > C\T - t\-' ^ (x^ , se(0,l). 

The question of global well-posedness of Davey-Stewartson II equation in H^{E?) 
(s > 1) is open. 

Unlike the success in nonlinear wave equations, the search for finite time blowup 
solutions for 3D Euler equations and other equations of fluids has not been suc- 
cessful. The well-known result is the Beale-Kato-Majda necessary condition [21]. 
There are results on non-existence of finite time blowup [52] [53] . 

12.7. Slow Collapse 

Since 1960's, Yudovich [209] had been promoting the idea of slow collapse. 

That is, although there is no finite time blowup, if the function's derivative grows 
to infinity in norm as time approaches infinity, then the function itself should gain 
randomness in space. 

The recent works of Feffcrman and Cordoba [52] [53] are in resonance with the 
slow collapse idea. Indeed, they found the temporal growth can be as fast as e®* 
and beyond. 
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12.8. Burgers Equation 

Burgers equation 

Ut + UUx = HUxx + fit, x) , 
was introduced by J. M. Burgers [37] as a simple model of turbulence. By 1950, E. 
Hopf [88] had conducted serious mathematical study on the Burgers equation 

and found interesting mathematical structures of this equation including the so- 
called Cole-Hopf transformation [88] [44] and Legendre transform. Moreover, Hopf 
investigated the limits fi and t ^ oo. It turns out that the order of taking the 
two limits is important. Especially with the works of Lax [113] [114], this led to 
a huge interests of studies on conservation laws for many years. Since 1992, Sinai 
[189] had led a study on Burgers equation with random data which can be random 
initial data or random forcing [58]. 

Another old model introduced by E. Hopf [87] did not catch too much atten- 
tion. As a turbulence model, the main drawback of Burgers equation is lack of 
incompressibility condition. In fact, there are studies focusing upon only incom- 
pressibility [74]. To remedy this drawback, other models are necessary. 

12.9. Other Model Equations 

A model to describe the hyperbolic structures in a neighborhood of a fixed point 
of 2D Euler equation under periodic boundary condition, was introduced in [133] 
[134] [131] [128]. It is a good model in terms of capturing the linear instability. 
At special value of a parameter, the explicit expression for the hyperbolic structure 
can be calculated [133] [131]. One fixed point's unstable manifold is the stable 
manifold of another fixed point, and vice versa. All are 2D ellipsoidal surfaces. 
Together, they form a lip shape hyperbolic structure. 

Another model is the so-called shell model [97] [181] [96] which model the 
energy transfer in the spectral space to understand for example Kolmogorov spectra. 

12.10. Kolmogorov Spectra and An Old Theory of Hopf 

The famous Kolmogorov —5/3 law of homogeneous isotropic turbulence [106] 
[107] [105] still fascinates a lot of researchers even nowadays. On the contrary, a 
statistical theory of Hopf on turbulence [89] [92] [90] [91] [196] has almost been 
forgotten. By introducing initial probability in a function space, and realizing the 
conservation of probability under the Navier-Stokes flow, Hopf derived a functional 
equation for the characteristic functional of the probability, around 1940, published 
in [89]. Initially, Hopf tried to find some near Gaussian solution [92], then he knew 
the work of Kolmogorov on the —5/3 law which is far from Gaussian probability, 
and verified by experiments. So the topic was not pursued much further. 

12.11. Onsager Conjecture 

In 1949, L. Onsager [167] conjectured that solutions of the incompressible Euler 
equation with Holder continuous velocity of order v > 1/3 conserves the energy, 
but not necessarily if < 1/3. In terms of Besov spaces, if the weak solution of 
Euler equation has certain regularity, it can be proved that energy indeed conserves 
[60] [47]. 
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12.12. Weak Turbulence 

Motivated by a study on 2D Euler equation [210], Zakharov led a study on 
infinte dimensional Hamiltonian systems in the spectral space [211]. Under near 
Gaussian and small amplitude assumptions, Zakharov heuristically gave a closure 
relation for the averaged equation. He also heuristically found some stationary 
solution to the averaged equation which leads to power-type energy spectra all 
of which he called Kolmogorov spectra. One of the mathematical manipulations 
he often used is the canonical transformation, based upon which he classified the 
kinetic Hamiltonian systems into 3-wave and 4-wave resonant systems. 

12.13. Renormalization Idea 

Renormalization group approach has been very successful in proving universal 

property, especially Fcigcnbaimi constants, of one dimensional maps [45]. Rcnormalization- 
group-type idea has also been applied to turbulence [168], although not very suc- 
cessful. 

12.14. Random Forcing 

As mentioned above, there are studies on Burgers equation under random forc- 
ing [58], and studies on the structures of incompressible vector fields [74]. Take 
any steady solution of 2D Euler equation, it defines a Hamiltonian system with the 
stream function being the Hamiltonian. There have been studies on such systems 
under random forcings [63] [62]. 

12.15. Strange Attractors and SBR Invariant Measure 

Roughly speaking, strange attractors are attractors in which dynamics has sen- 
sitive dependence upon initial data. If the Cantor sets proved in previous sections 
are also attractors, then they will be strange attractors. In view of the obvious 
fact that often the attractor may contain many different objects rather than only 
the Cantor set, global strange attractors are the ideal cases. For one dimensional 
logistic-type maps, and two dimensional Henon map in the small parameter range 
such that the 2D map can viewed as a perturbation of a ID logistic map, there have 
been proofs on the existence of strange attractors and SBR (Sinai-Bowen-Ruelle) 
invariant measures [23] [24] [25]. 

12.16. Arnold Diffusions 

The term Arnold diffusion started from the paper [12]. This paper is right 
after Arnold completed the proof on the persistence of KAM tori in [11]. If the 
degree of freedom is 2, the persistent KAM tori are 2 dimensional, and the level 
set of the perturbed Hamiltonian is 3 dimensional, then the 2 dimensional KAM 
tori will isolate the level set, and diffusion is impossible. In [12], Arnold gave 
an explicit example to show that diffusion is possible when the degree of freedom 
is more than 2. Another interesting point is that in [12] Arnold also derived an 
integral expression for a distance measurement, which is the so-called Melnikov 
integral [158]. Arnold's derivation was quite unique. There have been a lot of 
studies upon Arnold diffusions [149], unfortunately, doable examples arc not much 
beyond that of Arnold [12] . Sometimes, the diffusion can be very slow as the famous 
Nekhoroshev's theorem shows [164] [165]. 



110 



12. MISCELLANEOUS TOPICS 



12.17. Averaging Technique 

For systems with fast small oscillations, long time dynamics are governed by 
their averaged systems. There are quite good estimates for long time deviation 
of solutions of the averaged systems from those of the original systems [13] [17]. 
Whether or not the averaging method can be useful in chaos in partial differential 
equations is still to be seen. 
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